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ABSTRACT 


The transition between the light and dark areas of the luminous intensity 
distribution of a headlamp needs to fulfil statutory regulations. These 
regulations restrict the sharpness of the transition. Projection headlamps 
generate a sharp cutoff line. Reducing the sharpness is done by adding 
a structure to the light exit surface of the lens with an amplitude of a 
few micrometers. A long and error-prone process was used to design 
and manufacture structured lenses. The process used so far and the pro- 
blems associated with it are explained. Afterwards a new, faster and less 
error-prone process is presented. An algorithm was developed to struc- 
ture surfaces within the computer-aided design software directly. The 
lenses were manufactured using injection moulding. The moulds were 
manufactured by an ultra precision turning machine. A description of 
the machine and the data interface is given. The process being used until 
now did not include a qualification of the moulds. An optical measuring 
instrument was used to scan the structured surface. The result of the scan 
is a point grid. To verify the structure by a photometric simulation a sur- 
face fitting was done. The algorithm used for fitting is presented. Until 
now, the design of the structure was done manually by a trial-and-error 
process. The time for a photometric simulation was reduced to make the 
design by a genetic algorithm possible. Therefore, the calculation of the 
B-spline basis functions was vectorized. The function of the new process 
is demonstrated with two examples. 

This research and development project was funded by the German Fe- 
deral Ministry of Education and Research (BMBF) within the “Innovati- 
ons for Tomorrow’s Production, Services, and Work” program (funding 
number 02P14A100ff) and implemented by the Project Management 
Agency Karlsruhe (PTKA). 


IH 


KURZFASSUNG 


Der Übergang des hellen in den dunklen Bereich einer automobilen 
Scheinwerfer-Lichtverteilung unterliegt gesetzlichen Auflagen. Die Auf- 
lagen beschränken, wie scharf der Übergang sein darf. Scheinwerfer-Pro- 
jektionssysteme erzeugen zunächst eine zu scharfe Hell-Dunkel-Grenze. 
Um diese aufzuweichen, wird die Lichtaustrittsseite der Linse mit einer 
streuenden Struktur versehen, deren Amplitude im Bereich weniger 
Mikrometer liegt. Die Auslegung und Herstellung der strukturierten 
Linsen erfolgte bisher durch eine lange und fehleranfällige Prozesskette. 
Die bisher verwendete Prozesskette und die damit verbundenen Pro- 
bleme werden kurz erläutert. Anschließend wird eine neue, schnelle 
und fehlerunanfälligere Prozesskette vorgestellt. Es wurde ein Algorith- 
mus entwickelt, um im Computer-Aided Design Linsenflächen direkt 
strukturieren zu können. Die Herstellung der Linsen erfolgte durch 
ein Spritzgussverfahren. Dazu wurden Werkzeugeinsätze benötigt. Die 
zur Herstellung verwendete Ultrapräzisionsdrehmaschine wird kurz 
vorgestellt und die Datenschnittstelle erklärt. In der bisher verwende- 
ten Prozesskette existierte keine Möglichkeit die Werkzeugeinsätze zu 
qualifizieren. Um die lichttechnische Wirkung der eingebrachten Struk- 
tur zu überprüfen, wurde eine optische Messmaschine verwendet. Das 
Ergebnis der Vermessung ist ein Punktgitter der strukturierten Ober- 
fläche. Um eine lichttechnische Simulation durchzuführen, wurde eine 
Fläche ausgehend von dem Punktgitter berechnet. Der dazu verwendete 
Algorithmus wird vorgestellt. Bisher erfolgte die Auslegung der Struktu- 
ren in einem manuellen „Trial and Error" Prozess. Um eine Auslegung 
durch genetische Algorithmen zu ermöglichen, wurde die Zeit für ei- 
ne lichttechnische Simulation reduziert. Dazu wurde die Berechnung 
der B-Spline Basis-Funktionen vektorisiert. Die Funktion der neuen 


KURZFASSUNG 


Prozesskette wird anhand von zwei Beispielen demonstriert. Dieses 
Forschungs- und Entwicklungsprojekt wurde durch das Bundesministe- 
rium für Bildung und Forschung (BMBF) im Programm „Innovationen 
für die Produktion, Dienstleistung und Arbeit von morgen” (Förder- 
kennzeichen 02P14A100ff) gefördert und vom Projektträger Karlsruhe 
(PTKA) betreut. 
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KAPITEL1 


EINLEITUNG 


Seit 2007 die ersten Serienfahrzeuge mit LED-Scheinwerfern vom Band 
liefen, haben diese LED-Scheinwerfer immer mehr Marktanteile gewon- 
nen [1,2]. Durch den Einsatz von LEDs mussten neue Konzepte für 
Scheinwerfer erarbeitet werden [3]. Da der Lichtstrom einer einzelnen 


LED zu gering ist, ist es notwendig, die Straße mit mehreren LEDs zu 
beleuchten. Dazu werden häufig Projektionssysteme verwendet (vgl. 
Abb.|1.1). Ein Projektionssystem sammelt über einen Reflektor oder ei- 
ne Primäroptik das Licht der LEDs und konzentriert es in der Ebene 


einer Blende [4]. Die Blendenebene wird durch eine Linse auf die Stra- 
fe abgebildet. Die Verwendung von Freiformlinsen ergeben Freiheiten 
bei Design und geometrischer Auslegung, wie z. B. die Reduzierung 
der Baugröße oder die Abbildung mehrerer LEDs durch eine Linse [5]. 
Durch die Kontur der Blende generieren Projektionssysteme eine sehr 
scharfe Hell-Dunkel-Grenze (HDG) in der Lichtverteilung, wie in Abb. 


1.2|dargestellt [6]. Eine scharfe HDG reduziert den visuellen Komfort 


fiir den Fahrer. Um Objekte zu erkennen, sucht der Mensch nach Kanten 
bzw. Kontrasten [7]. Eine scharfe HDG führt dazu, dass der Fahrer sei- 
nen Blick ständig auf die HDG fokussiert. Um die HDG aufzuweichen, 
wird die Lichtaustrittsseite der Linse mit einer licht-streuenden Struktur 


überlagert [6]. Da die Amplitude der Strukturen im Bereich weniger Mi- 
krometer ist, werden diese als Mikrostrukturen bezeichnet. Eine weiche 


4 


HDG ist in Abb. 1.2\dargestellt. Die „Economic Commission for Europe’ 


(ECE) stellt gesetzliche Vorgaben, damit ein Scheinwerfer zugelassen 
wird EI Unter anderem gibt die Vorschrift an, wie die Schärfe der HDG 
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Abbildung 1.1: BiLED-Scheinwerfermodul 


Abbildung 1.2: Vergleich zwischen einer scharfen und einer weichen HDG. 


zu bestimmen ist und definiert den legalen Bereich für den Scharfegrad. 
Bei der Auslegung von Mikrostrukturen muss einerseits die Schärfe 
der HDG reduziert, andererseits darf durch das gestreute Licht der Ge- 
genverkehr nicht geblendet werden [6]. Daher werden Möglichkeiten 


benötigt, die Mikrostrukturen gezielt auslegen zu können. 


1.1 MOTIVATION UND UMFELD DER ARBEIT 


In dieser Arbeit liegt der Fokus auf Kunststofflinsen. Sowohl die Aus- 
legung als auch die Fertigung mikrostrukturierter Linsen erfolgt bei 
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Auslegung 


NG Freigabe Herstellung a 


Herstellung 


Abbildung 1.3: Bisher verwendete Prozesskette zur Herstellung mikrostrukturierter Lin- 
sen. 


der Firma HELLA GmbH & Co. KGaA durch einen langen und fehleran- 
fälligen Prozess. Bei der Entwicklung eines neuen Scheinwerfers wird 
die Schärfe der HDG zunächst nicht berücksichtigt. Das System wird 
mit einer nicht-strukturierten Linse ausgelegt. Wenn die Vorgaben von 
Kunde und Gesetz erfüllt werden, wird durch Strukturierung der Linse 
die Schärfe der HDG reduziert. Die Prozesskette, die zur Herstellung 


der strukturierten Linsen angewandt wird, ist in Abb. |1.3|dargestellt. 


Die Konstruktion der optischen Flächen, wie die nicht-strukturierte Lin- 
se, erfolgt in dem Computer-Aided Design (CAD) Programm CATIA 
[9]. Die analytisch beschriebene Linsenfläche wird in einer STEP-Datei 


abgespeichert und in dem Programm Rhinoceros (Rhino) geladen [10} 
11]. In Rhino wird auf Basis der Fläche ein Punktgitter berechnet, dass 
die Fläche repräsentiert und in einer ASCII-Textdatei abgespeichert 112]. 
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Die ASCH-Datei mit dem Punktgitter wird in dem Programm Mathe- 
matica eingeladen 113]. Mit Mathematica wird die Struktur hinzugefügt, 
indem die Punkte um wenige Mikrometer verschoben werden. Die ver- 
schobenen Punkte werden in einer ASCII-Datei abgespeichert und in 
das Programm Rhino geladen. Durch eine Fit-Funktion des Programms 
Rhino wird eine Fläche auf Basis des Punktgitters berechnet. Die Berech- 
nung dauert ca. 25 min auf derzeit üblichen Prozessoren. Die berechne- 
te Fläche wird in eine IGES-Datei abgespeichert [14]. Die IGES-Datei 
wird in CATIA geladen. Durch eine vom Unternehmen implementierte 
Schnittstelle wird die Flache in ein HELLA eigenes STEP-Format abge- 
speichert. Die STEP-Datei wird in dem unternehmenseigenen lichttech- 


nischen Simulationsprogramm zusammen mit den restlichen optischen 
Flachen des Scheinwerfermoduls geladen. Das Programm simuliert die 
Lichtverteilung einschließlich des Schärfegrades der HDG. Liegt der 
Schärfegrad außerhalb des Soll-Bereiches oder verursacht die Struk- 
tur zu viel Blendung, so werden erfahrungsgestiitzt die Parameter der 
Struktur in Mathematica angepasst. In manuellen Korrekturschleifen 
wird die Struktur so angepasst, dass die Vorgaben erfüllt werden. Zur 
Fertigung wird ausgehend von der berechneten Fläche in Rhino eine 
Rhino-3ds-Datei exportiert und zur Herstellung von Stahlwerkzeugen 
verwendet. Die Stahlwerkzeuge dienen als Formen für das Spritzgießen 
der Linsen [15]. Die abgeformten Linsen werden in einem gefertigten 
Scheinwerfermodul (vgl. Abb. 1.1) eingebaut. Die Lichtstärkeverteilung 
(LSV) des Moduls wird vermessen. Wenn die HDG nicht in Ordnung ist 
und Simulation und Vermessung voneinander abweichen, wird zunächst 


durch Änderung der Spritzguss-Parameter versucht, eine Verbesserung 
zu erreichen. Kann keine Verbesserung erreicht werden, werden die 
Parameter der Struktur in Mathematica verändert und erneut Stahlwerk- 
zeuge gefertigt. Diese Korrekturschritte werden solange wiederholt, bis 
die HDG der gemessenen LSV in Ordnung ist. 


Durch die vielen manuellen Schritte ist der Prozess langwierig. Die 
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Abbildung 1.4: Die zwei vertikalen und horizontalen Streifen stellen den lokalen Verlust 
der Mikrostruktur durch die Ubertragung von Rhino zu Catia dar. 


vielen Konvertierungen und Dateiformate können für eine fehlerhafte 
Informationsweitergabe sorgen. Bei einer ersten Analyse zeigte sich, 
dass die Mikrostruktur durch die Übertragung von Rhino zu CATIA 
lokal verschwindet, wie in Abb./1.4|dargestellt. Die Verwendung von Ma- 
thematica erfordert viel Erfahrung und lässt sich auf Grund der großen 
Freiheiten, die Mathematica bietet, schlecht standardisieren. Durch die 
Verwendung des kommerziellen Programms Rhino ist eine „Black Box” 
in der Prozesskette vorhanden, dessen Algorithmen nicht bekannt sind. 
Es existiert keine Möglichkeit zu überprüfen, ob die Struktur erfolgreich 
in das Stahlwerkzeug eingebracht wird. Dadurch entstehen teure Kor- 
rekturschleifen, in denen wiederholt Stahlwerkzeuge gefertigt werden 
müssen. Für das BiLED-Modul (vgl. Abb. wurden insgesamt acht 
Strukturen ausgelegt und für jede Struktur Stahlwerkzeuge gefertigt, 
bis die HDG des Moduls in Ordnung war. Durch die hohe Anzahl an 
Korrekturen kann die Laufzeit der Projekte überschritten werden und es 
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werden Ressourcen verschwendet. 

Im Rahmen des BMBF-Förderprojektes „effiziente Photonikproduktion 
durch intelligente Technologie” (ePiTec) ist das Ziel dieser Arbeit die 
Verbesserung der Prozesskette [16,17]. Die neue Prozesskette soll die 
Anzahl an Korrekturschleifen verringern bzw. idealerweise gänzlich 


verhindern. Dazu sind die folgenden Ziele zu erreichen: 


1. Die Auslegung von Mikrostrukturen soll direkt in dem CAD- 
Programm CATIA oder in dem lichttechnischen Simulationspro- 
gramm ermöglicht werden. Durch einen Projektpartner wird die 
Fertigungstechnik zur Herstellung der Stahlwerkzeuge verbessert. 
In Zusammenarbeit mit dem Projektpartner soll eine neue Daten- 
schnittstelle zwischen Auslegung und Herstellung implementiert 
werden. 


2. Ein weiterer Projektpartner ermöglicht die optische Vermessung 
der Stahlwerkzeuge. Das Ergebnis der Messung ist ein Punktgitter. 
Um zu überprüfen, ob das Werkzeug in Ordnung ist, soll eine 
lichttechnische Simulation auf Basis des Punktgitters ermöglicht 
werden. 


3. Die Auslegung der Mikrostrukturen soll durch den Einsatz von 
Optimierungsalgorithmen automatisiert werden. 


1.2 STAND DER TECHNIK 


Bezüglich der Auslegung von mikrostrukturierten Linsen zum Aufwei- 
chen der HDG existieren kaum wissenschaftliche Veröffentlichungen, 
jedoch einige Patente. HOLTZ und SCHMIDT beschreiben ein mathe- 
matisches Modell, dass einen Zusammenhang zwischen der Rauigkeit 
der Linsenoberfläche und der Schärfe der HDG bzw. eines Blendwer- 
tes enthält [18]. Es werden Herstellungsverfahren genannt, mit denen 
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die gewünschte Rauigkeit erreicht werden kann. FISCHER stellt ein Ver- 
fahren vor, bei dem die Schärfe der HDG durch punktförmige Defekte 
innerhalb der Linse reduziert wird [19]. Die Defekte werden durch einen 
Laser eingebracht. Um die Position der Defekte zu bestimmen, wird die 


Veränderung der HDG nach dem Einbringen eines Defektes vermessen. 
Mit den Informationen wird ein neuronales Netz trainiert, um die op- 
timalen Positionen der Defekte bestimmen zu können [20]. HAMKENS 
berichtet über die Herstellung von Linsen durch eine Kombination aus 


Spritzguss- und Press-Verfahren und erwähnt in diesem Zusammenhang 
Strukturen, die der Oberfläche eines Golfballs ähnlich sind [21]. BONITZ 
ET AL. überlagern zwei sinusförmige Strukturen mit Amplituden im 
Bereich 0,1 um bis 10 pm [22]. Die Strukturen sind so ausgerichtet, dass 
die Wellenbäuche und -täler konzentrische Kreise um die optische Achse 


bilden. KIESEL verwendet regelmäßige Strukturen der Größe 3 um bis 
30 um und Strukturen, die durch Sand- oder Kugelstrahlen hervorgeru- 
fen werden [23 . Von KIESEL ET AL. werden Strukturen erläutert, deren 
Amplitude z durch Kosinus-Funktionen der Form z(u) = A: cos* (u) 


gegeben sind [24]. Die Streuwirkung wird durch die Parameter A und e 


sowie Überlagerung mehrerer Funktionen eingestellt. Mit diesen Funk- 
tionen soll eine exakte Einstellung des Schärfegrades möglich sein. 


Wie für Patente üblich, werden die beanspruchten Produkteigenschaften 
erläutert. Es wird nicht erklärt, wie die Parameter der Eigenschaften für 
ein konkretes Produkt bestimmt werden. 


1.3 STRUKTUR DER ARBEIT 


In der Arbeit wird eine neue Prozesskette vorgestellt. Die Kapitel be- 
schäftigen sich mit den einzelnen Schritten zur Herstellung mikrostruk- 
turierter Linsen. In jedem Kapitel wird das notwendige theoretische 
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Hintergrundwissen erläutert und die Ergebnisse bzw. umgesetzten Ver- 
besserungen dargestellt. 


In Kapitel 2 werden die mathematischen Grundlagen beschrieben, mit 


denen im CAD Flächen dargestellt werden. Es werden Algorithmen ent- 
wickelt, um die Flächen auf modernen Prozessoren schnell auswerten 
zu können. Die Entwicklung von Scheinwerfern umfasst lichttechnische 
Simulationen. Die verwendete Methodik zur Simulation wird erläutert. 


Das Kapitel |3|beschreibt die neue Prozesskette. Es wird erklärt, wie 


eine glatte Fläche strukturiert wird und mit welchen Parametern die 
Strukturen definiert werden. Simulationsergebnisse, die die Streuwir- 
kung der Strukturen demonstrieren, werden präsentiert. Die verwendete 
Fertigungs- und Messtechnik wird beschrieben. Bei der Rückführung 
von Punktegittern in analytisch beschriebene Flächen existieren viele 
Möglichkeiten. Zunächst wird eine Übersicht über verschiedene An- 
sätze gegeben und anschließend das Verfahren im Detail erläutert, das 
in dieser Arbeit verwendet wird. In Kapitel 4|werden die Anforderun- 


gen an eine Scheinwerfer-Lichtverteilung, die für diese Arbeit relevant 
sind, dargestellt. An zwei Beispielen wird die Funktionsweise der neuen 
Prozesskette demonstriert. Um die Auslegung der Strukturen zu automa- 


tisieren wird in Kapitel|5|\die Anwendung von genetischen Algorithmen 


untersucht. In Kapitel|6|werden die Ergebnisse zusammengefasst und 


weitere Verbesserungsmöglichkeiten aufgezeigt. 


KAPITEL 2 


RECHNERGESTÜTZTE AUSLEGUNG 
OPTISCHER FLÄCHEN 


Bei der Entwicklung von automobilen Scheinwerfern ist die Erstellung ei- 
nes virtuellen Prototypen unabdingbar [25} 26]. Um Entwicklungskosten 
und -zeiten zu reduzieren, muss ein Computermodell des Scheinwer- 
fers erstellt, eine lichttechnische Simulation durchgeführt und bewertet 


werden, bevor ein erster physischer Prototyp gefertigt wird. Dazu wird 
eine Möglichkeit benötigt, die Geometrie des Scheinwerfers abzubil- 
den, diese mit optischen Eigenschaften zu versehen und zu simulieren. 
Die in dieser Arbeit verwendeten Modelle und Methodiken werden im 


Folgenden erklärt. 


2.1 RAY-TRACING 


Zur Simulation von Scheinwerfer-Lichtverteilungen wird das sogenann- 
te vorwärts Ray-Tracing verwendet [27]. Das Licht wird in Form von 


Strahlen modelliert und ausgehend von der Lichtquelle durch den op- 
tischen Aufbau verfolgt, bis die Strahlen eine zu beleuchtende Fläche 
erreichen oder den Simulationsaufbau verlassen [28]. In Abbildung/2.1 
sind schematisch verschiedene Vorgänge des Ray-Tracing dargestellt. 
Abhängig von dem Material der Fläche, auf die ein Strahl trifft, finden 


unterschiedliche Beeinflussungen statt. Ein Strahl kann unter anderem 
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LED Blende Linse 
Abbildung 2.1: Schematische Darstellung des Ray-Tracing. 


spiegelnd reflektiert (A), absorbiert (B), fresnel reflektiert (C), diffus 
reflektiert (D) oder gebrochen (E) werden [29]. 


2.1.1 MONTE CARLO-METHODE 


Die Interaktionen zwischen Strahlen und Flächen lassen sich in zwei Ka- 
tegorien einteilen. Eine Kategorie umfasst die Interaktionen, bei denen 
es nur einen Ergebnisstrahl gibt. Dazu zählen die spiegelnde Reflexi- 
on, die Absorption und die Brechung. Die zweite Kategorie umfasst 
die Interaktionen, bei denen mehrere Ergebnisstrahlen möglich sind. 
Bei Fresnel-Relexionen wird in Abhängigkeit der Polarisation und des 
Auftreffwinkels ein Teil des Lichtes reflektiert und der restliche Teil in 
das Medium hineingebrochen. Bei diffusen Reflexionen (z. B. nach dem 
Lambertschen Gesetz [30]) wird das Licht in einen Halbraum gestreut. 
Der in dieser Arbeit verwendete Ray-Tracer der Firma HELLA nutzt die 


Monte Carlo-Methode, um Strahl-Flächen-Interaktionen mit mehreren 
Ausgangsstrahlen zu handhaben [31,32]. 
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Trifft ein Strahl auf eine brechende Fläche, so werden entsprechend 
der Fresnelschen Formeln der Reflexionsgrad und der Transmissi- 
onsgrad bestimmt. Der Strahl wird mit der Wahrscheinlichkeit, die 
dem Reflexionsgrad entspricht, reflektiert und mit der Wahrscheinlich- 
keit, die dem Transmissionsgrad entspricht, gebrochen. Bei diffusen 
Reflexionen, die einem Lambert-Strahler entsprechen, wird die aus- 
gehende Richtung des Strahls aus einer zweidimensionalen Kosinus- 
Wahrscheinlichkeitsverteilung gezogen. 


Ein Nachteil der Monte Carlo-Methode ist das entstehende Simulati- 
onsrauschen [33]. Durch die auf Wahrscheinlichkeitsverteilungen basie- 


renden Strahl-Flächen-Interaktionen können einige Detektorzellen zu 
viel oder zu wenig Licht enthalten. Die Größe des Simulationsrauschens 
ist näherungsweise antiproportional zur Wurzel der Strahlanzahl. Um 
das Simulationsrauschen zu halbieren, muss die bei der Simulation ver- 
wendete Strahlanzahl vervierfacht werden. Besonders problematisch 
ist das Rauschen in Bereichen geringer Lichtströme der simulierten 
Lichtstärkeverteilungen. Diese Bereiche werden im Vergleich zu heller 
ausgeleuchteten Bereichen von nur wenigen Strahlen getroffen, sodass 
das Rauschen besonders stark ausfällt. Daher müssen zur Simulation 
einer auswertbaren HDG, abhängig von dem konkreten Aufbau des 
Scheinwerfermoduls, mehrere hundert Millionen Strahlen simuliert wer- 
den. 


2.1.2 SCHNITTPUNKTSBERECHNUNG 


Bei dem Ray-Tracing müssen wiederholt die Schnittpunkte zwischen 
einem Strahl und einer Fläche gefunden werden. Ein Strahl ist gegeben 


durch [28] 


Rt)=Ro+Rat, (2.1) 
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wobei Ry den Startpunkt und Ra den Richtungsvektor angibt. Um den 
Schnittpunkt mit einer Fäche Sin, v) zu berechnen, muss die Gleichung 

S(u,v) = R(t) (2.2) 
S(u,v) = Ro +Rg-t (2.3) 

gelöst werden. Dazu wird die Funktion 
F(u,v,t) = S(u,v) — Ro— Rat (2.4) 


definiert. Um die Nullstelle von F zu bestimmen, wird das Newton- 
Raphson-Verfahren verwendet [34]. Das Verfahren berechnet ausgehend 
von einem Startwert ho = (uo, Vo, to) iterativ die Nullstelle. Dazu wird 
die Jacobi-Matrix 


S 
J(u,0) = 2 Sy (u,v) 35Sy(u,0) Ray (2.5) 
En S 


aufgestellt. Durch Lösen des linearen Gleichungssystems 
J(h) - Ah = -F(h) (2.6) 


werden die Korrekturen Ah bestimmt. Der neue Lösungsvektor ist gege- 
ben durch [34]: 


neu mg Hatt + Ah (2.7) 


Das Verfahren wird iteriert, bis die Konvergenz erreicht ist. 


Die Suche nach Schnittpunkten durch das Newton-Raphson-Verfahren 
kann durch den Einsatz einer Hierarchie von Begrenzungsvolumen be- 
schleunigt werden [33]. In Abbildung 2.2]ist eine Hierarchie dargestellt. 
Die Fläche wird in kleine Teilflächen unterteilt und für jede Teilfläche 


ein Begrenzungsvolumen berechnet. Zum Finden der Schnittpunkte 
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Abbildung 2.2: Schematische Darstellung einer Hierarchie von Begrenzungsrahmen (Be- 
grenzungsvolumen) einer Kurve (blau). Die vom Strahl getroffenen Begrenzungsrahmen 
sind grün markiert. 


wird die Hierarchie der Begrenzungsvolumen von außen nach innen 
durchquert und nach Schnittpunkten zwischen Strahl und Begrenzungs- 
volumen gesucht. Wenn ein innen liegendes Begrenzungsvolumen vom 
Strahl geschnitten wird, wird nach dem Schnittpunkt auf der Fläche 
selbst gesucht. Durch die Verwendung von Begrenzungsvolumen kann 
die Verwendung von numerisch aufwendigen Verfahren reduziert wer- 
den. Die Berechnung von Schnittpunkten mit Ebenen ist schneller als die 
Suche nach Schnittpunkten auf Freiformflächen mit, z. B., dem Newton- 
Raphson-Verfahren. 


2.2 FLÄCHENDARSTELLUNG 


In der Automobilindustrie sind Freiformflächen nicht mehr wegzuden- 
ken. In den 1960er Jahren entwickelte Pierre Bezier die sogenannten 
Bezier-Kurven und -Flächen, um Freiformflächen modellieren zu kön- 
nen 36]. Nachteil der Bezier-Kurven und -Flächen ist der globale 
Einfluss der Kurven- bzw. Flächenparameter. Die Änderung eines Pa- 
rameters ändert den Verlauf der gesamten Geometrie [37]. Um lokale 
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Änderungen an der Geometrie vornehmen zu können, ist eine Flächen- 
beschreibung durch „Basis Splines” (B-Splines) besser geeignet [3839]. 


2.2.1 DAS KOORDINATENSYSTEM 


Das in dieser Arbeit verwendete Koordinatensystem ist in Abb. 2.3 


dargestellt. Die z-Achse zeigt aus Sicht eines Autofahrers in die Vor- 


Abbildung 2.3: Ausrichtung des Koordinatensystems. 


wärtsrichtung. Die y-Achse zeigt vertikal nach oben. Die x-Achse bildet 
mit den anderen Achsen ein rechtshändiges Koordinatensystem und 
zeigt damit aus Fahrersicht von rechts nach links. In einigen Kapiteln 
wird die Lichtrichtung in Winkeln angegeben. Beide Winkel werden 
bezüglich der z-Achse angegeben. Der Winkel « verläuft horizontal. Po- 
sitive Werte zeigen nach rechts. Der Winkel ß verläuft vertikal. Positive 
Werte zeigen nach oben. 


2.2.2 BASIS SPLINES 


Die Notation zur Beschreibung der Splines ist entnommen. Sei 
U = {ug,... ‚ur} der sogenannte Knotenvektor, bestehend aus einer 
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monoton steigenden Sequenz reeller Zahlen u;. Die u; werden als Knoten 
bezeichnet. Sind die Abstände zwischen den u; nicht äquidistant, so 
wird der Knotenvektor als nicht-uniform bezeichnet. Die B-Spline Basis- 
Funktionen vom Grad p sind definiert durch [40]: 


1 falls u; < u < Hi 
N; u) = a Er 2.8 
iolu) l ansonsten we) 
u—u; Ui+p+1— 4 
N; (u) = ———_N,; „_ı(u) + ————_Nja1 pU 2.9 
ipl ) Wisp — Ui i,p 1( ) Wit ps1 — Hin i+1,p 1( ) ( ) 


Eine B-Spline-Kurve basiert auf n + 1 Kontrollpunkte P; € R und ist 
gegeben durch [40]: 


n 
u)=),N;,(u)P; mitu € [uo ur] (2.10) 
i=0 


Zwischen der Anzahl der Knoten r + 1 und der Anzahl der Kontroll- 
punkte n + 1 gilt die Beziehung r = n + p + 1. Zur Veranschaulichung 
von Gl. sind in Abb. 2.4| zwei B-Spline-Kurven dargestellt. Die 
Kontrollpunkte P; spannen ein Polygon auf, dass die B-Spline-Kurve 


einhüllt. Abhängig von den Knoten wird der Parameterbereich defi- 
niert, in dem die Basis-Funktionen ungleich null sind (siehe Gl. (2.9)) 
und dadurch bestimmt, wie stark sich die einzelnen Kontrollpunkte auf 
den Verlauf der Kurve auswirken. Ist ein Knoten p-mal vorhanden, so 
entstehen Knicke in der Kurve (siehe blaue Kurve in Abb.]2.4). 


Durch Verwendung eines zweiten Knotenvektors V = {vo,... ‚vs} und 
einer Kontrollpunktmatrix Pj € RÊ wird eine B-Spline-Fläche definiert 
durch [40]: 


n m 


= LL Nina) Po (2.11) 


i=0 j=0 
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e Pi --- Kontrollpolygon == C(u) 4 | Knoten 


Abbildung 2.4: Vergleich von zwei B-Spline-Kurven mit den gleichen Kontrollpunkten 
und unterschiedlichen Knotenvektoren. Die äußeren Knoten sind vierfach in dem Knoten- 
vektor vorhanden. 


Hierbei gibt p den Grad in u-Richtung und q den Grad in v-Richtung der 
Fläche an. In Abb.|2.5]ist eine B-Spline-Fläche dargestellt. Im Gegensatz 
zur B-Spline-Kurve hat eine B-Spline-Fläche eine Kontrollpunktmatrix 
anstatt eines Vektors. Beide Geometriedarstellungen (Kurve und Fläche) 


basieren jedoch auf Knotenvektoren. Das hat zur Folge, dass z. B. eine 
Änderung eines Knotens in dem u-Knotenvektor einer Fläche die Fläche 
entlang der kompletten v-Richtung ändert. 


Die Ableitung in u- und v-Richtung einer B-Spline-Fläche berechnen 
sich aus [40] 


EES (2.12) 
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wobei die Ableitungen der Basis-Funktionen bestimmt werden durch 


[40]: 


(k-1) (k-1) 
N® = p- Na Nii p-1 


Uitp— Uji Win" Yi41 


(2.13) 


e Pij —— Kontrollpolygon Eu S(u,v) 


y(u, v) 
y(u, v) 


0 1 2 3 4 5 0 1 2 3 4 5 
x(u, v) x(u, v) 


Abbildung 2.5: Darstellung einer B-Spline-Fläche (rechts) und die dazugehörigen Kon- 


trollpunkte und das Kontrollpolygon (links). Die Knotenvektoren in u- und v-Richtung 
sind [0, 0,0,0,%,%1,1,1,1]. 
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Ni+1,3 | Ni+2,3 | Ni+2,3 


Ni-4,2 Ni+1,2 | Ni+2,2 | Ni+2,2 


Nj-3,1 


Ni-4,0 | Ni-3,0 


Ni-41 Nun | Ni+2,1 | Ni+2,1 


Ni+1,0 | Ni+2,0 | Ni+2,0 


Abbildung 2.6: Grafische Darstellung der Abhängigkeit gegeben durch Gl. und Gl. 
für kubische (p = 3) Basis-Funktionen. 


2.2.3 KONVENTIONELLE BERECHNUNG DER 
BASIS-FUNKTIONEN 


Konventionelle Berechnung der Basis-Funktionen Zur effizienten Berech- 
nung einer B-Spline-Kurve eignet sich die direkte Auswertung von Gl. 
(2.10) nicht, da für ein bestimmtes u ein Großteil der Basis-Funktionen 
null ist. Nach Gl. und existieren höchstens (p +1) Basis- 
Funktionen, die ungleich null sind. Für ein u € [u;,u;,1[ sind das die 
Basis-Funktionen N;—p,p bis N; p. Diese Eigenschaft ist in Abb. 2.6]für ku- 
bische Basis-Funktionen (p = 3) veranschaulicht. Die grünen und blau- 


en Pfeile deuten die rekursive Abhängigkeit der einzelnen Summanden 
von Gl. an. Die Basis-Funktionen außerhalb der roten Umrandung 
sind null, dass heißt, es müssen nur die Basis-Funktionen innerhalb 
berechnet werden. Bei der Bestimmung aller nicht-verschwindenden 
Basis-Funktionen ergeben sich durch Gl. redundante Berechnungen. 
Um die redundanten Berechnungen zu vermeiden werden die Notatio- 


nen 


left[j] =U- Ui+1-j und right[j] = Wii — u (2.14) 


eingeführt. Damit wird der Algorithmus 2.1 definiert, der in einigen 
48]. Zur Be- 
rechnung der Basis-Funktionen einschließlich deren Ableitungen kann 


quelloffenen Programmbibliotheken verwendet wird [42 
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der Algorithmus /2.2| verwendet werden. Dieser vermeidet ebenfalls 


redundante Berechnungen und wird von den quelloffenen Programmbi- 
bliotheken [44448] verwendet. 


Programmauflistung 2.1: Berechnung der nicht-verschwindenden Basis-Funktionen 


1 procedure BasisFuns (i,u,p, U) 


2 Nple1 

3 nor j=l 

4 left[j] — u- Uļli +1- j] 

5 right[j] + Ufi + j] — u 

6 saved + 0 

7 for r+0,j-1 

8 temp + N[r] / (right[r + 1] + left[j — r]) 
9 N[r] < saved + right[r + 1] - temp 

0 saved + left[j — r] - temp 

1 end for 


2 NIT < saved 
3 end for 
4 return N 


5 end procedure 


Programmauflistung 2.2: Berechnung der nicht-verschwindenden Basis-Funktionen und 
deren n-fachen Ableitungen [40] 


1 procedure DersBasisFuns (i,u,p,n,U) 


2 ndu[0][0) <1 

S ror j= ibe 

4 left[j] — u- Uļli +1- j] 
5 right[j] + Ufi + j] — u 

6 saved + 0 
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Kor 
ndufj][r] < right[r + 1] + left[j — r] 
temp + ndu[r][j — 1] / ndufj] [r] 
ndulr][j] +- saved + right[r + 1] - temp 
saved + left[j — r] - temp 
end for 
ndulj][j] + saved 
end for 
for j+ 0,p // Load basis functions 
ders[0][j] = ndufj][p] 
end for 
for r+-0,p 
al = © 
= Il 
a[0][0] — 1 
for k¢1,n // Loop to compute nth derivative 
d 4&0 
ilk = p= Is 
pk=-p-k 
if r>k 
a[s2][0] <- a[s1][0] / ndu[pk + 1] [rk] 
d < als2][0] - ndu[rk] [pk] 


end if 
if rk > —1 
lei 
else 
jil s= -rk 
end if 
if r—1< pk 
zf) 
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else 
[B= fr 
end if 
for j< j1,j2 
als2][j] < (a[s1][j] — a[s1][j — 1]) / ndu[pk + 1][rk + j] 
d + d+als2]|j] -ndu[rk + j] [pk] 
end for 
if r<pk 
a|s2][k] + —a[s1][k — 1] / ndu[pk + 1][r] 
d < d +a[s2]{k] - ndulr] [pk] 
end if 
ders[k][r] +- d 
i=l Gil <= nl 
end for 
end for 
en 
for k+ k,n 
for j+ 0,p 
ders[k][j] + r - ders[k] [j] 
end for 
la) 
end for 
return ders 


60 end procedure 


2.2.4 BERECHNUNG DURCH AVX2-INSTRUKTIONEN 


Um eine schnelle Simulation der lichttechnischen Auswirkung von opti- 


schen Flächen zu ermöglichen, müssen die B-Spline-Flächen und deren 
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erste Ableitungen schnell berechnet werden (vgl. Kap. 2.1). Die Algorith- 
men|2.1}und [2.2|/berechnen ausschließlich die nicht-verschwindenden 
Basis-Funktionen (roter Bereich in Abb. |2.6), jedoch sequentiell. Seit 
1966 können bestimmte Prozessoren die gleichen Operationen auf 


mehreren Daten durch sogenannte „Single Instruction, Multiple Data”- 
Anweisungen (kurz SIMD) gleichzeitig ausführen [49]. Durch Ein- 
führung von „Advanced Vector Extensions’-Instruktionen können 
Operationen auf vier Gleitkommazahlen mit doppelter Genauigkeit 
angewendet werden 1501. Diese Anweisungen können die Berechnung 
der Basis-Funktionen beschleunigen. Anstatt die Basis-Funktionen 
sequentiell zu berechnen, können mit SIMD-Instruktionen die Basis- 


Funktionen innerhalb einer Zeile aus Abb 2.6! gleichzeitig berechnet 


werden. Dadurch werden auch Basis-Funktionen berechnet, die null 
sind (N;_32 N;_31, Ni-21), wobei sich aufgrund der Gleichzeitigkeit 
dennoch eine Beschleunigung ergibt. 


In dieser Arbeit liegt der Fokus auf kubische B-Spline-Flächen. Kubi- 
sche Flächen sind kriimmungsstetig. Damit lassen sich tangentenstetige 
Lichtstärkeverteilungen erzeugen. Es wird ein Algorithmus entwickelt, 
um Gl. und durch AVX-Instruktionen berechnen zu können. 
Dazu werden die Notationen [51] 


1 
Ka eee o 2.15 
CLip itp — u; ( ) 

1 
Din (2.16) 

= Ui+p+1 Hoi 

C3ip = Ui CA (2.17) 
Chip = —Uitp+ ` Cip (2.18) 
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eingeführt und Gl. umgeformt zu: 
Ni p(t) = (Crip: U +349) Nip-1t (&2ip u + csip) Mis (2.19) 


Das Symbol © gibt die elementweise Multiplikation an. Die Berechnung 
der Basis-Funktionen N;_3,, bis Nj, von Grad p < 3 ist gegeben durch: 


Nip Clip u C3,i,p Nip-1 
Ni-1p| _ | |¢1i-1p u C3,i-1,p N;-1,p-1 
© + © 
Ni-2,p C1i-2,p u C3,i-2,p Ni-2,p-1 
Nj-3,p C1i-3,p u C3,i-3,p Nj-3,p-1 
N amaaa 
simd_fmadd 
EE 
simd_mul 
C2 ip u Chip Ni41,p-1 
Ginet el. le Nip-1 
=. Co j— u C4 j— Nj-1 p— 
simd add 2,i—2,p 4,i—2,p i—1,p—1 
C2,i-3,p u C4,i-3,p Ni-2,p-1 
simd_fmadd 
BEE 
simd_mul 


(2.20) 


Durch Gl. (2.20) wird eine Zeile aus Abb 2 ëbasierend auf der darunter 
liegenden Zeile berechnet, wobei die Berechnung der einzelnen Elemen- 


te jeder Zeile durch SIMD-Instruktionen stattfindet. Die Anweisung 
simd_add steht für die gleichzeitige Addition, die Anweisung simd_mul 
für die gleichzeitige elementweise Multiplikation von zwei Vektoren 
mit je vier Gleitkommazahlen [50, 52]. Die Anweisung simd_fmadd (,,fu- 
sed multiply-add“) bezeichnet eine vereinte Multiplikation-Addition 


[53]. Diese hat den Vorteil, dass, im Gegensatz zur sequentiellen Mul- 


tiplikation und anschließenden Addition, nur eine Rundung stattfin- 
det und die Berechnung schneller ist [54]. Algorithmus|2.3| verwendet 
SIMD-Anweisungen zur Berechnung der nicht verschwindenden Basis- 


Funktionen vom Grad drei. 
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Programmauflistung 2.3: Berechnung der nicht-verschwindenden kubischen Basis-Funk- 
tionen durch SIMD-Anweisungen 


1 procedure getBasisFuns (i, u, C11, C21, €31, C41, 


2 C12, C22, C32, C42, C13, C23, C33, C43) 
3 s+-i-3 

A Uayx + {u, u, u, u} 

5 

6 // Contains IN; 30, N; 2,07 N; 1,0, N;o} 

Z Nau = {0, 0, 0, 1} 

8 

9 / Compute {N; 31, NS an, Niii Nis} 

10 Nau +getNextDegreeBasisFuns(s, c11, C21, C31, 

11 C41, Uavx, Navx) 
12 // Compute {Nj_39, Ni-22, Ni-12, Niaz} 

133 Nayx <getNextDegreeBasisFuns(s, c12, C22, C32, 

14 C42, Uavx, Navx) 
is // Compute {N 33, N; 23, N; 13, Nia} 

16 Navx + getNextDegreeBasisFuns(ßs, cı3, C23, C33, 

17 c43, Uavx, Navx) 


18 return Nayx 

19 end procedure 

20 procedure getNextDegreeBasisFuns(s, c1, C2, C3 

21 C4, Uavx, Navx) 

=æ 7/ Compute IN oy, Nite Nig Ni ax by eyelie 
permutation 

23 Na AVX di: permuteLeft(N AVX) 

24 // Compute {N; 2%, Nr Nik, UN 

25 Nzayx + setLastElementZero(N> avx) 


27 // Load {clsl, cıls+1], cıls+2], cils +3]} 
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C1,AVX 7 load AVX(c1, s) 


// Load {cls], c2[ls +1], c2[s +2], co[s +3] } 
C2,avx + loadAVX(c, s) 


// Load {czls], c3[s +1], c3[s + 2], c3[s + 3] } 
C3,AVx 7 loadAVX (cz, s) 


// Load {c4ls], cals+1], cals+ 2], cals + 3]} 
C4, AVX 7 load AVX (ca, s) 


// Compute cy:u+c3 
CLeft  simd_fmadd lu aux, C1,avx, C3,AVX) 
// Compute Co-u+c4 
CRight L simd_fmadd(u AVX, C2,AVX, C4, AVX) 


Navx + simd_add (simd_mul (Giese ‚N AVX) 5 
simd_mul(cright, N2,avx)) 
return Nayx 


47 end procedure 
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Zur Berechnung der Ableitung wird mit Gl. (2.15) bis (2.18) die Gl. (2.13) 


umformuliert zu: 


N 
GP CLip Ni p-1 Crip Ni+1,p-1 
N ipl CLi-1,p Ni-1,p-1 C2,i-1,p Nip-1 
nN) El © N: + Sch © N; 
i—2,p Li—2,p i—2,p—1 2,i—2,p i-1,p—1 
NI C1i-3,p Ni-3,p-1 C2,i-3,p Ni-2,p-1 
i—3,p Ne a en N Seene 
simd_mul simd_mul 
KE 
simd_add 
N, m 
simd_mul 
(2.21) 


In Algorithmus/2.4ist die Berechnung durch SIMD-Anweisungen von 
Gl. (2.21) für kubische Basis-Funktionen dargestellt. 


Programmauflistung 2.4: Berechnung der nicht-verschwindenden kubischen Basis-Funk- 
tionen und deren ersten Ableitungen durch SIMD-Anweisungen 


1 procedure getBasisDerivat(i, u, cq, C21, C31, C41, 


2 C12, C22, C32, C42, C13, C23, C33, C43) 
3 si-—3 

4 Uayx {u, u, u, u} 

5 

6 EE {N; Sir N; 2,0, N; 1,0, Nio} 

7 Nayx = {0, 0, 0, 1} 

8 

o // Compute JN: an, Noi Nin Nia} 

10  Nayx +getNextDegreeBasisFuns(s, c11, C21, C31, C41, 
11 uavx, Navx) 

12 

3 // Compute INS 35, Niaz: Ni-12, Nio} 

14 Nayx + getNextDegreeBasisFuns(s, c12, C22, C32, C42, 
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uavx, Navx) 


// Compute IN: aa Nj-23, Ni-13, Ni3} 
Navx, Héi + getNextDegreeBasisFuns Derivative ( 
S, C13, C23, C33, C43, UAVX, Navx) 


return Nayx, Héi 


end procedure 


procedure getNextDegreeBasisFunsDerivative ( 
S, C1, C2, C3, C4, Uayx, Navx) 
// ‘Compute (N, 5, Nicie Nip Nese} by eyelic 
permutation 
Nz aux 7 permuteLeft(N AVX) 
// Compute {N; au Nat Nix, 0} 
Nz avx + setLastElementZero(N> ayx) 


// Load {c1[s], cils +1], cıls +2], cıls+3]} 
C1,AVX 7 load AVX(c1, s) 


// Load {c2[s], co[s +1], co[s +2], c2[s + 3]} 
C2, AVX 7 load AVX (co, s) 


// Load {czls], c3[s +1], c3[s + 2], c3[s + 3] } 
c3,avx + loadAVX(c3, s) 


// Load {cals], cals +1], cals+ 2], cals + 3]} 
C4, AVX © load AVX (cq, s) 


// Compute cy:u+c3 
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42 CLeft *— simd_fmadd(u AVX, C1,AVX, C3, AVX) 

43 // Compute co-u+cy4 

44 CRight  Simd_fmadd(uAvx, C2,avx, C4,Avx) 

45 

46 // Compute derivative 

47 nw. + simd_add (simd_mul(cı,avx, Navx) , 
48 simd _mul(cz, AVX, N2 avx)) 
a NO) + simd_mul(NU),, {3,3,3,3}) 

50 

a Navx + simd_add(simd_mul(crer, Navx), 
52 simd_mul(cpignt, N2,avx)) 
53 return Nayx, Wu? 

54 end procedure 


Um eine B-Spline-Fläche und die ersten Ableitungen in u- und v- 
Richtung bestimmen zu können, müssen die Gleichungen und 
(2.12) ausgewertet werden. Beide Gleichungen lassen sich als Vektor- 
Matrix-Vektor-Produkt umschreiben [40]: 


S.(u,v) = N(u) - Pe - Ñ (v) (2.22) 


Hierbei steht c € {x,y,z} für eine der drei Raumrichtungen und P- für 
die jeweilige Kontrollpunktmatrix. Die Vektoren sind gegeben durch: 


Ni- -3,3(u) Nj-33(2) 

oun — | N-z) | o= | N-20) 

an Ni-1,3(u) an Nj-13(2) Es 
Ni3(u) Nato) 


Bei der Auswertung einer Fläche muss Gl. (2.22) für jede Raumkoor- 
dinate berechnet werden. Die Umsetzung von Gl. (2.22) durch SIMD- 
Anweisungen ist schwierig, da es keine SIMD-Anweisung gibt, die alle 
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vier Zahlen innerhalb eines AVX-Registers aufsummiert. Diese Art der 
Addition wird bei dem Produkt P; - N(v) benötigt. Nach der Multiplika- 
tion des Vektors mit einer Zeile der Matrix muss die Addition ausgeführt 
werden. Schneller ist die Berechnung, in dem Gl. umgeschrieben 
wird zu: 


S.(u,v) = (Nu) . ÑT (%0), Pe), (2.24) 


Hierbei bezeichnet (, je das Frobenius-Skalarprodukt [55]: 
(A Be = LY anbi (2.25) 
i k 


Nach Gl. (2.24) wird eine Matrix 


NN = N(u) : NI to 
Ni-33°Nj-33 Ni-33°Nj-23 Ni-33°Nj-13 Ni-33°Nj3 
Ni-23°Nj-33 Ni-23°Nj-23 Ni-23°Nj-13 Ni-23°Nj3 
Niwa? Nias  Ni-13: Nos Near rm, Neorg Nj3 
Nig-Nj-33  Ni3-Nj-23 Ni3-Nj-13 Nig: Nijs 
(2.26) 


berechnet. Die Matrix NN wird fiir die Berechnung der Flache in je- 
der Raumrichtung benötigt. Jede Zeile der Matrix wird in einem AVX- 
Register gespeichert. Nach der elementweisen Multiplikation mit jeweils 
einer Kontrollpunktmatrix P. ergibt sich die Matrix: 


NNO Pc= 


NNoo : Pei-3,j-3 NNoi ' Pci-3,j-2  NNo2- Pei-3,j-1 NNos- Pc i-3,j 
NNio- Pei-2,j-3 Nun: Pei-2j-2 NNi2- Pei-2j-1 Na: Pc i-2 
NN: Pei-1,j-3  NNa1  Pei-1,j-2 Ma: Bau NNa: Poi-1j 
NN30°Poij-3  NNai-Poij-2 Ma: Bai  NNa3-Peij 
(2.27) 
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Die Multiplikationen innerhalb einer Zeile werden parallel durch SIMD- 
Anweisungen ausgeführt. Nach der Multiplikation ist jede Zeile aus 
Matrix in einem AVX-Register gespeichert. Nach Gl. müs- 
sen die Elemente aus addiert werden. Dazu werden erst die vier 
Register, die die Matrix enthalten, und anschließend die Einträge 
des verbleibenden Registers aufsummiert. Eine Pseudocode-Darstellung 


des Frobenius-Skalarproduktes ist in Algorithmus 2.5|gegeben. 


Programmauflistung 2.5: Berechnung des Frobenius-Skalarproduktes durch SIMD- 
Anweisungen. 


1 procedure Frobenius (N,avx, Na Aus, Pc) 
2 // Nyavx contains 

{N;-33(u), Ni-23(u), Ni-13(u), Ni3(u)} 
3 // Ny avx contains 


{N;_33(0), Nj-23(v), N;-1,3(0), Nato) 


// Load elements of N,ayx fourfold in an AVX 
register each 
6 Ni3,avx + {Ni-33(u), Ni-33(u), Ni-33(u), Ni-33(u)} 
7  Ni-2avx + {Ni-23(u), Ni-23(u), Ni-23(u), Ni-23(u)} 
8 Ni1avx + {Ni-13(u), Ni-13(u), Ni-13(u), Ni-13(u)} 
a Niavx + {Nis(), Ni3(u), Ni3(u), Na3(u)} 


11 // Mulitply with N,avx, to get NN 
12  N;_3,avx + simd_mul(N;_3,avx, No,avx) 

GD Nj-2,avx + simd_mul(N;_2,avx, No,avx) 

14 Nj-1,avx + simd_mul(N;_1avx, No,avx) 

15 N; avx SE simd_mul(N; avx, No avx) 


7 // Multiply with E 
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18 NP;_3,avx + simd_mul(N;_3,avx, {Pei-3,j-3, Pci-3,j—-2, 


19 een, 
22 NPi-2,avx + simd_mul(Nj-3,avx, {Pei-2,j-3, Pei-2,j-2, 
21 Pos 2, el 
2  NPi-1,avx + simd_mul(N;-1,avx, {Pci-1j-3, Pei-1,j-2, 
23 ee) 


a Ni aux + simd_mul(N;avx, {Poij-3, Poij-2» Peij-1 Peij}) 
25 


26 // Sum up registers 

27 Sumavx <- simd_add (simd_add (NP;_3,avx, NP;_2,avx) D 

28 simd_add (NP;_1,avx, NP; avx)) 

29 = SuMskalar + Sumayx|0] + Sumayx[1] + Sumayx[2] + Sumavx[3] 
30 return SumMstalar 

31 end procedure 


Die Ableitungen der Fläche lassen sich durch Gl. (2.24) berechnen, indem 
einer der beiden Vektoren durch seine Ableitung ausgetauscht wird: 


2 Sc(u,o) = (Zë NI foi, Di (2.28) 
2 S.(u,o) = (New) INT), Be). (2.29) 


Zur Auswertung der Fläche sowie deren Ableitungen müssen für jede 
Raumrichtung die Gleichungen (2.24), (2.28) und (2.29) ausgewertet 


werden. 


2.2.5 VERGLEICH ZWISCHEN KONVENTIONELLER UND 
AVX-BASIERENDER BERECHNUNG 


Zum Testen der erreichten Verbesserung wird ein Computer mit einem 
Intel Xeon E5-1650 v3 verwendet [56]. Um die Laufzeit der konventionel- 
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Tabelle 2.1: Vergleiche der Leistung zwischen konventioneller und SIMD-basierter Aus- 
wertung von B-Spline-Flächen. Die Angaben zu den Implementierungen (Imp.) beziehen 
sich auf den Anhang. 


Berechnung Imp. | Laufzeit (s) | Beschleunigung 
Basis-Funktionen = bg 3,6 
6.2 4,28 
Basis-Funktionen 6.3 42,0 73 
& erste Ableitung 6.4 5,8 ` 
6.5 113,1 
Auswertung der Fläche 56 244 4,6 
Auswertung der Fläche | [6.7 248,4 59 
& erste Ableitungen 6.8 42,0 7 


len mit den AVX-basierten Algorithmen vergleichen zu können, ist eine 
konkrete Implementierung notwendig. Der Quellcode der Implementie- 


rungen ist im Anhang |6|aufgefiihrt. Der Quellcode wird mit Microsoft 


Visual Studio 2012 kompiliert und zum überprüfen der Verbesserung 
108 mal aufgerufen. In Tabelle 2.1lsind die Ergebnisse zusammengefasst. 
Die größte Beschleunigung mit einem Faktor von 7,3 ergab sich bei der 
Berechnung der Basis-Funktionen einschließlich den ersten Ableitungen. 
Bei der Auswertung der Fläche einschließlich Ableitungen bleibt unge- 
fähr eine Verbesserung von einem Faktor sechs übrig. Die Funktionen 
zur Auswertung der Fläche bekommen die u- und v-Parameter überge- 
ben und müssen die dazugehörigen i- und j-Indizes der Knotenvektoren 
bestimmen. Bei den Berechnungen der Basis-Funktionen wurden die 
Indizes vorab bestimmt. Daher ergab sich nur ein Faktor sechs bei der 
Flächen-Auswertung. 


Neben der parallelen Berechnung durch SIMD-Anweisungen basieren 
die Laufzeitunterschiede auf weiteren Punkten. Die Konstanten zur Be- 
rechnung der Basis-Funktionen (Gl. (2.15) bis (2.18)) werden bei der 
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SIMD-Variante vorab, bei der bisher verwendeten Implementierung bei 
jedem Funktionsaufruf neu berechnet. Die SIMD-Variante ist ausschließ- 
lich für kubische Flächen implementiert und kann nur die jeweilige 
Funktion und ihre erste Ableitung berechnen. Die konventionelle Vari- 
ante funktioniert für beliebige Grade und beliebige Ableitungen. Das 
ist möglich, weil vom Grad abhängende Schleifen und Verzweigungen 


verwendet werden (siehe Algorithmen 2.1/und 2.2|sowie Implementie- 


rungen 6.5lund 6.7). 


Einige Prozessor-Instruktionen werden vom Prozessor in Teilaufgaben 
unterteilt. Die Teilaufgaben von aufeinander folgenden Instruktionen 
können gleichzeitig ausgeführt werden, wenn für die Ausführung unter- 
schiedliche Komponenten des Prozessors verwendet werden (sogenann- 
tes „Pipelining“) 1571. Um den Vorteil der Pipeline nutzen zu können, 


müssen zu jedem Zeitpunkt die zukünftigen Instruktionen bekannt sein. 
Enthält das Programm konditionale Sprünge, so versucht der Prozessor 
zu erraten, ob ein Sprung stattfindet (sogenanntes „Branch Prediction”) 


[58]. Wurde falsch geraten, enthält die Pipeline ungültige Instruktionen. 


Die Pipeline muss geleert und mit den korrekten Instruktionen gefüllt 
werden, wodurch ein Zeitverlust entsteht [57]. Durch Vermeidung von 
Verzweigungen kann die Laufzeit eines Programmes reduziert werden 
[59,160]. Daher wurde bei der Entwicklung der SIMD basierten Auswer- 
tung von B-Spline-Flächen auf Schleifen „for“ und Verzweigungen „if“ 


verzichtet, wie der Vergleich von Algorithmus 2.1mit!2.3ljund Algorith- 


mus 2.2mit|2.4 zeigt. 
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PROZESSKETTE ZUR HERSTELLUNG 
MIKROSTRUKTURIERTER LINSEN 


In diesem Kapitel wird erläutert, welche Anforderungen an Schein- 
werfer-Lichtverteilungen bestehen, wie mikrostrukturierte Flächen in 
CAD-Programmen dargestellt, wie die Strukturen in Stahlwerkzeuge 
eingebracht und wie die dazu benötigten Fertigungsdaten aufbereitet 
werden. Es wird die Messapparatur vorgestellt, mit denen die Oberflä- 
che der Stahlwerkzeuge flächig vermessen werden und Algorithmen 
vorgestellt, mit denen Flächen auf Basis der gemessenen Punktwolken 
bestimmt werden können. Der Inhalt dieses Kapitels stammt zu Teilen 
aus der Veröffentlichung [61]. 


3.1 ANFORDERUNG AN LICHTSTÄRKEVERTEILUNGEN 


Es existieren zwei Regelungen, die in den meisten Ländern die Anforde- 
rungen an Scheinwerfern stellen. Eine Regelung wird von der „Economic 
Comission for Europe“ (ECE) ausgegeben, die Zweite ist die „Federal 
Motor Vehicle Safety Standard 108” (FMVSS) [8, |62]. Beide Regelun- 
gen stellen Anforderungen an die LSV eines Scheinwerfers, sodass der 


Scheinwerfer bestimmte Lichtstärke-Werte in vorgegebenen Bereichen 
mindestens haben muss, bzw. nicht überschreiten darf. Es wird in beiden 
Regelungen die gleiche Definition für die Schärfe der HDG, auch ge- 
nannt Gradient oder AK31, verwendet. Die Gradient-Kurve G(ß) wird 
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Abbildung 3.1: Darstellung des Zusammenhangs zwischen der Lichtstärke I, und der 
Gradienten-Kurve G(ß) für eine scharfe (blau) und weiche (orange) HDG. 


aus der Lichtstärke I, bei einem vertikalen Winkel f berechnet EI 62]: 


G(B) = log(I,(B)) — log(v(B + 0,1°)) (3.1) 


Der Maximalwert der Kurve (3.1) wird als „Gradient“ bezeichnet. Der 
Zusammenhang aus Gl. ist in Abb.|3.1|visualisiert. 


Beide Vorschriften geben an, wie die Werte der Lichtstärke zu messen 
sind. Die ECE-Regelung unterscheidet zwischen minimaler und maxi- 
maler Schärfe. Beide Werte werden bei einem horizontalen Winkel von 
& = —2,5° gemessen. Die minimale Schärfe ist durch eine Messung der 
Beleuchtungsstärke in einem Abstand von 10m mit einem Detektor der 
ungefähr einen Durchmesser von 10 mm hat oder in einem Abstand von 
25 m mit einem Detektor der ungefähr einen Durchmesser von 30 mm 
hat, zu bestimmen. Die vorgegebenen Kombinationen aus Abstand und 
Detektorgröße ergeben einen unterschiedlichen Raumwinkelbereich. Da- 
her können die zwei Messkonfigurationen unterschiedliche Ergebnisse 
liefern. Die maximale Schärfe ist durch eine Messung der Beleuchtungs- 
stärke in einem Abstand von 25m mit einem Detektor der ungefähr 
einen Durchmesser von 30 mm hat, zu bestimmen. Beide Messungen 
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sollen mit einer vertikalen Schrittweite von 0,05° durchgeführt werden. 
Die minimale Schärfe muss über einen Wert von 0,13 und die maximale 
Schärfe unter einem Wert von 0,4 liegen. 


Die FMVSS fordert eine Messung bei o = —2,5° oder a = 2,0°. Der De- 
tektor muss einen Durchmesser von 10 mm haben und in einem Abstand 
von 10 m zum Scheinwerfer platziert sein. Der Gradient darf einen Wert 
von 0,13 nicht unterschreiten. Eine obere Grenze ist in der Regelung 
nicht gegeben. 


Die ECE-Regelung definiert unter anderem, dass die Lichtstärke bei 
(x = 0°, B = 0°), bezeichnet als „HV“, mindestens 50 cd und höchstens 
625 cd betragen darf. Bei (x = 2,5°, ß = 1°), bezeichnet als „BR“, müssen 
mindestens 50 cd und höchstens 1750 cd vorhanden sein. 


3.2 DESIGN 


Die in dieser Arbeit verwendete CAD-Software und die Ray-Tracing- 
Software basieren auf „Non-Uniform Rational B-Spline“-Flächen 
(NURBS) [40]. Daher sollte das Ergebnis der Strukturierung wieder eine 
NURBS-Fläche sein. Eine kubische B-Spline Fläche (vgl. Kap.|2.2.2) ist 
ein Spezialfall von NURBS, bei dem nicht alle Freiheitsgrade verwendet 


werden. 


3.2.1 VERFAHREN ZUR STRUKTURIERUNG 


Das Ziel der Strukturierung ist eine lokale Verkippung der Flächennor- 
malen, um die Abbildungsqualität der Linse zu reduzieren. Dadurch 
wird eine scharfe HDG aufgeweicht. Die zu strukturierende Fläche wird 
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Seitenansicht Draufsicht 


Linse 


d Extrapolation 
Punktgitter 


Abbildung 3.2: Schematische Darstellung des ersten Schrittes der Strukturierung: Gene- 
rierung eines äquidistanten Punktgitters. 


um wenige Millimeter extrapoliert (typischerweise einen Millimeter), 
wie in Abb.\3.2/dargestellt. Die Extrapolation dient dazu, Randeffekte 
bei der späteren Flächenberechnung unwirksam zu machen (siehe un- 
ten). Es wird ein äquidistantes Punktgitter in der Ebene senkrecht zur 
optischen Achse generiert. Der typische Punktabstand im Gitter beträgt 
50 um bis 200 um. Das Punktgitter wird in die extrapolierte Fläche pro- 
jiziert, indem die Gitterpunkte als Startpunkte und die optische Achse 
als Richtungsvektor von Strahlen verwendet werden. Die Schnittpunkte 
zwischen Strahlen und Fläche werden, wie in Kap.|2.1.2|/beschrieben, be- 
rechnet. In Abb.\3.3\ist das projizierte Punktgitter schematisch dargestellt. 
Ausgehend von dem projizierten Punktgitter werden zwei Verfahren 


angewendet, um eine Strukturierung der Fläche zu erreichen. 


Bei dem ersten Verfahren werden die Punkte entlang der Normalen 
der Fläche unterschiedlich verschoben. Die Verschiebung ist in Abb. 3.4 
schematisch dargestellt. Der individuelle Versatz der Punkte wird in Ab- 
hängigkeit der HDG des Scheinwerfermoduls bestimmt. Die detaillierte 
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Seitenansicht Draufsicht 


Abbildung 3.3: Schematische Darstellung des zweiten Schrittes der Strukturierung: Pro- 
jektion des Punktgitters. 


Seitenansicht Draufsicht 


Abbildung 3.4: Schematische Darstellung des dritten Schrittes des ersten Verfahrens der 
Strukturierung: Verschiebung des Punktgitters entlang der Normalenvektoren. 
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Erklärung, wie der Versatz bestimmt wird, ist in Kap. gegeben. 
Aus dem Punktgitter wird eine Fläche bestimmt, indem eine interpolie- 
rende kubische B-Spline-Fläche berechnet wird. Die Berechnung ist in 
Kap.|3.2.2/erlautert. In Abb.|3.5/ist eine kubische Fläche aufgezeigt, die 
das Punktgitter aus Abb.|3.4linterpoliert. Die Parameter zur Erstellung 


der Struktur sind für Darstellungszwecke deutlich übertrieben gewählt. 
Wie in Abb.|3.6|dargestellt, werden die extrapolierten Bereiche durch 
zuschneiden entfernt. 


Seitenansicht Draufsicht 


Abbildung 3.5: Schematische Darstellung des vierten Schrittes des ersten Verfahrens der 
Strukturierung: Berechnung einer interpolierenden Fläche. 


Die zweite Variante zur Strukturierung verwendet das in die Fläche 
projizierte Punktgitter aus Abb. Anstatt die Punkte wie bei der 
ersten Variante zu verschieben, wird eine interpolierende Fläche von 
dem Punktgitter berechnet. Die berechnete Fläche ist näherungsweise 
identisch mit der Ausgangsfläche, wird aber mit einer größeren An- 
zahl an Kontrollpunkten (vgl. Gl. (2.11)) beschrieben. Die Abb.|3.7|zeigt 
eindimensional und schematisch die Fläche und ihre Kontrollpunkte. 


Anstatt die Punkte des Gitters vor der Flächenberechnung zu verschie- 
ben, werden direkt die Kontrollpunkte verschoben. Bei dem ersten Ver- 
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Seitenansicht Draufsicht 


ÆN 


Abbildung 3.6: Schematische Darstellung des fünften Schrittes der Strukturierung: Zu- 
schneiden der Fläche. 


fahren besteht ein direkter geometrischer Zusammenhang zwischen der 
Verschiebung der Gitterpunkte und der Fläche. Durch die Interpolati- 
on der verschobenen Punkte, wird die Fläche immer um den gleichen 
Betrag verändert wie die Punkte. Bei dem zweiten Verfahren ist diese di- 
rekte Abhängigkeit nicht gegeben. Durch die Verschiebung der Kontroll- 
punkte wird ebenfalls eine lokale Änderung verursacht, jedoch wird die 
Fläche nicht um den gleichen Betrag verschoben wie der Kontrollpunkt. 
Die zweite Methode hat den Vorteil der schnelleren Strukturierung 
der Fläche, wenn wiederholt die gleiche Fläche mit unterschiedlichen 
Strukturen versehen werden soll. Bei der ersten Methode muss für jede 
Strukturierung ein lineares Gleichungssystem (siehe Kap.|3.2.2) gelöst 


werden, bei der zweiten Methode nur einmalig. In Kap. Giwird die Ausle- 


gung der Strukturen durch Optimierungsalgorithmen erläutert. Bei dem 
Optimierungsprozess müssen fortlaufend neue Strukturen bestimmt 
werden. Die zweite Methode ist ungefähr zwanzig mal schneller bei der 
Strukturierung als die erste Methode. 
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=e Kontrollpolygon —— Strukturierte Fläche 
e Punktgitter der Ausgangsflache 


Abbildung 3.7: Schematische Darstellung des zweiten Verfahrens. Links: Interpolierende 
Fläche des Punktgitters. Rechts: Strukturierung der Fläche durch Verschieben der Kon- 
trollpunkte. 


3.2.2 BERECHNUNG EINER INTERPOLIERENDEN B-SPLINE-FLÄCHE 


Die zu interpolierenden Punkte eines Gitters der Größe (n+1)x (m+1) 
sind durch Ou gegeben. Um eine interpolierende Fläche S(u,v) zu 
erhalten, müssen die Knotenvektoren U und V sowie die Kontrollpunkte 
Pij bestimmt werden. Als Knoten werden die x- und y- Koordinaten der 
nicht verschobenen Gitterpunkte gewählt (vgl. Abb.|3.3). Dadurch wird 
die Fläche so parametrisiert, dass Parameterraum und Realraum ähnlich 


sind, also 


sel y (3.2) 
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Das hat Vorteile bei der Projektion von Punkten auf Flächen, wie in Kap. 
3.5|beschrieben. Zwischen der Anzahl an Knoten (r + 1) und der Anzahl 
an Kontrollpunkten (n + 1) gilt der Zusammenhang 


r+1=(n+1)+p+1 (3.3) 


Damit die äußeren Kontrollpunkte die Endpunkte der Fläche sind, müs- 
sen die äußeren Knoten (p + 1)-mal vorhanden sein. Dadurch resultiert, 
dass nicht jeder Gitterpunkt in dem Knotenvektor vorhanden sein kann. 
Die Koordinaten des zweiten und des vorletzten Punktes werden nicht 
in den Knotenvektor übernommen. Der Grad der Fläche wird in beide 
Richtungen auf Grad drei gesetzt (p = q = 3). Dadurch ergibt sich eine 
krümmungsstetige Fläche. Beispielsweise wird für eine Linse mit einem 
extrapolierten Durchmesser von 70 mm und einem Gitterabstand von 
0,1 mm der Knotenvektor 


U = {-35; —35; —35; —35; —34,8; —34,7; —34,6; —34,5; ..., 
34,5; 34,6; 34,7; 34,8; 35; 35; 35; 35} 


verwendet. Fiir eine Linse mit einem extrapolierten Durchmesser von 
80mm und einem Gitterabstand von 0,2 mm wird der Knotenvektor 


U = {—40; —40; —40; —40; —39,6; —39,4; —39,2; —39,0; ..., 
39,0; 39,2; 39,4; 39,6; 40; 40; 40; 40; } 


verwendet. Die äußeren Kontrollpunkte sind identisch mit den Gitter- 


punkten: 
Pro E Oxo Poy, = Orm fürk=0,...,n (3.4) 
Poy = Go, B, = On für l = 0, Br ltt (3.5) 


Zu bestimmen sind die inneren Kontrollpunkte für k = 1,...,n — 1 und 
l = 1,...,m — 1. Damit die Fläche die Punkte interpoliert, müssen die 
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Kontrollpunkte so bestimmt werden, dass die Gleichung 
Qk) = 5(2,,0)) = LE Nip (Hk) Nig OP; (3.6) 


erfüllt wird. Die Symbole a, und o geben die Parameter an, an denen 
die Interpolation stattfindet, bzw. an denen die Flache die Gitterpunkte 
bertihrt. Die Parameter werden identisch mit den Knoten u € U und 
v € V gewählt: 


üg = Uk+p-1 fürk=2,...,n—2 (3.7) 
Ül = Vl+g—1 für l = 2, seg = 2, (3.8) 


wodurch sich ein tridiagonales Gleichungssystem (siehe unten) aufstel- 
len lässt [40]. Die Verschiebung der Indizes um p — 1 ergibt sich durch 


die (p + 1)-fache Multiplizität der äußeren Knoten. Die Parameter für 
die Gitterpunkte mit k = 1,n — 1 und ! = 1,n — 1 werden durch den 
Mittelwert der benachbarten Parameter bestimmt: 


ak = “tai S "EL fürk = 1n-1 (3.9) 
5 = Ur su fir! =1,m—1 (3.10) 


Um die Kontrollpunkte zu bestimmen, sodass Gl. erfüllt wird, wird 
die Gleichung umgeschrieben zu [40]: 


n 


Ori A Nal (EN HIGH Bi) - Keng GAR: (3.11) 
i=0 


i=0 

mit Rii = -N Pj (3.12) 
Die Berechnung der interpolierenden Fläche wird durch mehrere 
Kurven-Interpolationen durchgeführt. Es werden m + 1 Kurven- 


Interpolationen durch die Punkte Qp,,...,Qn, (für | = 0,...,m) aus- 
geführt, um die R; zu erhalten. Weitere n + 1 Kurven-Interpolationen 
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durch R;o Pan 


von [40]): 


ak = Ny_13 (Ak) 


by = Nga (ñp) 


DESIGN 


+, Rim (für i = 0, . . ., n) bestimmen die P; j. Mit (modifiziert 


Ck = Ma (x) (3.13) 


wird für jede Kurven-Interpolation ein Gleichungssystem aufgestellt: 


mit 
a by 
0 a 
0 0 
0 0 
0 0 


Cy 


a3 


Qi 
Q2 
Q3 


Ro 
Rı 
R2 
R3 
=M. . (3.14) 
Rn-2 
Rn-1 
Ry 
M — 
0 0 0 0 
0 0 0 0 
0 0 0 0 
An—2 Dn-2 Cn-2 0 
N, 3,3 (iy 1) An—-1 Dunn Cn-1 


(3.15) 


Hierbei wird verwendet, dass an einem inneren Knoten nur drei Basis- 
Funktionen ungleich null sind. Da die Kontrollpunkte am Rand mit den 
Gitterpunkten übereinstimmen, wird Gl. (3.14) umgeformt zu (modifi- 


45 


PROZESSKETTE ZUR HERSTELLUNG MIKROSTRUKTURIERTER LINSEN 


ziert von [40]): 


Qi — a1Ro Rı 
Q2 R2 
Q3 R3 
=M | | (3.16) 
Qn-2 Rn-2 
Qn—1 — Cn-1Rn Rn-1 
mit 
by cy N33 (ii; ) 0 > 0 0 0 
ag bo C2 0 SSES 0 0 0 
0 a3 b3 C3 0 0 0 
M= . 
0 0 0 0 ++: An_2 ba—2 Cn-2 
0 0 0 O +++ Ny-33(@n-1) an1 bna 
(3.17) 


Mit zwei Umformungsschritten werden die Elemente N33 (#1) und 
N„-3,3 (@n—1) aus der ersten und letzten Zeile der Matrix entfernt. Da- 
durch ergibt sich ein tridiagonales Gleichungssystem, das mit dem 
Tridiagonalmatrix-Algorithmus gelöst werden kann [63]. 


3.2.3 BESTIMMUNG DER AMPLITUDEN 


Die Amplituden werden aus Normalverteilungen (Gauß-Kurven) ge- 
nerierf! Dazu wird das Punktgitter (vgl. Abb. 3.2) in Zellen eingeteilt. 
Die Einteilung ist für eine Zellengröße von 2x2 schematisch in Abb. 


3.8 dargestellt. Jedem Punkt aus der Zelle wird eine Normalverteilung 


1 Die Idee entstammt der alten Prozesskette (vgl. Kap. 1). 
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Abbildung 3.8: Einteilung der Gitterpunkte (vgl. Abb. 8.2) in Zellen der Größe 2x2 zur 
Zuordnung von Normalverteilungen. 


zugeordnet. Die Standardabweichungen der Normalverteilungen sind 
in jeder Zelle identisch. Der Versatz der Punkte (vgl. Kap. 3.2.1) wird 
aus den Normalverteilungen bestimmt. Die Parameter, über die das 


Streuverhalten der Struktur eingestellt wird, sind die Zellengröße, die 
Standardabweichungen o der Normalverteilungen und der Gitterab- 
stand. Um die Auswirkung der Parameter zu demonstrieren, wird im 
CAD eine planare Fläche der Größe 10 mm mal 10 mm strukturiert und 
in der Simulationsumgebung mit parallelem Licht bestrahlt. Die von den 
strukturierten Platten ausgehenden LSVen sind in den Abb.|3.9][3.11lund 
3.13|dargestellt. Die Abb.3.10/[3.12]und [3.14] zeigen ausschnittsweise die 
Strukturen. Die zum Erstellen der Strukturen verwendeten Parameter 


sind unter den einzelnen Abbildungen angegeben. Hierbei geben Ax 
und Ay die Auflösung der Punktgitter an. Die Verwendung von Zellen 
der Größe 1x1 ergeben Streuverhalten, die Normalverteilungen ähnlich 
sind (vgl. Abb. 3.9). Durch die Vergrößerung auf Zellen der Größe 2x2 
wird das Streuverhalten rechteckiger (vgl. Abb.|3.11). Um die Streuung 
in vertikaler oder horizontaler Richtung asymmetrischer zu gestalten, 


47 


PROZESSKETTE ZUR HERSTELLUNG MIKROSTRUKTURIERTER LINSEN 


Lichtstärke (wilk. Einh.) 


—2 0 2 a(°) -2 0 2a(°) 


()Ax=Ay=100pm co =0,5pm (b) Ax = Ay=100pm ez 0,7 um 


BU) 
2 
0 
-2 
-2 0 2a(°) 
(c) Ax = Ay=100pm o=1pm (d) Ax = Ay=100pm o=2pm 


Abbildung 3.9: Streuverhalten strukturierter Platten für Zellen der Größe 1x1. 
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Strukturhöhe (um) 


0.45 0.45 

= 0.15 => 0.15 

E E 

£ £ 

> -0.15 > -0.15 
-0.45 -0.45 


-0.45 -0.15 0.15 0.45 


-0.45 -0.15 0.15 0.45 
x (mm) 


x (mm) 
(a) Ax = Ay=100pm co =05pm 
0.45 


(b) Ax = Ay = 100pm o = 0,7 um 


-0.45 


-0.45 -0.15 0.15 0.45 
x (mm) 


-0.45 -0.15 0.15 0.45 
x (mm) 


(c) Ax = Ay = 100um o=1pm (d) Ax = Ay=100pm o=2pm 


Abbildung 3.10: Strukturhöhen der Platten für Zellen der Größe 1x1. Die grauen Linien 
kennzeichnen die Zellengrößen. 
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Lichtstärke (wilk. Einh.) 


-1 0 1 a(°) -1 0 Lol 


(a) Ax = Ay = 100 pm (b) Ax = Ay = 100 pm 


7 (95 oh D oh 
"Age 05)” pas 


-1 0 1 a(°) -1 0 Lol 


(c) Ax = Ay = 100um (d) Ax = Ay = 100 pm 


SI OV ft Na 
"Am atb "äi. oJ” 


Abbildung 3.11: Streuverhalten strukturierter Platten für Zellen der Größe 2x2. 
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Strukturhöhe (um) 


=. 
E 0.15 
> 
-0.35 
-0.45 -0.05 0.35 
x (mm) 
(b) Ax = Ay = 100 pm 
auf ZK sl ee 
"KR o5)” SEET Me 


-0.45 -0.05 0.35 
x (mm) 


-0.45 -0.05 0.35 
x (mm) 


(c) Ax = Ay = 100um (d) Ax = Ay = 100 pm 
En 10 A y= 1 0 m 
"Ae ay” Ab)" 
Abbildung 3.12: Strukturhöhen der Platten für Zellen der Größe 2x2. Die grauen Linien 


kennzeichnen die Zellengrößen. Die Farbskala ist unterhalb von —1,5 um und oberhalb 
von 1,5 um logarithmisch. 
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Lichtstärke (wilk. Einh.) 


e H 
(a) Ax = 100 um, Ay = 500 pm (b) Ax = 100 um, Ay = 500 pm 
o=1lpm 10 
T= pm 
0 1 
pe) p (°) 


1.5 
1.0 
0.5 
0.0 
—0.5 
—1.0 


1.5 
1.0 
0.5 
0.0 
—0.5 
—1.0 


-1.5 =1:5 
"or 
-0.5 0 0.5 -0.5 0 0.5 
(c) Ax = 500 um, Ay = (d) Ax = 500 um, Ay = 
100 pm 100 pm 
go=1pm Il 0 
Ge 


Abbildung 3.13: Streuverhalten strukturierter Platten für Punktgitter mit unterschiedli- 
chen Auflösungen in x- und y-Richtung. 
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Strukturhöhe (um) 


y (mm) 
y (mm) 


(a) Ax = 100 um, Ay = 500 pm (b) Ax = 100 um, Ay = 500 pm 


= SF 


y (mm) 
y (mm) 


2 
-1.75 -0.75 0.25 1.25 


-1.75 -0.75 0.25 1.25 
x (mm) 


x (mm) 


(c) Ax = 500 pm, Ay = 100 pm (d) Ax = 500 pm, Ay = 100 pm 


go=1pm 1 0 
T= pm. 
0 1 


Abbildung 3.14: Strukturhöhen der Platten für Punktgitter mit unterschiedlichen Auf- 
lösungen in x- und y-Richtung. Die grauen Linien kennzeichnen die Zellengrößen. Zur 


besseren Sichtbarkeit sind nur die Zellengrenzen in der Richtung der größeren Ausdeh- 
nung eingezeichnet. 
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können die Gitterabstände in x- und y-Richtung unterschiedlich gewählt 
werden (vgl. Abb. B-13). Lichtstärkeverteilungen von Scheinwerfern mit 
Matrix-LED-Systemen haben neben horizontalen auch vertikale HDGs 
[64-66]. Um je nach Anforderungen die HDGs unterschiedlich aufzu- 
weichen, können die Strukturen mit asymmetrischem Streuverhalten 


verwendet werden. 


3.3 FERTIGUNG 


Die Linsen werden im Spritzgussverfahren hergestellt [15]. Der thermo- 


plastische Kunststoff PMMA wird erhitzt und in den Hohlraum eines 
Stahlwerkzeuges gespritzt [67]. Nach dem Abkühlen nimmt das Materi- 
al einen festen Zustand an. Durch die Form des Werkzeuges bzw. des 


Werkzeugeinsatzes wird die Linse und ihre Mikrostruktur vorgegeben. 
In Abb. sind die in dieser Arbeit verwendeten Werkzeughälften 
dargestellt. Die Werkzeugeinsätze mit der strukturierten Lichtaustritts- 


Abbildung 3.15: Links: Werkzeughälfte mit Einsatz für die Lichtaustrittsseite der Linse 
sowie spritzgegossene Linsen. Rechts: Werkzeughälfte mit Lichteintrittseite der Linse. 


seite können gewechselt werden. Ein Werkzeugeinsatz ist in Abb.|3.16 


dargestellt. Die Einsätze werden durch konventionelle Verfahren ohne 
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Abbildung 3.16: Werkzeugeinsatz mit Mikrostruktur. 


Struktur und mit 50 um Aufmaß hergestellt. Die Struktur wird durch 
eine Drehmaschine mit einem Diamantschneidewerkzeug in den Einsatz 
eingebracht [68]. Um die Bearbeitung eines Stahleinsatzes durch ein Dia- 
mantschneidewerkzeug zu ermöglichen, muss der Diamant durch Ultra- 
schall in Schwingung versetzt werden [69]. Die Schwingung unterbricht 
mit einer Frequenz von der Größenordnung im 100 Kilohertz-Bereich 
den Schneidvorgang, wodurch der Verschleiß des Schneidwerkzeugs 


soweit reduziert wird, dass eine Bearbeitung eines Stahlbauteils erst er- 
möglicht wird [70-72]. In Abb.|3.17|ist die im Förderprojekt verwendete 
Drehmaschine und ein Diamantschneidwerkzeug dargestellt. 


Im Rahmen des Förderprojektes wurde von der Firma Innolite ein neu- 
artiger Steuerungsansatz für Drehmaschinen entwickelt, um eine höhere 
Präzision bei der Fertigung zu erreichen [74, 75]. Die Drehmaschine 


fertigt dabei auf Basis eines polaren Punktgitters, wie für eine kleine 
Punktanzahl schematisch in Abb. |3.18| dargestellt. Eine typische Auf- 
lösung des Punktgitters für eine mikrostrukturierte Linse mit einem 
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Abbildung 3.17: Links: Im Projekt verwendete Drehmaschine „Innolite IL 300” Kei Rechts: 
Bauteil und Diamantwerkzeug auf der Drehmaschine [74]. 


Abbildung 3.18: Schematische Darstellung eines polaren Punktgitters. 
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Mahr 


Aspheric 3D 


MarForm MFU 200 


Abbildung 3.19: Links: In dem Projekt verwendete Messmaschine Mahr MFU200 eil 
Rechts: Werkzeugeinsatz auf Messmaschine [61]. 


Durchmesser von 70 mm beträgt 0,05° und 5 um. Dadurch ergeben sich 
ca. 50 Millionen Punkte. Diese Punkte werden aus der strukturierten Lin- 
senfläche generiert und in dem Dateiformat „Hierarchical Data Format 
5“ (HDF5) an die Fertigung übergeben [76]. 


3.4 VERMESSUNG 


Um zu überprüfen, ob die Fertigung der mikrostrukturierten Werkzeu- 
geinsätze erfolgreich ist, werden die Einsätze optisch vermessen. Für 
die Messung der Einsätze wird das System MFU200 der Firma Mahr 
verwendet, wie in Abb. 3.19]dargestellt [77]. Das Messsystem ist eine 
zylindrische Koordinaten-Messmaschine, bestehend aus zwei linearen 
(X, Z) und einer drehbaren Messachse (C). Die Klemmvorrichtung zur 


57 


PROZESSKETTE ZUR HERSTELLUNG MIKROSTRUKTURIERTER LINSEN 


Aufnahme des Werkzeugeinsatzes kann durch vier Achsen verschoben 
und gekippt werden, um den Einsatz zentral und orthogonal zum Mess- 
sensor zu platzieren. Zur Bestimmung des Zenits kann eine zusätzliche 
lineare Achse (Y) verwendet werden. Das Messverfahren basiert auf 
Weißlichtinterferometrie. Der Messsensor hat einen Arbeitsbereich von 


+20° zur Flächennormalen bei einem Abstand von (300 +5) pm. Um 
die äußeren Bereiche der Einsätze, mit einer Steilheit größer 20° zu ver- 
messen, kann der Messsensor gekippt werden (B-Achse). Im Fokus des 
optischen Systems beträgt der Strahldurchmesser 2 um. 


Zur Zentrierung des Einsatzes wird entlang zweier Linien gemessen, wie 
in Abb.|3.20 dargestellt. Der Scheitelpunkt wird über Ausgleichspolyno- 


Abbildung 3.20: Bewegung des Messkopfes zur Justierung (A) und Messung (B) [e1]. 


me bestimmt. Die Vermessung des Einsatzes erfolgt durch die Messung 
von konzentrischen Kreisen. Die X- und Z-Achsen fahren nach Soll- 
Daten eine Position an, die zu einem vorgegebenen Radius gehört. An 
der Stelle wird um die C-Achse mit konstanter Geschwindigkeit gedreht 
und währenddessen mit konstantem Winkelabstand Messwerte aufge- 
nommen. Die in dieser Arbeit verwendete Schrittweite beträgt 0,1° und 
25 um. Bei den verwendeten Linsen mit 70mm Durchmesser ergeben 
sich Punktgitter mit ca. 5 Millionen Punkten. Diese Messauflösung stellt 
einen Kompromiss zwischen Detailgrad und Messzeit, die hierbei bei 
mehreren Stunden lag. 
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3.5 FLÄCHENRÜCKFÜHRUNG 


Um die Verwendbarkeit der gefertigten Werkzeugeinsätze zu überprü- 
fen, ist ein Soll-Ist-Vergleich zwischen Messdaten und Fertigungsdaten 
nicht sinnvoll [78} 79]. Bei dem Vergleich können zwar grobe Unterschie- 
de ausgemacht und daraus eventuelle Probleme bei der Fertigung abge- 
leitet werden, jedoch ist das photometrische Ergebnis der abgeformten 
Linse entscheidend. Daher muss auf Basis des gemessenen Punktgitters 
eine lichttechnische Simulation durchgeführt werden. Damit kann die 
LSV der abgeformten Linse bewertet werden. Zur Bestimmung der LSV 
durch eine Ray-Iracing-Simulation wird Information über die Fläche 
zwischen den Punkten benötigt. Daher muss eine Flächendarstellung 
des Punktgitters bestimmt werden. 

Eine Interpolation der Punkte, wie in Kap. 3.2]beschrieben, ist unge- 
eignet. Die Interpolation würde das Messrauschen mit in die Fläche 
übernehmen und benötigt für jeden Messpunkt einen Kontrollpunkt 
(vgl. Kap. 3.2.2). Durch die Ubernahme des Messrauschens wiirde die 
Abbildungsqualität der Linse in der Simulation schlechter sein als sie 


tatsächlich ist. Eine Fläche mit ca. 5 Millionen Kontrollpunkten ist im 
Ray-Tracing nicht gut handhabbar. Zum Beispiel würde die Berechnung 
einer Hierarchie von Begrenzungsvolumen (vgl. Kap. 2.1.2) und das an- 
schließende Ray-Iracing vergleichsweise lange dauern. Daher muss eine 


bessere Darstellung der Fläche bestimmt werden, die das Punktgitter 
bestmöglich approximiert. 


3.5.1 ÜBERSICHT VERSCHIEDENER VERFAHREN 
ZUR RÜCKFÜHRUNG 


Wie bei der Interpolation müssen Kontrollpunkte, Knoten und eventuel- 
le Parameter für die Messpunkte bestimmt werden. Für das Anpassen 
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einer Fläche an ein Punktgitter existieren verschiedene Ansätze. PIEGEL 
und TILLER berechnen eine Fläche, indem entlang zweier Richtungen 
mehrere Kurvenanpassungen durchgeführt werden. Das entspricht dem 
Vorgehen aus Kapitel[3.2.2} jedoch mit Methode der kleinsten Quadra- 
te HO, 80]. Die Knoten und Parameter werden vorab berechnet. Für 


den Algorithmus müssen die Messdaten gitterbasiert vorliegen. Diese 
Methode wurde durch HEMMERICH im Bezug auf mikrostrukturierte 
Flächen hin untersucht [81]. In der Arbeit von Hemmerich wurde die Flä- 
chenrückführung erfolgreich für kartesische Punktgitter durchgeführt. 


Eine Anwendung auf polare Gitter ist möglich, indem die Kurvenanpas- 
sungen in radialer und azimuthaler Richtung verlaufen. Entsprechend 
verlaufen die u- und v- Richtungen der resultierenden Fläche bzw. des 


Kontrollpolygons radial und azimuthal, wie in Abb. 3.21! dargestellt. 


Bei dem Importieren trennt CATIA eine so parametrisierte Fläche in 
zwei Flächen auf. Das ist zwar kein direktes Problem, jedoch zeigt sich, 


SLYN \ 
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Abbildung 3.21: Kontrollpunkte (rot) und Kontrollpolygon (blau) einer B-Spline Fläche. 
Die Lücke zwischen uo und u10 dient der Veranschaulichung. 


dass CATIA polar-parametrisierte Flächen beim Einladen umrechnet. 
Beim Umrechnen besteht die Möglichkeit der Datenverluste. Ein wei- 
terer Nachteil ist die ungleiche Verteilung der Kontrollpunkte. Durch 
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das polare Gitter sind im Inneren der Fläche die Abstände der Kon- 
trollpunkte geringer als im Randbereich. Eine ausreichende Anzahl an 
Kontrollpunkten im Randbereich erfordert eine unnötig hohe Anzahl an 
Kontrollpunkten im Inneren. 


BECKER ET AL. führen eine zweistufige Flächenanpassung durch, die 


für nicht-strukturierte Gitter verwendet werden kann [82]. In der ersten 
Stufe wird im Gegensatz zu der Methode von PIEGEL und TILLER ein 
komplexerer Algorithmus zur Bestimmung der Parameter verwendet. 
Durch Aufstellen und Lösen eines Gleichungssystems werden die Po- 
sitionen der Kontrollpunkte bestimmt. In dem zweiten Schritt wird als 
Flächendarstellung ,, NURBS” (Non-Uniform Rational B-Spline) verwen- 


det [40]. Durch NURBS ergeben sich weitere Freiheitsgrade, sogenannte 


Gewichte. Mit einem Gradientenverfahren werden die Gewichte und die 
Positionen der Kontrollpunkte optimiert. Die Methode ist für nominelle 
Daten ausgelegt, die nicht durch ein Messrauschen verunreinigt sind. 


ZHANG ET AL. wählen aus den Messdaten einige Punkte aus [83]. Die 
Punkte dienen als Stützstellen zur Bestimmung einer interpolierenden 
B-Spline-Fläche. Um die Fläche an die Messdaten anzupassen, wird das 
Newton-Raphson-Verfahren angewendet. Bei dem Verfahren werden 
die Positionen der Stützstellen verschoben, bis die Differenz zwischen 
Fläche und Messdaten minimal ist. Es werden keine Informationen 
gegeben, wie die Stützstellen bzw. ihre Anzahl gewählt werden sollen. 


In der Arbeit von CARLSSON wird eine NURBS-Flächenrückführung 
durch die Gauss-Newton-Methode untersucht [84]. Variiert werden Pa- 
rameter, Kontrollpunkte und Gewichte. Die Anzahl der Messpunkte 
beträgt ca. 500 und die Größe der Kontrollpunktmatrix der Fläche 7x7. 
Bei dem Verfahren wird die Jacobi-Matrix aufgestellt. Die Größe der 
Matrix wächst in beiden Dimensionen linear mit der Anzahl an Mess- 


punkten. Da bei der Rückführung der vermessenen Werkzeugeinsätze 
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die Anzahl an Messpunkten im Bereich mehrerer Millionen liegt, ist die 
Methode nicht anwendbar. 


RANDRIANARIVONY und BRUNNEN nutzen beschränkte nichtlineare 
Optimierung, um neben den Positionen der Kontrollpunkte und der 
Gewichte auch die Knoten einer NURBS-Fläche zu optimieren [8586]. 
Es werden keine Informationen gegeben, wie die Parameter der Mess- 


punkte bestimmt oder die Anzahl der Kontrollpunkte gewählt werden 
soll. 


BRUJIC ET AL. stellen eine Methode vor, um eine bereits angepasste 
Fläche noch besser anzupassen [87]. Zur Berechnung der neuen Kon- 
trollpunkte wird die Funktion 


MI, | 2 nom, | 2 
f= X (Sa) rot YO (By - 3,0) 
k=0 i=0 j=0 
(3.18) 
n m = 0 
We 
i=0 j= 


bezüglich den neuen Kontrollpunkten Pj minimiert. Der erste Term 
bestraft den Abstand zwischen der neuen Fläche S, beschrieben durch 
die Pj, und den Messpunkten Ok. Die Parameter tix, Öp ergeben sich 
durch Projektion der Messpunkte auf die Startflache, beschrieben durch 
die Kontrollpunkte po. Der zweite Term bestraft den Abstand zwischen 
Flache und den neuen Kontrollpunkten, wodurch die neue Flache glatter 
verläuft. Die Stärke der Glättung kann durch den Parameter y eingestellt 
werden. In Situationen, wo die Form der Startflache in Bereichen einer 
diinnen Messpunktdichte erhalten bleiben soll, kann der dritte Term 
verwendet werden. Der dritte Term bestraft den Abstand zwischen den 
Kontrollpunkten der neuen Fläche zu denen der Startfläche. Wie stark 
sich Start und Endfläche ähnlich seien sollen, kann durch den Parameter 
ö eingestellt werden. 
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3.5.2 ITERATIVE GEOMETRISCHE APPROXIMATION 


Zur Rückführung wird der „Iterative Geometrische Approximations”- 
Algorithmus verwendet [88]. Der Algorithmus hat die Vorteile, dass er 
ohne Aufstellen bzw. Lösen eines Gleichungssystems auskommt, keinen 


zusätzlichen Term zum Glätten (z. B. Biegeenergie) benötigt und ermög- 
licht, die Approximations-Genauigkeit lokal beeinflussen zu können. 
Zur Flächenrückführung wird eine Ausgangsfläche, beschrieben durch 
B-Splines, benötigt. Der Algorithmus berechnet den Abstand zwischen 
Messdaten und Fläche und verschiebt die Kontrollpunkte entsprechend 
des Differenzvektors. 


Die Ausgangsfläche ist durch $ (u,v) und die Fläche nach jedem Ite- 
rationsschritt durch 5% (u,v) beschrieben. Der Index « gibt die Anzahl 
der Iterationsschritte an. Die Messdaten werden durch On, k=0,...,N 
beschrieben. Fiir jeden Messpunkt wird die orthogonale Projektion auf 
die Fläche Silo) (u,v) bestimmt, bzw. die Parameter ü E al”, die den pro- 
jizierten Punkt beschreiben. Wie die orthogonale Projektion bestimmt 


wird, ist in Kap. 3.5.3|gegeben. Wie in Kap.|2.2.2}geben us und o: die Ele- 
(a) (a) 
k 1% 


mente der Knotenvektoren der Fläche an. Zu den Parametern ü 
werden die Indizes s und t so bestimmt, dass u, < a”) < Us+ı und v; < 
al”) < v;41 gilt. Der projizierte Punkt Silo! (a ; a) berechnet sich aus 
den nicht verschwindenden Basis-Funktionen Nip (a), i =S—p,...,8 
und Nj (0 ),j = t-q,...,t und damit aus den Kontrollpunkten 
e i= s—p,...,8,] = t—q,...,t. Die Abhängigkeit ist in Abb. 
3.22|veranschaulicht. Zur Anpassung der Fläche an die Punkte werden 


die Differenzvektoren Q; — 5% Ciel a) bei jedem Iterationsschritt 
berechnet. Die Differenzvektoren werden mit den Basis-Funktionen ge- 
wichtet und zu den Kontrollpunkten hinzu addiert. Dadurch wird die 
Fläche iterativ näher an die Punkte herangezogen. Ein Iterationsschritt 


ist in der Auflistung)3.1|beschrieben. 


63 


PROZESSKETTE ZUR HERSTELLUNG MIKROSTRUKTURIERTER LINSEN 


Programmauflistung 3.1: Iterative Geometrische Approximation [ss] 


ı for k+-0,N 
ig, ö + getProjectedPoint (QIk]) 


2 

8 s, t4 getSpan (ak, Ok) 

ı forie-s-p,5 

5 fOr j= t= D 

6 Wans) = Wnum [illj] + Ni, pa) Nig (Ox) - 
(SI) 5) (ii, ö%)) 

7 Waenlilli] — Waenlilli] + Maul Neal 

8 end for 


9 end for 
10 end for 
1 for i¢O,n 


2 for j+0,m 
3(4+1) pa) ` Well 
Du Fig "Wl 


14 end for 
15 end for 


Die Iteration wird so lange durchgeführt, bis die Anderung der mitt- 
leren Abweichung zwischen Fläche und Punktgitter kleiner als 1 nm 
beträgt. Es wird angenommen, dass ab einem Unterschied von Inm 
die Abweichungen keine lichttechnische Auswirkung mehr haben. Wie 
beschrieben, benötigt der Algorithmus eine Ausgangsfläche. In dieser 
Arbeit wird dafür die hoch-aufgelöste unstrukturierte Fläche (vgl. Kap. 
3.2.1| „zweite Variante”) verwendet. 
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Abbildung 3.22: Ein bikubisches Kontrollpolygon [ss]. Die schwarzen unterbrochenen 
Linien geben die Positionen der inneren Knoten an. Das schwarze Kreuz markiert den 
projizierten Punkt. Die blau markierten Kontrollpunkte beeinflussen die Fläche an der 
Stelle mit dem schwarzen Kreuz. 


3.5.3 ORTHOGONALE PROJEKTION AUF FLÄCHEN 


Zur Projektion der Punkte auf die Fläche wird ein Algorithmus verwen- 
det, der die Krümmung der Fläche nutzt, um die Projektionspunkte zu 
bestimmen [89, 90]. Der Punkt, dessen Projektion auf die Fläche g (u,v) 
gesucht wird, ist durch ¥ gegeben. Der Algorithmus benötigt einen 
Startpunkt fo. Der Startpunkt po stellt eine erste Vermutung des Projekti- 
onspunktes dar und ist im Parameterraum durch die Koordinaten (ü,ö) 
gegeben. Die Tangentialvektoren von der Fläche an der Stelle jig sind 5} 
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Abbildung 3.23: Schematische Darstellung der orthogonalen Projektion (modifiziert von 


ei 


und 52. In Abb. 3.23list der Aufbau veranschaulicht. Der Normalenvek- 
tor ist H = S, x Sp. Um den geratenen Projektionspunkt fu genauer zu 


bestimmen, wird der Differenzvektor X — po als Linearkombination aus 
Tangentialvektoren und Normalenvektor dargestellt [89]: 


X — Po = M51 + A252 +vn (3.19) 


Durch das Lösen des Gleichungssystems (3.19) ergeben sich die Parame- 
ter An, A2 und v. Die zweifachen Ableitungen der Fläche sind gegeben 
durch 


9? 2 - ` 2 
Sap le?) So = 502 (M2) 


(3.20) 


Mit den Koeffizienten der ersten und zweiten Fundamentalform 


8; = (8,5;) und hy = (Sit) (3.21) 
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berechnet sich die Krümmung «x in Richtung des Tangentialvektors = 
A151 St A252 aus [89] 


X2. (hyA;A; 
x= H ch d (3.22) 
Bye (gijAiAj) 
In der Ebene, die die Punkte fo, X und den Vektor ñ enthält, befindet 
sich der Krümmungskreis mit Krümmung x und dem Mittelpunkt 


> 


Ze = Po + =. (3.23) 


Die orthogonale Projektion des Punktes ¥ auf den Kriimmungskreis 


ergibt den Punkt 


> > 


X-% 


d Cet nn, 
HT RE] 


(3.24) 


Mit dem Projektionspunkt 7 wird die Änderung Ah der Parameter be- 
rechnet aus [89]: 


ET 
Ah = sgn ((&%,q — Bal, el (3.25) 
K 


Hierbei gibt sgn() die Vorzeichenfunktion an. Die Parameter werden 
aktualisiert durch [89]: 


ū + ū + Ah (3.26) 
oc DL Ah (3.27) 
po + 5(n,ö) (3.28) 


Die Schritte bis werden solange wiederholt, bis die ge- 
wünschte Genauigkeit erreicht ist. Als Konvergenzkriterium wird ver- 
wendet, dass der Kosinus des Winkels zwischen ñ und X — po kleiner 
als 1078 oder der Abstand zwischen X und jo kleiner als 1 x 1071? mm 
ist. Die Werte wurden empirisch gefunden. Bei größeren Werten endete 
die Konvergenz des Algorithmus vermutlich auf Grund der endlichen 
Rechengenauigkeit. 
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3.5.4 UNTERSUCHUNG DER KONVERGENZ 


Zur Überprüfung der Konvergenz wird eine Linse mit einem Radius 
von 35mm mit einer Mikrostruktur versehen. Für die Erstellung der 
Mikrostruktur wird ein quadratisches Gitter mit einem Punktabstand 
von 100 um verwendet. Die Amplituden werden aus einer Normalver- 
teilung mit einer Standardabweichung von 0,5 pm bestimmt. Von der 
strukturierten Linsenfläche wird ein polares Punktgitter mit einer Auf- 
lösung von 0,1° und 25 pm abgeleitet. Diese Auflösung entspricht der 
Messauflösung (vgl. Kap. 3.4). Zum Starten der Flächenrückführung 


wird eine unstrukturierte Linsenfläche verwendet, deren Kontrollpunk- 
te einen lateralen Abstand von AP,,, = 100 um haben und damit der 
Auflösung der strukturierten Linse entsprechen. In Abb. (blaue 
Kurve) ist die mittlere orthogonale Abweichung zwischen Punktgitter 
und Fläche in Abhängigkeit der Iterationsschritte dargestellt. Der Algo- 


ORauschen = 100 nm 
ORauschen = 200 nm 


ORauschen = 0 nm 
—— rauschen = 100 nm, AP,,, = 254m 


10° 


Mittlere absolute 
Abweichung (um) 


T T T T T 
0 200 400 600 800 1000 
Iterationsschritte 


Abbildung 3.24: Konvergenzverhalten des Algorithmus Groe eine mikrostrukturierte 
Linse. 
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rithmus erreicht nach 1100 Iterationsschritten eine mittlere Abweichung 
von 0,5nm, wobei die Konvergenz noch nicht erreicht wird. Theoretisch 
ist eine Konvergenz bis zur Rechnergenauigkeit möglich, da es sich hier 
um nominelle Daten handelt. Für weitere Tests wird das Punktgitter 
mit einem Rauschen versehen, indem jeder Punkt entlang der optischen 
Achse der Linse verschoben wird. Die Verschiebung wird dabei aus 
einer Normalverteilung generiert. Die Standardabweichung der Vertei- 
lung beträgt bei einem Versuch 100 nm und bei einem anderen Versuch 


200 nm. Das Konvergenzverhalten bei der Rückführung ist in Abb. 3.24 


dargestellt. Der Algorithmus konvergierte gegen eine mittlere ortho- 
gonale Abweichung von 65nm bzw. 130 nm. Wird die Auflösung der 
Ausgangsfläche auf einen Abstand von 25 um erhöht, so reduzierte sich 
die mittlere orthogonale Abweichung für das Punktgitter mit 100 nm 
Rauschen auf 22nm (rote Kurve in Abb. 3.24). Das bedeutet, dass ein 
Teil des Rauschens durch den Algorithmus in die Fläche übernommen 


wurde. Durch den geringeren Abstand der Kontrollpunkte von 25 um 
anstatt 100 um sind sechzehn Mal so viele Freiheitsgrade vorhanden. 
Um die Auswirkungen von einer Uberanpassung zu demonstrieren, 
wird die nominelle strukturierte und alle zurückgeführten Linsenflä- 
chen in einer lichttechnischen Simulation für ein Scheinwerfermodul 
(vgl. Kap. 4) verwendet. Die simulierten Gradienten-Kurven (vgl. Kap. 
8.1) sind in Abb. |3.25/dargestellt. Die Gradienten-Kurve der nominel- 
len Linsenfläche stimmt im Rahmen des Simulationsrauschens mit den 


Gradienten-Kurven der Rückführungen überein, die auf Basis der Flä- 
che mit einer Auflösung von 100 um durchgeführt wird. Die Simulation 
mit der 25 um-Ausgangsfläche zeigt auf Grund der Überanpassung 
einen weicheren Gradienten. Daher ist es wichtig, die Auflösung der 
Ausgangsfläche so gering wie möglich zu halten, um kein potentielles 
Messrauschen mit in die Fläche zu übernehmen. 


Bei den Rückführungen von realen Messdaten wird die Auflösung der 


nominellen Fläche verwendet. Wie in Kap.|4|gezeigt wird, ist die Auf- 
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lösung im Rahmen dieser Arbeit immer ausreichend. Sollte bei einer 
Rückführung die Auflösung nicht ausreichend sein, so kann durch Hin- 
zufügen von weiteren Spalten oder Zeilen (bzw. hinzufügen weiterer 
Knoten) in der Kontrollpunktmatrix die Anzahl an Freiheitsgrade der 
Fläche und damit die Genauigkeit der Rückführung erhöht werden. 


— Nominell 
Ban —— Rauschen = 0 nm 
e 0.24 —— ORauschen = 100 nm 
E —— ORauschen = 200 nm 
0.14 Orauschen = 100 nm, 
— APy,y=25um 
0.0 “a T T T T 
-1.00 -0.75 -0.50 -0.25 0.00 


BO) 


Abbildung 3.25: Simulierte Gradienten-Kurven in Abhängigkeit eines künstlichen Mess- 
rauschens und der Auflösung der Ausgangsfläche. 
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ANWENDUNG DER PROZESSKETTE 


In Abb. |4.1ist eine Übersicht über die neue Prozesskette dargestellt. An- 


hand von zwei verschieden Mikrostrukturen wird die Funktion der neu- 


en Prozesskette demonstriert. Eine Mikrostruktur wurde für den ECE-, 
die zweite für den FMVSS-Bereich ausgelegt. Üblicherweise sind die 
HDGs von Scheinwerfersysten für den FMVSS-Bereich deutlich weicher. 
Die Mikrostrukturen wurden manuell, d h. in einem „trial-and-error”- 
Verfahren, von einem Optikingeniuer ausgelegt. Die Funktionweise der 
automatischen Korrekturschleife (vgl. Abb. 4.1) wird in Kap.\5lerläutert 
und demonstriert. Der Inhalt dieses Kapitels stammt zu Teilen aus der 
Veröffentlichung [61]. Es werden die zwei Mikrostrukturen dargestellt, 


zu denen jeweils zwei Werkzeugeinsätze gefertigt wurden. Die Ergebnis- 
se der Flächenrückführungen und die darauf basierenden Simulationen 
werden präsentiert. 


4.1 VERMESSUNG VON LICHTSTÄRKEVERTEILUNGEN 


Zur Vermessung von Lichtstärkeverteilungen existieren mehrere Mög- 
lichkeiten [92]. Eine vergleichsweise schnelle Variante ist die indirekte 
Messung durch eine Bildverarbeitungs-Anlage (vgl. Abb. |4.2). Der 
Scheinwerfer wird in einem abgedunkelten Raum in 10 m Entfernung 


vor einer Messwand platziert. Das Scheinwerferlicht verursacht eine 
Leuchtdichte auf der diffus reflektierenden Messwand. Mit einer Leucht- 
dichtemesskamera wird ein Leuchtdichtebild aufgenommen. Über die 
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Abbildung 4.1: Neue Prozesskette zur Herstellung mikrostrukturierter Linsen. 


Leuchtdichte kann auf die LSV des Scheinwerfers zurückgerechnet 


werden. 


4.2 FUR ECE-GEBIETE 


Zur Uberpriifung der neuen, in Kap. B|vorgestellten Prozesskette wird 
das in Abb. |4.3/dargestellte Scheinwerfermodul BiLED verwendet. Das 
Modul besteht aus insgesamt fiinf LED-Lichtquellen, die das Abblend- 
licht erzeugen. Der aus drei Kammern bestehende Reflektor reflektiert 
das Licht tiber eine Blende in Richtung der Linse. Die Linse bildet die 
Blendenebene auf die Strafse ab und erzeugt hierdurch die HDG. Durch 
Vermessungen des BiLed-Moduls mit einer nicht-strukturierten Linse 
auf einer BV-Anlage, zeigte sich, dass oberhalb der HDG Streulicht von 
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Scheinwerfer 


Abbildung 4.2: Schematische Darstellung einer Bildverarbeitungs-Anlage. 


ca. 75 cd auftritt. Daher werden zu allen Simulationsergebnissen 75 cd 
hinzugefügt. Die von dem Modul ohne strukturierte Linse simulierte 
LSV sowie die Gradienten-Kurve sind in Abb.|4.4 dargestellt. 


Die ausgelegte Mikrostruktur ist in Abb. 14.5| dargestellt. Wie in Kap. 
3.2.1]beschrieben, basiert die Struktur auf einem Punktgitter, wobei hier 
eine Auflösung von 100 um verwendet wird. Das BiLed-Modul erzeugt 
durch zusätzliche LEDs ein Fernlicht. Die Fernlicht-LEDs sind in Abb. 
4.3|nicht dargestellt. Das Fernlicht wird wie das Abblendlicht durch die 
Linse gebrochen. Wie in Abb.\4.5lerkennbar, hat die Struktur in einem 


horizontalen Band eine größere Amplitude. Die unterschiedliche Aus- 
prägung der Struktur ist nötig, um die HDG ausreichend aufzuweichen, 
ohne das Fernlicht des Systems zu stark zu beeinflussen. Die simulierte 
LSV und die Gradienten-Kurve sind in Abb. 14.6a| dargestellt. Wie in 
Kap. |3lbeschrieben, wird ein polares Punktgitter mit ca. 50,4 Millionen 


Punkten exportiert und zum Herstellen von zwei Werkzeugeinsätzen 
durch eine Diamantdrehmaschine verwendet. Die Werkzeugeinsätze 
werden vermessen, wodurch sich zwei polare Punktgitter mit je ca. 5 
Millionen Punkten ergeben. Für jedes gemessene Punktgitter wird eine 
Flächenrückführung (vgl. Kap. B.5.2) durchgeführt. Die Abweichungen 
der rückgeführten Fläche zu den Messdaten ist in Abb. 4.7 dargestellt. 
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Reflektor LED-Lichtquellen Blende 


Abbildung 4.3: Aufbau des BiLED-Moduls, das zum Testen der neuen Prozesskette ver- 
wendet wird. 
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Abbildung 4.4: Simulation der LSV und der Gradienten-Kurve des BiLED-Moduls 
mit nicht-strukturierter Linse. Die gestrichelte Linie in der LSV gibt die Position der 
Gradienten-Kurve an. 


Die mittlere Abweichung zwischen Fläche und Punktgitter beträgt 32 nm 
für Werkzeugeinsatz 1 und 47 nm für Werkzeugeinsatz 2. Wie in Abb. 


4.7|zu erkennen ist, sind in den Messdaten von Werkzeugeinsatz 2 kreis- 


förmige Abweichungen mit einer Größe von ca. 0,1 um enthalten. Diese 
können entweder durch einen Fertigungs- oder Messfehler entstehen. 
Die Farbskala von Abb. |4.7|hat einen maximalen Wert von 10 pm, da 
vereinzelt Ausreißer in der Größenordnung von Mikrometern in den 


Messdaten vorhanden sind. Die Ausreißer können durch Verunreinigun- 
gen (z. B. Staub) des Werkzeugeinsatzes beim Messvorgang verursacht 
worden sein. 


Die simulierte LSV der rückgeführten Flächen sowie die dazugehöri- 
gen Gradienten-Kurven sind in Abb. /4.6b| und Abb. 14.6c! dargestellt. 
Bei beiden Werkzeugeinsätzen ist der Gradient von der rückgeführten 


Fläche um 0,06 größer als nominell ausgelegt. Aus beiden Werkzeuge- 
insätzen wurden Linsen im Spritzguss-Verfahren abgeformt, in einem 
Scheinwerfermodul verbaut und die erzeugte LSV photometrisch auf 
einer Bildverarbeitungsanlage vermessen. Die Gradienten-Kurven der 
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x (mm) 


54372 -1 05 -0.235 0 025 05 1 2 345 
z (um) 


Abbildung 4.5: Nominelle Mikrostruktur zur Aufweichung der HDG für den ECE-Bereich. 
Die Farbskala ist unter —0,5 um und über 0,5 ym logarithmisch. 


gemessenen LSVen sind in den Abb./4.6b|und Abb. /4.6c|dargestellt. Die 
Gradienten-Kurven der gemessenen LSVen passen sehr gut zu denen 
der rückgeführten Flächen. Die Abweichung der Gradienten-Kurven zu 


der nominellen Simulation kann durch eine nicht korrekte Abformung 
der Struktur beim Spritzgussprozess oder eine zu geringe Strukturhöhe 


im Werkzeugeinsatz verursacht sein. In Abb.|4.8)werden nominelle und 
gemessene Strukturhöhe des Werkzeugeinsatzes 1 verglichen. Die Struk- 
turhöhe entspricht in der Nähe des Zenits sehr gut den Vorgaben. In 
der Nähe der horizontalen Beschnittkanten ist die gemessene Struktur- 
höhe geringer als vorgegeben. Durch die geringere Strukturhöhe wird 
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(a) Nominelle Simulation. 
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(b) Werkzeugeinsatz 1. 


— Nominell ——Rückf. 
I, (cd) = Abgeformte Linsen 
10! 10? 10? 10* G(B) 
0.4 
0.2 
-1 
= | = a(°) 0.0 o 
3 -2 -1 0 -1.0 =05 oof) 


(c) Werkzeugeinsatz 2. 


Abbildung 4.6: Simulation der LSV und der Gradienten-Kurve des BiLED-Moduls mit 
einer strukturierten Linse ausgelegt für den ECE-Bereich im Vergleich zu gemessenen 
Gradienten-Kurven von abgeformten Linsen. Die gestrichelte Linie in den LSVen gibt die 
Position der Gradienten-Kurve an. 
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10 10-7 10° 10° 10! 
Abweichung (um) 


Abbildung 4.7: Abweichung zwischen Flächenrückführung und Messdaten für die zwei 
Werkzeugeinsätze mit Struktur für ECE-Gebiete. Links: Werkzeugeinsatz Nr. 1, Rechts: 
Werkzeugeinsatz Nr. 2. 


weniger Licht gestreut, sodass der Gradient einen größeren Wert hat. 
Ähnliche Unterschiede zeigen sich auch bei dem Werkzeugeinsatz 2. Die 
Abb./4.9|zeigt ein Foto einer strukturierten Linse. 
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z (um) 
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Abbildung 4.8: Vergleich der Strukturhöhen zwischen den nominellen und den gemesse- 
nen Daten für Werkzeugeinsatz 1. 


Abbildung 4.9: Foto einer Linse mit der ECE-Struktur. 
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4.3 FÜR FMVSS-GEBIETE 


Die für den FMVSS-Bereich ausgelegte Mikrostruktur ist in Abb.|4.10]dar- 
gestellt. Die Struktur basiert auf einem Gitter mit einer Auflösung von 
100 pm. Wie bei der ECE-Struktur ist die Strukturierung auf der Linsen- 
fläche inhomogen. Zwei Bereiche der Linse sind nicht strukturiert, um 


y (mm) 


x (mm) 


-10 -5 2 -1 -0.5 -0.25 0 025 05 1 2 5 10 
z (um) 


Abbildung 4.10: Nominelle Mikrostruktur zur Aufweichung der HDG für den FMVSS- 
Bereich. Die Farbskala ist unter —0,5 um und über 0,5 um logarithmisch. 


die Beeinflussung des Fernlichtes gering zu halten. In den Bereichen mit 
Struktur ist die Amplitude der Struktur größer als bei der ECE-Struktur. 
Die größere Amplitude ist notwendig, um einen geringeren Gradient 
mit einem Wert von 0,2 zu erreichen. Die nominelle simulierte LSV und 
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die Gradienten-Kurve sind in Abb.!4.11a dargestellt. Analog zur Herstel- 
lung der Linsen mit ECE-Struktur, werden zwei Werkzeugeinsätze (Nr. 


3 & 4) gefertigt und vermessen. Das Ergebnis der Flächenrückführung 
bzw. die Abweichung zwischen Punktgitter und Fläche sind in Abb. 


4.12 dargestellt. Für die Werkzeugeinsätze 3 und 4 beträgt die mittlere 


orthogonale Abweichung zwischen Fläche und Punktgitter 35 nm und 
37 nm. Im Bereich des mittleren Streifens der Struktur zeigen sich lokal 
größerer Abweichungen als bei der ECE-Struktur, wobei die Strukturhö- 
he selbst größer ist. In Abb. 4.11blund Abb./4.11c|sind die simulierten 
LSVen und Gradienten-Kurven auf Basis der rückgeführten Flächen 


dargestellt. In den gleichen Abbildungen sind die Gradienten-Kurven 
von gemessenen LSVen basierend auf abgeformten Linsen aufgezeigt. 
Die Gradienten-Kurven der abgeformten Linsen passen sehr gut zu der 
nominellen Simulation. Die Gradienten-Kurve der Flächenrückführun- 
gen weichen mit einem etwas größeren Gradienten ab. In Abb. 4.13]sind 
die nominellen und gemessenen Strukturhöhen des Werkzeugeinsat- 
zes 3 verglichen. Im Bereich um den Zenit der Linse wird die nominelle 
Strukturvorgabe gut erfüllt. In dem Bereich nahe der Schnittkante der 
Linse zeigen sich starke Abweichungen, bzw. eine Korrelation zwischen 
Soll- und Ist-Struktur ist kaum noch feststellbar. Da die Strukturhöhe 
am Randbereich um ca. einen Faktor 10 geringer ist als im Inneren, sind 
die Auswirkungen der Abweichung auf die Gradienten-Kurve gering. 
Die abgeformten Linsen zeigen einen Einfall in der Lichteintrittsseite 
der Linse, visualisiert in Abb. Der Einfall hat eine defokussieren- 
de Wirkung auf das Licht, wodurch der Gradient der HDG verringert 
werden kann. Durch eine zu geringe Strukturhöhe, wird der Gradient 
größer, wie in den simulierten Gradienten-Kurven der rückgeführten Flä- 
chen dargestellt. Der Einfall der Lichteintrittsseite kann der zu geringen 
Strukturhöhe durch die Defokussierung entgegenwirken, weshalb die 
Gradienten-Kurven der abgeformten Linsen mit denen der nominellen 


Simulation übereinstimmen. 


81 


ANWENDUNG DER PROZESSKETTE 


0.1 


0.0 | — 
-1.0 -0.5 0.0 8 C) 


-3 -2 -1 0 40 


(a) Nominelle Simulation. 
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(b) Werkzeugeinsatz 3. 
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(c) Werkzeugeinsatz 4. 


Abbildung 4.11: Simulation der LSV und der Gradienten-Kurve des BiLED-Moduls mit 
einer strukturierten Linse, ausgelegt für den FMVSS-Bereich im Vergleich zu gemessenen 
Gradienten-Kurven von abgeformten Linsen. Die gestrichelte Linie in den LSVen gibt die 
Position der Gradienten-Kurve an. 
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10-7 10 10-! 10° 10! 
Abweichung (um) 
Abbildung 4.12: Abweichung zwischen Flächenrückführung und Messdaten für die zwei 


Werkzeugeinsätze mit Struktur für FMVSS-Gebiete. Links: Werkzeugeinsatz Nr. 3, Rechts: 
Werkzeugeinsatz Nr. 4. 


| — Nominell —— Rückf. Messdaten 


z (um) 


Abbildung 4.13: Vergleich der Strukturhöhen zwischen den nominellen und den gemesse- 
nen Daten für Werkzeugeinsatz 3. 
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Abbildung 4.14: Foto einer Fresnel-Reflexion einer Deckenleuchte an der Lichteintrittsseite 
der Linse. Die Verzerrung der rechteckigen Leuchte entsteht durch den Einfall der nominell 
planaren Fläche. 
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KAPITEL 5 


OPTIMIERUNG DER AUSLEGUNG 
DER MIKROSTRUKTUREN 


Die Auslegung der Strukturen bisher, einschließlich denen aus Kapitel 


4 erfolgt in manuellen Korrekturschleifen. Ein Optikingenieur erstellt 
eine Struktur, führt eine lichttechnische Simulation durch und bewertet 
die Gradienten-Kurve. Je nach resultierender Gradienten-Kurve wird 
z. B. die Amplitude der Struktur vergrößert oder verringert. Um eine 
Struktur in angemessener Zeit auszulegen, wird viel Erfahrung benötigt. 
Dieses Kapitel befasst sich mit der Reduzeriung der Simulationszeit von 
Scheinwefersystemen mit strukturierten Linsen und der Auslegung von 
solchen Strukturen durch Anwendung von Optimierungsalgorithmen. 
Der Inhalt dieses Kapitels stammt zu großen Teilen aus der Veröffentli- 


chung [93]. 


5.1 REDUZIERUNG DER SIMULATIONSSZEIT 


Um eine auswertbare, also nicht zu verrauschte Gradienten-Kurve zu 
erhalten, muss mit mehreren hundert Millionen Strahlen simuliert wer- 
den. Eine Simulation mit 500 Millionen Strahlen des BiLED-Moduls (vgl. 
Abb. dauerte auf dem Test-Computer mit einem Intel Xeon E5-1650 
v3 ca. eine Stunde [56]. Für einen Optimierungsprozess, bei dem stän- 
dig neue Strukturen generiert und simuliert werden, ist eine Stunde 
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Linse: Licht- 
Reflektor eintrittsseite 


LED 
Strahlen 
Blende Detektor 


Abbildung 5.1: Schematische Darstellung des Aufbaus zur Vorbereitung des Optimie- 
rungsprozesses 


zu lang. Zum Beispiel würde die Simulation von 1000 verschiedenen 
Strukturen ca. 40 Tage dauern. Zur Reduzierung der Simulationszeit 


wird der Aufbau aus Abb. |5.1/verwendet. Der Aufbau entspricht dem 


kompletten Scheinwerfermodul, ohne die Lichtaustrittsfläche der Linse. 
Ein Strahlendetektor ist hinter der Lichteintrittsseite der Linse platziert. 
Der Detektor speichert die Informationen der Strahlen die ihn treffen. Zu 
den Informationen gehören Schnittpunkt, Richtungsvektor, Lichtstrom 
und Wellenlänge. Die abgespeicherten Strahlen werden als Lichtquel- 
le für den zweiten Schritt der Vorbereitung verwendet. Die Strahlen- 
Lichtquelle modelliert das Scheinwerfermodul ohne Lichtaustrittsseite 
der Linse. 


Bei der Simulation des BiLED-Moduls mit 500 Millionen Strahlen er- 
reichen ca. 400 Millionen Strahlen den Detektor. Wird die Strahlen- 
Lichtquelle zur Simulation des Moduls verwendet und nur der Einfluss 
der Lichtaustrittsseite simuliert, so reduziert sich die Simulationszeit mit 
den 400 Millionen Strahlen auf ca. 10 Minuten. Im zweiten Schritt wer- 
den die für die HDG relevanten Strahlen herausgefiltert. Dazu wird eine 
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Freiformfläche 
ohne Mikrostruktur 


Strahlen” 
Detektor 
(grün) 


Lichtquelle a (°) 
(gelb) N 


Lichtstärkeverteilung 


Abbildung 5.2: Filterung der Strahlen: Es werden nur die Strahlen behalten, die in einen 
Bereich in der Nähe der HDG treffen (weißes Rechteck). 


Simulation mit der unstrukturierten Linsenfläche durchgeführt. Es wird 
die Linsenfläche verwendet, die die gleiche Kontrollpunktdichte hat wie 


die später erstellte strukturierte Fläche (vgl. Kap|3.2.1|,,zweites Verfah- 


ren”). Dadurch haben unstrukturierte und strukturierte Linsenflächen 
eine ähnliche Parametrisierung. Werden die Strahlen durch die unstruk- 
turierte Linsenfläche in einen vorher definierten Bereich in der Nähe der 
HDG gebrochen (vgl. Abb. 5.2), so werden die Strahl-Informationen und 
die uv-Koordinaten der Schnittpunkte gespeichert. Das Ergebnis der 


Filterung ist eine Strahlen-Lichtquelle, die nur die Strahlen enthält, die 
die HDG erzeugen bzw. die Gradienten-Kurve beeinflussen. Ein Versuch 
zeigt, dass nach der Filterung ca. 10 Millionen Strahlen übrig bleiben. Die 
Simulation des BiLED-Moduls mit der gefilterten Strahlen-Lichtquelle 
benötigt ca. 18s. 


Bei einer Standard-Simulation wird zu jeder Zeit die Information ge- 
speichert, in welchem Material sich der jeweilige Strahl befindet. Durch 
numerische Ungenauigkeiten kann es passieren, dass ein Strahl, der sich 
in einem Material A befindet auf eine Flache trifft, die die Materialien B 
und C voneinander trennt. Dieser Fall stellt einen Simulationsfehler da, 


87 


OPTIMIERUNG DER AUSLEGUNG DER MIKROSTRUKTUREN 


sodass die weitere Verfolgung des Strahls verworfen wird. Da nur eine 
optische Fläche in dem Simulationsaufbau während der Optimierung 
vorhanden ist, kann auf die Überprüfung verzichtet werden. Um die 
Schnittpunkte bei einer Simulation schnell berechnen zu können, wird 
vor dem Ray-Tracing eine Hierarchie von Begrenzungsvolumen aufge- 
baut (vgl. Kap. P-1.2). Die Berechnung der Volumen dauern für eine 
strukturierte Fläche einige Sekunden. Wie oben beschrieben, werden bei 
der Filterung die uv-Koordinaten der Schnittpunkte gespeichert. Anstatt 
die Begrenzungsrahmen zu verwenden, wird bei der Optimierung direkt 
auf der Fläche mit dem Newton-Verfahren (vgl. Kap. /2.1.2) nach den 
Schnittpunkten gesucht. Die uv-Koordinaten werden als Startpunkte für 


die Iteration verwendet. 

Bei dem Newton-Verfahren müssen die B-Spline-Flächen und deren 
Ableitungen ausgewertet werden. Während der Optimierung werden 
dafür die SIMD-Implementierungen (vgl. Kap. 2.2.4) verwendet. Durch 
die Verwendung des „zweiten Verfahrens “ aus Kap. 3.2.1|zur Strukturie- 


rung der Flächen reduziert sich die Strukturierung auf eine Verschiebung 
der Kontrollpunkte. Es wird kein Gleichungssystem gelöst, wodurch 
Zeit eingespart wird. Bei der Optimierung müssen kontinuierlich Struk- 
turen erstellt, simuliert und die Gradienten-Kurven bewertet werden. 
Diese drei Schritte werden zwölf-fach parallel ausgeführt, sodass auf 
jedem logischen Prozessorkern je eine Struktur bewertet wird. Durch 
die genannten Maßnahmen können, unter Verwendung einer Strahlen- 
Lichtquelle mit 10 Millionen Strahlen, zwölf Strukturen in 24s erstellt 
und bewertet werden. Das entspricht umgerechnet 2s pro Struktur. 


5.2 GENETISCHE ALGORITHMEN 


Um eine Optimierung durch Algorithmen zu ermöglichen, muss de- 
finiert werden, was eine gute von einer schlechten Mikrostruktur un- 
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terscheidet. Dazu wird eine Soll-Gradienten-Kurve definiert. In dieser 
Arbeit werden die simulierten Gradienten-Kurven der Strukturen, aus- 
gelegt durch Optikingenieure, als Soll-Kurven verwendet. Als Qualitäts- 
maß wird die quadratische Abweichung zwischen der vorgegebenen 
und den simulierten Gradienten-Kurven verwendet: 


N 
p=- 3 (Gist(Bi) — Gson(Bi))” (5.1) 

i= 
Da das Ergebnis einer Simulation eine diskretisierte Gradienten-Kur- 
ve ist, gibt N die Anzahl an Stiitzstellen an. Im Zusammenhang mit 
genetischen Algorithmen beschreibt Gl. die Qualitäts- bzw. Fitness- 
Funktion, die tiberlicherweise maximiert wird. Daher wird als Qua- 
litätsmaß das Negative der Summe aus Gl. verwendet. Ziel der 
Optimierung ist es, eine Struktur zu finden, die ® maximiert. Die ideale 
Struktur würde ® = 0 erfüllen. Je nach Anzahl an Parametern, die bei 
der Optimierung verwendet werden, sind auf Grund der stochastischen 
Strukturen, lokale Minima in der Funktion ® zu erwarten. Die Ableitung 
der Qualitätsfunktion nach den Parametern der Mikrostruktur (z. 
B. Gitterabstand oder Amplitude) ist schwierig zu bestimmen. Daher 
wird ein Algorithmus benötigt, der global nach Minima sucht und un- 
abhängig von den Ableitungen der Qualitätsfunktion ist. Sogenannte 


genetische Algorithmen erfüllen diese Kriterien 94 95]. 


Die genetischen Algorithmen basieren auf Phänomenen, wie sie in der 
Natur vorkommen. Dazu gehören die natürliche Selektion und die gene- 


tische Evolution. In Abb.|5.3\ist schematisch der Ablauf eines genetischen 


Algorithmus dargestellt. Zu Beginn wird eine Sammlung (Population) 
von verschiedenen mikrostrukturierten Linsenflächen (Individuen) be- 
stimmt, deren Parameter (Gene) zufällig gewählt werden. Es wird der 
Fitness-Wert ® jedes Individuums bestimmt. Für die Mikrostrukturen 
wird dazu die Gl. ausgewertet. Die Fitness wird so interpretiert, 
dass ein größeres ® eine bessere Fitness darstellt. Wird ein Abbruch- 
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Abbruch 
Kriterium 
erreicht? 


Zufällige Bewertung 
Population der Fitness 


Selektion 


Rekombination 


Mutation 


Neue Population 


Abbildung 5.3: Schematische Darstellung des Ablaufes eines genetischen Algorithmus. 


kriterium erfüllt, so endet der Algorithmus. Abbruchkriterien können 
eine maximale Anzahl an Iterationen (Generationen) oder ein erreichter 
Fitness-Wert sein. 


Ist kein Abbruchkriterium erfüllt, so wird eine neue Population be- 
stimmt. Zur Bestimmung der neuen Population werden die Individuen, 
ihrem Fitness-Wert entsprechend, sortiert. Eine vorher definierte Anzahl 
an Individuen mit den besten Fitness-Werten werden unverändert in 
die neue Population übernommen (sog. Elitismus). Die restlichen neuen 
Individuen werden auf Basis der bestehenden Population bestimmt. Für 
die sog. Turnier-Auswahl (Selektion) werden zwei Individuen zufällig 
aus der bestehenden Population, einschließlich der bereits übernomme- 
nen Elite, gewählt [96]. Die Fitness beider Individuen wird verglichen. 


Das Individuum mit der besseren Fitness „gewinnt“ das Turnier, wird 
aus der alten Population entfernt und wird für die Rekombination ver- 
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wendet. Über eine zweite Turnier-Auswahl wird ein weiteres Individu- 
um für die Rekombination bestimmt. Bei der Rekombination werden 
zwei neue Individuen (Nachkommen) aus den zwei zuvor gewählten 
Individuen bestimmt. Hier wird die sog. „Intermediate Recombinati- 
on“ angewendet [96]. Bei dieser Methode werden die neuen Parameter 


durch gewichtete Linearkombination der alten Parameter bestimmt. 
Seien 7 = Jun, v1, ..-, ul und © = {wo,wı,...,wj} die Parameter der 
beiden gewählten Mikrostrukturen (Individuen). Bei der Rekombination 
werden die neuen Parameter wie folgt bestimmt: 


Programmauflistung 5.1: „Intermediate Recombination” [96] 


1 for i+ 0,1 

2 do 

3 a <randomValue(—p, 1+ p) 
4 b + randomValue(—p, 1+ p) 
5 t 4 avi + (1 — a)wi 

6 s¢ bw; + (1 — bi: 

7 while outsideBounds(t,s) 

8 


Oi = i 
9 Wj — Ss 
10 end for 


Es werden zwei Zufallszahlen zwischen —p und 1+ p bestimmt. Durch 
lineare Kombination werden die neuen Parameter t und s bestimmt. 
Sind die Parameter außerhalb eventuell gesetzter Schranken, wird der 
Vorgang wiederholt, bis die neuen Parameter innerhalb der Schranken 
liegen. Wird p = 0 verwendet, werden die neuen Parameter nur inner- 
halb des Intervalls [v;, wi] der alten Parameter bestimmt. Dadurch kann 
eine vorzeitige Konvergenz in Bereichen lokaler Minima stattfinden. Um 
das zu verhindern, wird p = 0,5 verwendet, sodass die neuen Parameter 
auch außerhalb des Intervalls [v;, w;] gewählt werden können. 
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Die Parameter der neuen Mikrostrukturen (Nachkommen) werden mu- 
tiert. Durch Mutation wird eine vorzeitige Konvergenz des Algorithmus 
zusätzich verhindert und die Suche nach globalen Minima ermöglicht. 
Mutation bedeutet, dass zufällig ausgewählte Parameter eines Individu- 
ums um zufällig gewählte Werte verschoben werden [97]. Die Mutation 
die hier verwendet wird, wählt von jeder Mikrostruktur einen Parameter 


aus und verschiebt diesen durch 
Wj — Wj TET (Wmaxi = Wmin,i) € (5.2) 


Die Symbole Wmax,j und Wminj geben die Schranken für den Parameter 
w; an. Das Symbol r gibt den Mutationsbereich an. Es wird r = 0,1 
verwendet [97]. Das Vorzeichen des Versatzes ist durch s gegeben, wobei 
s € {-1,1} zufällig bestimmt wird. Sollte das neue w; außerhalb der 
Schranken liegen, so wird w; auf den Wert der jeweiligen Schranke 
gesetzt. Die Schritte Selektion, Rekombination und Mutation werden 
wiederholt, bis die neue Population mit Individuen gefüllt ist. Es werden 
die Fitness-Werte der neuen Population bestimmt und das Abbruch- 
Kriterium überprüft. Die Iteration wird durchgeführt, bis das Abbruch- 
Kriterium erreicht ist. 


5.3 STOCHASTISCHE STRUKTUREN 
5.3.1 BILED-MODUL 


Zum Testen der Optimierung wird das BiLED-Modul verwendet. Es 


wird die erste Methode aus Kap.|3.2.1)zur Strukturierung angewendet. 


Die Mikrostruktur basiert auf einem Punktgitter mit einem Abstand von 
100 um. Die Amplituden werden aus Zellen der Größe 2x2 gezogen (vgl. 
Kap. 3.2.3). Die Einträge der Hauptdiagonalen sind die Parameter, die 
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— Bel-ı‘, 0.5°] 
— pe[-15*% 1°] 
—— PEl-1.25°, 13] 
--- Soll 


Abbildung 5.4: Vergleich zwischen der Soll- und den resultierenden Gradienten-Kurven 
fiir die drei Optimierungen mit unterschiedlich gefilterten Strahlen-Lichtquellen. 


der Algorithmus optimieren soll. Die Eintrage der Nebendiagonalen 
werden auf null gesetzt, also: 


T= & a (5.3) 


Bei der Optimierung beträgt die Populationsgröße 36, wobei eine Elite, 
die 12 Individuen umfasste, verwendet wird. Die Soll-Gradienten-Kurve 


ist in Abb. 5.4\ dargestellt und entspricht der simulierten Gradienten- 
Kurve des BiLED-Moduls mit ECE-Struktur. Es werden drei Optimie- 
rungen durchgeführt, um die Auswirkungen des Filterbereiches (vgl. 
Kap.5.lund Abb.[5.2) und der Strahlanzahl zu untersuchen. Die Filter- 
bereiche in vertikaler Richtung sind [—1°, 0,5°] mit 6 Millionen Strahlen, 
[-1,25°, 1°] mit 8 Millionen Strahlen und [—1,5°, 1.5°] mit 10 Millionen 
Strahlen. Der horizontale Bereich ist [—5°, 1°]. Die Auflösung des Detek- 
tors, der in der Simulation verwendet wird, beträgt 0,05° vertikal und 


0,25° horizontal. Die Auflösung in horizontaler Richtung wird niedrig 
gewählt, um das Simulationsrauschen zu verringern. Das ist möglich, 
da benachbarte horizontale Zellen in der Nähe der HDG einen sehr 
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ähnlichen Wert haben. Die nach 20 Generationen erzielten Gradienten- 
Kurven der jeweiligen Strukturen mit der maximalen Fitness sind in 
Abb.|5.4|dargestellt. Alle drei Kurven passen gut zur Zielvorgabe. Die 
Abweichung der blauen Kurve in dem Bereich f > 0° entsteht durch 


die gefilterte Strahlen-Lichtquelle. Bei der blauen Kurve werden nur die 
Strahlen verwendet, die auf einen Bereich unterhalb von 0,5° auftreffen. 
Da die Struktur das Licht auch in negativer Richtung ablenkt, kann die 
Gradienten-Kurve in diesem Bereich nicht korrekt simuliert werden. In 


Abb.\5.5list die maximale und die durchschnittliche Fitness pro Gene- 


ration aufgetragen. Die Kurven der maximalen Fitness erreichten einen 


— BEel-1?, 0.5°] — BeE[-1.5°, 1°] — BeE[—-1.25°, 1°] 
-0.01 0.0 
-0.02 -05 
5 
£ -0.05 le -1.0 
Ki 
-0.1 
-1.5 
-0.2 
2.0 
0 10 20 0 10 20 
Nr. Generation Nr. Generation 


Abbildung 5.5: Konvergenzverhalten des genetischen Algorithmus bei der Optimierung 
von stochastischen Strukturen für das BiLED-Modul. Links: Maximale Fitness pro Genera- 
tion. Rechts: Durchschnittliche Fitness pro Generation. 


Wert zwischen —0,014 bis —0, 012. Die Kurven der durchschnittlichen 
Fitness erreichten einen Wert von ca. —0,1. Die Gradienten-Kurven 
sind durch die Monte Carlo-Simulation immer verrauscht, weshalb eine 
exakte Anpassung sehr unwahrscheinlich ist. Der Unterschied zwischen 
maximaler und durchschnittlicher Fitness bedeutet, dass die Population 
noch nicht komplett konvergiert ist. Da die Gradienten-Kurven jedoch 
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sehr gut zu den Vorgaben passen, hätten weitere Iterationen keine 
merkliche Verbesserung mehr gebracht. Die simulierten LSVen sind in 
der Abb. 5.6 dargestellt. Die Farbskala zu den LSVen ist in willkürlichen 


BE[-1°, 0.5°] 


2 
I 
I 
I, (wilk. Einh.) 
M oor 
10° 10! 10? 10? 


BE[-1.5°, 1.5°] BE [-1.25°, 1°] 


2 2 
I 
Z 0 z 0 
a a 
_2 2 


Abbildung 5.6: Simulierte LSVen für das BiLED-Modul mit jeweils der stochastischen 
Struktur mit der maximalen Fitness. Die gestrichelte Linie zeigt die Position, an der die 
Gradienten-Kurve bestimmt wurde. 


Einheiten angegeben, da bei der Generierung der gefilterten Strahlen- 
Lichtquelle ein Skalierungsfaktor verloren geht. Da die Berechnung 
der Gradienten-Kurve invariant bezüglich der Skalierung der LSV ist, 
werden keine absoluten Lichtstärkewerte benötigt. Die Parameter der 
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Abbildung 5.7: Simulation des kompletten BiLED-Moduls mit einer Mikrostruktur erstellt 
mit den Parametern oz 1. Die gestrichelte Linie markiert die Position zur Auswertung 
der Gradienten-Kurve. Rechts: Vergleich zwischen Soll- und Ist-Gradienten-Kurve. 


Strukturen mit dem besten Fitness-Wert sind: 


0,67 0 0,66 0 
aa) og Ber Sg Ro 


0,6 0 
01-1.25,1] = 0 0,16 pm 


Die Parameter der Strukturen sind nahezu identisch. Da die Opti- 
mierung nur auf zwei Parameter basiert, war eine Ähnlichkeit der 
Ergebnisse zu erwarten. Bei zwei Parametern ist die Dimension des 
Parameterraums so klein, dass die Wahrscheinlichkeit, das der Opti- 
mierungsalgorithmus bei unterschiedlichen Läufen unterschiedliche 
Minima findet, sehr gering ist. Die Simulation der Mikrostruktur mit 
dem kompletten Scheinwerfer-Modul ist in Abb.|5.7|dargestellt. Die Si- 
mulation des kompletten Moduls reproduziert sehr gut die Gradienten- 


Kurve der Optimierung. Dieser Versuch zeigt, dass Mikrostrukturen 
durch Optimierungsalgorithmen ausgelegt werden können. Bei der Aus- 
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LED-Lichtquellen Primäroptik Linse 


Abbildung 5.8: Aufbau des MLS-BiLED-Scheinwerfermoduls. 


wahl des Filterbereiches zur Generierung der Strahlen-Lichtquelle muss 
darauf geachtet werden, dass der Bereich um ca. 0,5° größer ist als der 
Bereich der Soll-Gradienten-Kurve, wie Abb.|5.4|zeigt. Bei einer flache- 
ren Soll-Gradienten-Kurve müsste der Bereich eventuell größer gewählt 
werden, da bei einer flacheren Kurve mehr Streuung benötigt wird. 


5.3.2 MLS-BILED-MODUL 


Für einen zweiten Test wird das MLS-BiLED-Scheinwerfermodul ver- 
wendet. Die Abb. |5.8| zeigt den Aufbau des Moduls. Für das Modul 
wurde eine Mikrostruktur von einem Optikingenieur unter Anwen- 
dung der alten Prozesskette (vgl. Abb. 1.3) ausgelegt. Die simulierte 
LSV und Gradienten-Kurve ist in Abb.|5.9\dargestellt. Durch den gene- 
tischen Algorithmus wird versucht, eine Mikrostruktur zu finden. Die 
Mikrostruktur soll eine Gradienten-Kurve erzeugen, die der Nominel- 
len, vom Optikingenieur ausgelegten, möglichst ähnlich ist. Dazu wird 
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Abbildung 5.9: Simulation des kompletten MLS-BiLED-Moduls mit einer Mikrostruktur 
ausgelegt durch einen Optikingenieur. Rechts: Vergleich zwischen der simulierten (grün) 
und der für die Optimierung verwendeten Soll-Gradienten-Kurve (rot). 


eine Soll-Gradienten-Kurve definiert (rote Kurve in Abb. 5.9), die der 
nominellen geglatteten Kurve entspricht. Im Vergleich zu der Struktur 
fiir das BiLED-Modul ist der Gradient geringer und die Gradienten- 
Kurve fällt erst bei höheren Winkeln auf Null ab. Durch die flachere 
Kurve wirkt sich das Simulationsrauschen stärker aus. Daher werden 


zwei Optimierungen durchgeführt. Bei der ersten Optimierung werden 
10 Millionen Strahlen in dem Bereich [—1.5°, 1.5°] und bei der zweiten 
Optimierung 40 Millionen Strahlen in dem Bereich [—2.5°, 2.5°] verwen- 
det. Die Größe der Population beträgt 36 Individuen, die Größe der Elite 
beträgt 12 Individuen. Die erzielten Gradienten-Kurven sind in Abb. 
5.10]dargestellt. Beide Kurven passen gut zur Zielvorgabe. Die blaue 
Kurve, basierend auf der Simulation mit 10 Millionen Strahlen, zeigt ein 
größeres Simulationsrauschen auf Grund der geringeren Strahlanzahl. 
Das Rauschen zeigt sich auch in einer etwas schlechteren Fitness (vgl. 
Abb. 5.11). Wie bei dem BiLED-Modul, ist die erreichte durchschnittliche 
Fitness der Population um ungefähr einen Faktor 10 schlechter als die 
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Abbildung 5.10: Vergleich zwischen der Soll- und den resultierenden Gradienten-Kurven 


für die zwei Optimierungen mit unterschiedlich gefilterten Strahlen-Lichtquellen. 


maximale Fitness. Weitere Iterationen hätten keine merkliche Verbes- 
serung gebracht, da Soll- und Ergebniskurve sehr ähnlich sind. In den 
Abb. sind die simulierten LSVen dargestellt. 


Die Parameter der Strukturen sind: 


104 0 1,03 0 
7-15, 1.5] = 0 0,19 mm O[-25, 2.5] = 0 0,32 pm. 


Die Simulation des kompletten Moduls mit der Struktur mit den Para- 


metern o a 2.5] ist in Abb.|5.13]dargestellt. Die Gradienten-Kurve der 


simulierten LSV ist stärker verrauscht als die der Optimierungen, da hier 
eine horizontale Detektorauflösung von 0,05° verwendet wird. Trotz des 
Rauschens ist erkennbar, dass simulierte und vorgegebene Gradienten- 
Kurven gut übereinstimmen. 


Die Auslegung der Strukturen für das MLS-BiLED-Modul bestätigen, 
dass eine Auslegung durch Optimierungsalgorithmen möglich ist. Beide 
Versuche verwendeten zwei Optimierungsparameter. Durch die geringe 
Anzahl an Parametern und die gesetzten Schranken wäre der Einsatz 


99 


OPTIMIERUNG DER AUSLEGUNG DER MIKROSTRUKTUREN 


— BE[-1.5°, 1.5°] — BEl[-2.5°, 2.5°] 
-0.01 0.0 
-0.02 
—0.5 
5 
E e 
£ -0.05 l 
—1.0 
-0.1 
-0.2 —1.5 
0 20 40 0 20 40 
Nr. Generation Nr. Generation 


Abbildung 5.11: Konvergenzverhalten des genetischen Algorithmus bei der Optimierung 
von stochastischen Strukturen für das MLS-BiLED-Modul. Links: Maximale Fitness pro 
Generation. Rechts: Durchschnittliche Fitness pro Generation. 


eines genetischen Algorithmus eventuell nicht nötig gewesen. Vermut- 
lich hätten durch ein lineares Absuchen des Parameterraums geeignete 
Strukturen gefunden werden können. Die Parameter, die zur Optimie- 
rung verwendet werden, stellen eine Art Streumaß dar. Je größer die 
Parameter gewählt werden, umso stärker wird die Flächennormale der 
Linse um ihre ursprüngliche Ausrichtung verkippt. Es findet eine zufäl- 
lige Streuung des Lichtes statt, wodurch der Gradient verringert wird. 
In weiteren Versuchen wird nach Strukturen gesucht, die keine zufäl- 
lige Verkippung der Flächennormalen haben, sondern bei denen die 
Verkippung direkter der Soll-Gradienten-Kurve angepasst werden. 
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Abbildung 5.12: Simulierte LSVen für das MLS-BiLED-Modul mit jeweils der stochasti- 
schen Struktur mit der maximalen Fitness. Die gestrichelte Linie zeigt die Position, an der 


die Gradienten-Kurve bestimmt wurde. 
0.25 
0.20 
0. 003 


ly (cd) B o 


Abbildung 5.13: Simulation des kompletten MLS-BiLED-Moduls mit einer Mikrostruktur 
erstellt mit den Parametern o ——— 2.50]. Rechts: Vergleich zwischen Ziel- und Ist-Gradien- 
ten-Kurve. 
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5.4 PERIODISCHE FREIFORM-STRUKTUREN 


Bei den folgenden Versuchen wird die zweite Methode (Verschiebung 
der Kontrollpunkte, vgl. Kap. 3.2.1) angewendet, um für das BiLED- 
Modul geeignete nicht-stochastische Strukturen zu finden. Anstatt das 


Punktgitter in Zellen einzuteilen, wird die Kontrollpunktmatrix der zu 
strukturierenden Fläche eingeteilt. Bei einem Versuch beträgt die Grö- 
ße der Zellen 1x10 (horizontal x vertikal), bei einem zweitem Versuch 
2x6. Die Parameter, die dem Optimierungsalgorithmus zur Verfügung 
stehen, sind Verschiebungen der Kontrollpunkte in vertikaler (y) und 
longitudinaler (z) Richtung. Die Verschiebung in y-Richtung gibt die 
Möglichkeit, eine Struktur mit steileren Flanken (größere Ablenkung) zu 
generieren, ohne die Amplitude der Strukturen vergrößern zu müssen. 
Es wird eine Fläche verwendet, deren Kontrollpunkte einen Abstand 
von 200 um in vertikaler und horizontaler Richtung haben. Der Opti- 
mierungsalgorithmus darf die Kontrollpunkte in vertikaler Richtung 


um +95 um und in longitudinaler Richtung um +1 pm verschieben. Es 
wird eine Strahlenlichtquelle mit 10 Millionen Strahlen aus dem Bereich 
[-3°, 1°] horizontal und [—5°, 1°] vertikal verwendet. Die Größe der 
Population beträgt 48 und die Größe der Elite 24. Die Gradienten-Kur- 
ven der jeweiligen Mikrostruktur mit dem höchsten Fitness-Wert sind 


in Abb. 5.14 dargestellt. Die vorgegebene Gradienten-Kurve wird sehr 


gut reproduziert. Auf Grund der höheren Anzahl an Parametern, wird 
eine Konvergenz erst nach ca. 200 Generationen erreicht (vgl. Abb.B.15). 
In Abb 5.16|sind die simulierten LSVen dargestellt. Im Vergleich zu den 
simulierten LSVen der stochastischen Struktur (vgl. Abb. 5.6) fallt auf, 
dass die Lichtstärke in dem Bereich (x > 0°,ß > 0°) größer ist. Die 
Parameter, mit denen die Strukturen erstellt worden sind, sind: 


APyixio = (2,9; —2,7; —13,6; —1,4; —3,1; 18,9; —19,2; —14,7; 
7,1; =1,7)" um 
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Abbildung 5.14: Vergleich zwischen der Soll- und den resultierenden Gradienten-Kurven 
fiir die zwei Optimierungen mit unterschiedlichen Freiheitsgraden der Mikrostrukturen. 


AP; ıxıo = (— 0,31; 0,58; 0,49; 0,41; 0,47; —0,10; —0,18; 0,14; 
— 0,22; 0,81)" pm 


an, — 67 0,05 —0,14 
44,2 —22,6 0,82 0,46 
<7, 9° 250-9 0,65 —0,60 
AP = fs 7 AP = D z 
pare 1 CA Keep | a 
41,8 —42,3 0,59 0,89 
10,2 34,6 0,24 —0,36 


Die simulierten LSVen und Gradienten-Kurven des kompletten Moduls 
sind in Abb. 15.17 dargestellt. Die Gradienten-Kurven passen sehr gut zu 
den Vorgaben. Die simulierten LSVen zeigen wellenförmige Modulatio- 
nen in dem Bereich 8 < —10°. Die Modulationen sind bei der Struktur 
mit der Zellengröße 1x10 stärker ausgeprägt als bei der Struktur mit der 
Zellengröße 2x6. Die Verwendung von Zellen, die in einer Richtung eine 
Größe von 1 haben, erzeugen eine eindimensionale Struktur, die das 
Licht nur in vertikaler Richtung ablenkt (vgl. Abb. 5-17). Zweidimensio- 
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Abbildung 5.15: Konvergenzverhalten des genetischen Algorithmus bei der Optimierung 
von periodischen Freiform-Strukturen. Links: Maximale Fitness pro Generation. Rechts: 
Durchschnittliche Fitness pro Generation. 


nale Zellen verursachen eine zweidimensionale Streuung. Dadurch wird 
die Ausprägung der Modulation verringert. 

Wie bei den Simulationen der Optimierung zeigt sich auch bei der Si- 
mulation des gesamten Moduls eine höhere Lichtstärke in dem Bereich 
(a > 0°, B > 0°). Die Lichtstärke übersteigt in diesem Bereich teilweise 


die gesetzlichen Höchstgrenzen. In Tab. 5.1|sind die Werte der Mess- 


punkte aufgezeigt, die überschritten werden. Die Tabelle vergleicht die 


gesetzlichen Höchstwerte mit den Simulationen aus den Kap. |4.2}|5.3.1 
und 5.4] Sowohl die stochastische Struktur, ausgelegt von einem Optik- 
ingenieur zum Testen der neuen Prozesskette, als auch die stochastische 


Struktur, bestimmt durch den genetischen Algorithmus, erfüllen die 
Vorgaben. Die periodische Freiform-Struktur mit der Zellengröße 1x10 
überschreitet beide Werte, die mit der Zellengröße 2x6 den HV-Wert. 


Dieses Kapitel zeigt, dass es möglich ist, Strukturen durch Optimie- 
rungsalgorithmen zu finden, die eine vorgegebene Gradienten-Kurve 
reproduzieren und deren Amplituden nicht auf Normalverteilungen 
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Abbildung 5.16: Simulierte LSVen mit jeweils der periodischen Freiform-Struktur mit der 
maximalen Fitness. Die gestrichelte Linie zeigt die Position, an der die Gradienten-Kurve 
bestimmt wurde. 


(vgl. Kap|5.3.1) basieren. Die Strukturen erzeugen jedoch wellenförmige 
Modulationen der LSV oder streuen soviel Licht über die HDG, dass 
gesetzliche Werte überschritten werden. Daher kann die Gradienten- 


Kurve nicht unabhängig als einzelnes Kriterium verwendet werden, 
um eine Struktur als gut oder schlecht zu bewerten. Um eine Ausle- 
gung zu ermöglichen, muss die Bewertungsfunktion angepasst 
werden, sodass diese die gesetzlichen Grenzen beachtet. Das Verhin- 
dern der wellenförmigen Modulationen ist schwieriger. Um diese zu 
verhindern, muss entweder die Bewertungsfunktion angepasst oder 
die Parameter der Strukturen so gewählt werden, dass diese nicht ent- 
stehen können. Für die Anpassung der Bewertungsfunktion muss die 
Modulation mathematisch beschrieben werden, was keine triviale Pro- 
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Tabelle 5.1: Vergleich von gesetzlichen Vorgaben für die LSV für ein Abblendlicht mit den 
Simulationsergebnissen. Die Angabe o7_1.25,1] bezieht sich auf Abb.|5.6 


o S gesetzliches 
Ren a EU) | Mas led) 
HV 0 0 625 
BR 2,5 1,0 1750 

stochastisch | 0_1251 

Bez. | nen Ger | i E 1x10 (cd) | 2x6 (cd) 
HV 598 578 724 1374 
BR 308 310 2096 1394 


blemstellung ist. Vermutlich ist das Verhindern der Modulation durch 
die Verwendung von zweidimensionalen Zellen einfacher. Wie der Ver- 
gleich zwischen den LSVen der Zellengröße 1x10 und 2x6 zeigt, können 
die Modulationen durch die Wahl der Zellengröße verringert werden. 
Bei eindimensionalen Zellen wird das Licht nur entlang einer Achse 
gestreut, wodurch die Modulationen entstehen. Durch die Verwendung 
von zweidimensionalen Zellen wird das Licht zweidimensional gestreut 
womit sich die Modulationen verhindern lassen. 
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Abbildung 5.17: Oben & Mitte: Simulation des kompletten BiLED-Moduls mit den nicht- 
stochastischen Mikrostrukturen. Die gestrichelten Linien in den LSVen geben die Position 
des Auswertebereiches für die Gradienten-Kurve an. Unten: Darstellung der Strukturen. 
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ZUSAMMENFASSUNG UND AUSBLICK 


In dieser Arbeit wurde die Prozesskette zur Auslegung und Herstellung 
von mikrostrukturierten Linsen verbessert. Die neue Prozesskette ist 
erneut in Abb.|6.1 dargestellt. Durch Entwicklung eines Strukturierungs- 


Algorithmus konnten die kommerziellen Programme Rhino und Mathe- 
matica entfernt werden (vgl. Kap. 1). Die Anzahl der zur Beschreibung 


I 


1 
' Manuelle | Automatische 1. ` Nein 


a Korrekturschleife 


SEN Freigabe ei 
Herstellung 


Abbildung 6.1: Neue Prozesskette zur Herstellung mikrostrukturierter Linsen. 


der Fläche verwendeten Datenformate konnte von fünf auf zwei redu- 
ziert werden. Dadurch reduziert sich neben der Prozesszeit auch die 
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Gefahr an Datenverluste durch Konvertierungsfehler. Einerseits koste- 
ten die vielen manuellen Schritte Zeit, andererseits benötigte Rhino ca. 
25 min um eine Fläche zu berechnen. Die Zeit für die Strukturierung 
einer Fläche konnte auf wenige Sekunden reduziert werden. Der Al- 
gorithmus zur Strukturierung wurde in das CAD-Programm CATIA 
integriert. Für die Optikingenieure unterscheidet sich die Funktion nicht 
von anderen CATIA-Funktionen. Die alte Prozesskette erforderte viel Er- 
fahrung, um alle notwendigen Programme bedienen zu können. Für die 
neue Prozesskette reichen CATIA-Kenntnisse aus. Die Funktionsweise 
wurde an zwei Beispielen erfolgreich demonstriert. 


Um eine automatisierte Auslegung der Strukturen zu ermöglichen, wur- 
de die Simulationszeit stark reduziert. Die Reduzierung basiert auf der 
Entwicklung eines auf moderne Prozessoren angepassten Algorithmus 
zu Berechnung der B-Spline Basis-Funktionen sowie der Vermeidung 
unnötiger Strahlverfolgungen bei der Simulation. Durch Anwendung 
von genetischen Algorithmen wurde die Auslegung der Mikrostruktu- 
ren automatisiert. An zwei Beispielen wurde gezeigt, dass vorgegebene 
Gradienten-Kurven reproduziert werden können. Bisher basierten die 
Amplituden der Strukturen auf Normalverteilungen. Mit dem Optimie- 
rungsalgorithmus konnten Strukturen gefunden werden, indem direkt 
die Kontrollpunkte der Flächen verschoben wurden. Die Strukturen 
erzeugen die vorgegebenen Gradienten-Kurven, überschreiten jedoch 
die gesetzlichen Blendwerte. Zukünftig könnten durch die Freiform- 
Strukturen eventuell die HDG bzw. der Verlauf der Gradienten-Kurve 
mit mehr Freiheitsgraden gestaltet werden. 


Die Prozesskette der Herstellung wurde um die Vermessung, Rückfüh- 
rung und lichttechnische Simulation der Werkzeugeinsätze erweitert. 
Dadurch wurde eine neue Analysemöglichkeit geschaffen, um eventuell 
auftretende Fehler frühzeitig zu erkennen. Damit kann die Anzahl an 
Werkzeugkorrekturen reduziert werden, wodurch Zeit und Ressourcen 
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eingespart werden. An vier Werkzeugeinsätzen wurde die Funktions- 
weise erfolgreich demonstriert. 


Ausblick 


Bei der Flächenrückführung der vermessenen Werkzeugeinsätze wurden 
die Positionen der Kontrollpunkte als Freiheitsgrade verwendet, wobei 
die Anzahl bei der Rückführung festgehalten wurde. Sollten zukünftig 
Abweichungen auftreten, die eine höhere Anzahl an Freiheitsgraden be- 
nötigt, so kann der Algorithmus angepasst werden. Zunächst würde die 
bestehende Anzahl verwendet werden. Finden sich Bereiche, in denen 
die Anpassung nicht genau genug ist, so kann durch das Hinzufügen 
weiterer Zeilen bzw. Spalten in der Kontrollpunktmatrix die Flexibilität 
der Fläche lokal erhöht werden. 


Um eine noch bessere Kontrolle über den Herstellungsprozess zu bekom- 
men, wäre eine Vermessung der strukturierten Linsenoberfläche wün- 
schenswert. Zum Zeitpunkt dieser Arbeit konnte seitens der Messtechnik 
keine Linsenoberfläche komplett vermessen werden. Durch langwellige 
Abweichungen in der Lichtaustrittsseite der Linse sowie eine Verkip- 
pung zwischen Lichteintritts- zu Lichtaustrittsseite konnten die Linsen 
auf der Messmaschine nicht ausgerichtet werden. Die Vermessung, Rück- 
führung und Simulation würde die Kontrolle des Spritzgussprozesses 
ermöglichen. Wird die Vermessung des Werkzeugeinsatzes verglichen 
mit der Vermessung der Linsenoberfläche, so lässt sich beurteilen, ob 
die Strukturen ordnungsgemäß aus dem Einsatz abgeformt wurden. 


Die Freiform-Strukturen erzeugen zu viel Blendung. Um das zu verhin- 
dern, muss eine neue Bewertungsfunktion gefunden werden, die sowohl 
die Gradienten-Kurve, als auch die gesetzlichen Blendwerte beachtet. 
Durch eine größere Versuchsreihe kann herausgefunden werden, wie 
viele Parameter minimal notwendig sind, um gewünschte Gradienten- 
Kurven zu generieren. Das Herabsenken der Anzahl an Parametern 
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reduziert die Optimierungszeit. Moderne Grafikkarten sind für schnel- 
le Ray-Tracing-Berechnungen ausgelegt 99]. Eventuell kann die 
Optimierungszeit durch die Berechnung auf Grafikkarten anstatt des 
Prozessors reduziert werden. 


Bei Strukturierung der Linsen durch den Optimierungsalgorithmus wird 
die ganze Linse mit der Struktur versehen. Die Optikingenieure struktu- 
rieren z. B. in Abhängigkeit des Bereiches, durch den das Fernlicht die 
Linse passiert. Denkbar ist, vor oder während des Optimierungsprozes- 
ses zu analysieren, durch welchen Bereich der Linse das meiste Licht auf 
die HDG trifft. Damit kann die Struktur nur in dem Bereich aufgetragen 
werden, in dem sie notwendig ist. Weitergehend könnte die Struktur in 
Form von Markenzeichen (z. B. Mercedes-Stern) aufgebracht werden, 
wodurch sich neue Design-Möglichkeiten ergeben. 


112 


LITERATUR 


[1] 


[4] 


[5] 


[8] 


[9] 


[10] 


E. F. Schubert. Light-Emitting Diodes. 2. Aufl. Cambridge 
University Press, 2006. 


R. Schedel. “LED-Scheinwerfer gewinnen Marktanteile”. ATZ - 
Automobiltechnische Zeitschrift 109.11 (2007), S. 1010-1013. 


M. Gotz, M. Kleinkes, S. Pietzonka und N. Schiermeister. 
“Künftige Entwicklungen von LED-Scheinwerfern”. ATZ - 
Automobiltechnische Zeitschrift 111.9 (2009), S. 620-627. 


M. Maier, J. Moisel und F. Herold. “Multibeam-Scheinwerfer in 
der Mercedes-Benz CLS-Klasse”. ATZ - Automobiltechnische 
Zeitschrift 117.2 (2015), S. 16-21. 


R. Lachmayer, M. Gotz, M. Kleinkes und W. Pohlmann. 
“LED-Technik im Scheinwerfer”. ATZ - Automobiltechnische 
Zeitschrift 108.11 (2006), S. 956-961. 


B. Wordenweber, J. Wallaschek, P. Boyce und D. Hoffman. 
Automotive Lighting and Human Vision. Springer-Verlab Berlin 
Heidelberg, 2007. 

E. B. Goldstein. Wahrnehmungspsychologie: Eine Einfiihrung. 
Spektrum Akademischer Verlag, 1997. 

Economic Comission for Europe. Regulation No. 123 - Uniform 


provisions concerning the approval of adaptive front-lighting systems 
(AFS) for motor vehicles. 2013. 


Dassault Systemes. CATIA V5-R2016. URL: https: 


/ /www.3ds.com/de/produkte-und-services/catia/produkte/ 


R. Anderl und D. Trippner. STEP STandard for the Exchange of 
Product Model Data. Vieweg+Teubner Verlag, 2000. 


113 


[11] 


[12] 
[13] 


[14] 


[15] 


[16] 


[17] 


[18] 


[19] 


[20] 


[21] 


[22] 


LITERATUR 


Robert McNeel & Associates. Rhinoceros 3D V5. URL: 


https:/ /www.rhino3d.com/ | 


B. Fortner. The Data Handbook. Springer-Verlag New York, 1995. 


Wolfram Research. Mathematica. URL: 
https:/ /www.wolfram.com/mathematica / 


U. S. Product Data Association. Initial Graphics Exchange 
Spcification IGES 5.3. 1996. 


J. Bliedtner und G. Gräfe. Optiktechnologie: Grundlagen - Verfahren - 
Anwendungen - Beispiele. Fachbuchverl. Leipzig im 
Carl-Hanser-Verlag, 2008. 


Bundesministerium für Bildung und Forschung. Industrie 4.0 - 
Innovationen für die Produktion von morgen. 2017. 


C. Wenzel, A. Beutler, D. Zimmermann, A. Stockfisch und 
S. Lowis. Abschlussbericht zum Verbundvorhaben: Effiziente 
Photonikproduktion durch intelligente Technologie (ePiTec). 2019. 


S. Holtz und J. Schmidt. “Schweinwerferlinse für einen 


Kraftfahrzeugscheinwerfer”. Pat. DE 10 2005 009 556 A1. 2005. 


J. Fischer. “Verfahren zum Herstellen einer Scheinwerferlinse für 
einen Kraftfahrzeugscheinwerfer”. Pat. WO 2009 012 735 A1. 
2009. 


F. Hutter, L. Kotthoff und J. Vanschoren. Automated Machine 
Learning. Springer International Publishing, 2019. 


J. H. Hamkens. “Verfahren zum Herstellen eines optischen 
Linsenelementes, insbesondere einer Scheinwerferlinse für einen 
Kraftfahrzeugscheinwerfer”. Pat. W0 2009 018 798 A2. 2009. 


R. Bonitz, R. Biertümpfel, W. Semar und M. Redey. “Optische 
Linse mit Weichzeicheneffekt”. Pat. EP 1 514 148 B1. 2008. 


114 


LITERATUR 


[23] M. Kiesel. “Zur Erzeugung einer definierten 
Overhead-Beleuchtung eingerichteter Fahrzeugscheinwerfer”. 
Pat. DE 10 2009 020 593 B4. 2017. 


[24] M. Kiesel, B. v. Blanckenhagen und E.-O. Rosenhahn. 
“Beleuchtungseinrichtung in Form eines Projektionsscheinwerfers 
für Kraftfahrzeuge”. Pat. DE 10 2008 023 551 B4. 2019. 


[25] G. G. Wang. “Definition and Review of Virtual Prototyping”. 
Journal of Computing and Information Science in Engineering 2.3 
(2002), S. 232-236. 


[26] J. C. Schaaf Jr. und F. L. Thompson. “System Concept 
Development with Virtual Prototyping”. Proceedings of the 29th 
Conference on Winter Simulation. 1997, S. 941-947. 


[27] J.D. Foley, A. van Dam, S. K. Feiner und J. F. Hughes. Computer 
Graphics: Principles and Practice. 2. Aufl. Addison-Wesley, 1996. 

[28] AG Glassner, Hrsg. An Introduction to Ray Tracing. Academic 
Press Ltd., 1989. 

[29] F. Pedrotti, L. Pedrotti, W. Bausch und H. Schmidt. Optik für 
Ingenieure: Grundlagen. Springer Berlin Heidelberg, 2007. 

[30] O. Reeb. Grundlagen der Photometrie. G. Braun Karlsruhe, 1962. 

[31] N. Metropolis und S. Ulam. “The Monte Carlo Method”. Journal of 

the American Statistical Association 44.247 (1949), S. 335-341. 

[32] N. Metropolis. “THE BEGINNING of the MONTE CARLO 

METHOD”. Los Alamos Science Special Issue 15 (1987), S. 125-130. 

[33] M. Pharr und G. Humphreys. Physically Based Rendering. 2. Aufl. 

Morgan Kaufmann Publishers Inc., 2010. 


[34] W.H. Press, S. A. Teukolsky, W. T. Vetterling und B. P. Flannery. 
Numerical Recipes in C: The Art of Scientific Computing. 2. Aufl. 
Cambridge University Press, 1992. 


115 


[35] 


[36] 


[37] 


[38] 


[39] 


[40] 


[41] 


[42] 


[43] 


[44] 


[45] 


LITERATUR 


P. E. Bezier. “Example of an Existing System in the Motor 
Industry: The Unisurf System”. Proceedings of the Royal Society of 
London. Series A 321 (1971), S. 207-218. 


C. Rabut. “On Pierre Bézier’s life and motivations”. 
Computer-Aided Design 34 (2002), S. 493-510. 

D. F. Rogers. An Introduction to NURBS: With Historical Perspective. 
Morgan Kaufmann Publishers Inc., 2001. 


I. J. Schoenberg. “Contributions to the Problem of Approximation 
of Equidistant Data by Analytic Functions”. Quarterly of Applied 
Mathematics 4.1 (1946), S. 45-99. 


R. F. Riesenfeld. “Applications of B-spline Approximation to 
Geometric Problems of Computer-aided Design.” Diss. Syracuse 
University, 1973. 

L. Piegl und W. Tiller. The NURBS Book. 2. Aufl. Springer-Verlag, 
1996. 


C. De Boor. A Practical Guide to Splines; rev. ed. Applied 
Mathematical Sciences. Springer, 2001. 


GNU Scientific Library (GSL). Version 2.5. 19. Juli 2019. URL: 


http://www.gnu.org/software/gsl/ | 


M. Galassi, J. Davies, J. Theiler, B. Gough und G. Jungman. GNU 
Scientific Library - Reference Manual, Third Edition, for GSL Version 
1.12 (3. ed.). 2009. 


geomdl. Version 5.2.8. 19. Juli 2019. URL: 


https:/ /pypi.org/project/geomdl/ 


O. R. Bingol und A. Krishnamurthy. “NURBS-Python: An 
open-source object-oriented NURBS modeling framework in 
Python”. SoftwareX 9 (2019), S. 85-94. 


116 


[46] 


[47] 


[48] 


[52] 


[53] 


[54] 


[55] 


[56] 


[57] 


LITERATUR 


Object Oriented Finite Element Solver (OOFEM). Version 2.5. 19. Juli 


2019. URL: http://www.oofem.org, 

B. Patzäk. “OOFEM - an object-oriented simulation tool for 
advanced modeling of materials and structures”. Acta Polytechnica 
52.6 (2012), S. 59-66. 


Overture - An Object-Oriented Toolkit for Solving Partial Differential 
Equations in Complex Geometry. Version 25. 19. Juli 2019. URL: 


http:/ /www.overtureframework.org/ 


R. M. Hord. The Illiac IV - The First Supercomputer. Springer-Verlag 
Berlin Heidelberg, 1982. 


C. Lomont. Introduction to Intel Advanced Vector Extensions. 2011. 


O. Abert, M. Geimer und S. A. Muller. “Direct and Fast Ray 
Tracing of NURBS Surfaces”. IEEE Symposium on Interactive Ray 
Tracing (2006), S. 161-168. 


Intel Intrinsics Guide. 19. Juli 2019. URL: https: 
/ /software.intel.com/sites/landingpage/IntrinsicsGuide/ 


Microprocessor Standards Committee. 754-2019 - IEEE Standard 
for Floating-Point Arithmetic. IEEE, 2019. 


P. Gepner. “Using AVX2 Instruction Set to Increase Performance 
of High Performance Computing Code”. Computing and 
Informatics 36 (2017), S. 1001-1018. 


S. Dreiseitl. Mathematik fiir Software Engineering. Springer Vieweg, 
2018. 


Intel® Xeon® Prozessor E5-1650 v3. 28. Marz 2020. URL: https: 
/ /ark.intel.com/content/www/de/de/ark/products/82765/ 
intel-xeon-processor-e5- 1650-v3-15m-cache-3-50-ghz.html 


G. Hager und G. Wellein. Introduction to High Performance 
Computing for Scientists and Engineers. CRC Press, Inc., 2010. 


117 


[58] 


[59] 


[60] 


[61] 


[62] 


[63] 


[64] 


[65] 


[66] 


[67] 


LITERATUR 


J. E. Smith. “A Study of Branch Prediction Strategies”. Proceedings 
of the 8th Annual Symposium on Computer Architecture. 1981, 
S. 135-148. 


Intel. Intel 64 and IA-32 Architectures Optimization Reference Manual. 
2016. 


R. Kapoor. Avoiding the Cost of Branch Misprediction. 2009. URL: 


https: / /software.intel.com/en-us /articles /avoiding-the-cost- 


of-branch-misprediction 


D. Zimmermann u.a. “A procedure for designing and 
manufacturing microstructured lenses used in automotive 
headlamps”. Advanced Optical Technologies 8.6 (2019), 483-489. 


National Highway Traffic Safety Administration. Federal Motor 
Vehicle Safety Standard 108. 2007. 


S. D. Conte und C. W. D. Boor. Elementary Numerical Analysis: An 
Algorithmic Approach. 3. Aufl. McGraw-Hill Higher Education, 
1980. 


C. Schmidt und F.-J. Kalze. “Dynamic Cut-Off-Line geometry as 
the next step in forward lighting beyond AFS”. International 
Symposium on Automotive Lighting. 2007. 


C. Schmidt, F.-J. Kalze und T. Irmscher. “Illumination Strategies 
for Dynamic Headlamp Functions like Adaptive and Vertical 
Cut-Off-Line”. International Symposium on Automotive Lighting. 
2009. 


J. Moisel, R. Ackermann und M. Griesinger. “Adaptive 
Headlights utilizing LED-Arrays”. International Symposium on 
Automotive Lighting. 2009. 


C. Brecher, C. Baum, B. Meiers, D. De Simone und R. Krappig. 
Kunststoffkomponenten fiir LED-Beleuchtungsanwendungen. Springer 
Vieweg, 2016. 


118 


[68] 


[69] 


[70] 


[71] 


[72] 


[73] 


[74] 


[75] 


[76] 


LITERATUR 


R. L. Rhorer und C. J. Evans. “Fabrication of Optics by Diamond 
Turning”. In: Handbook of Optics. Hrsg. von M. Bass. New York: 
McGraw-Hill, Inc., 1995. Kap. 41, S. 1-13. 


T. Moriwaki und E. Shamoto. “Ultraprecision Diamond Turning 
of Stainless Steel by Applying Ultrasonic Vibration”. CIRP Annals 
40.1 (1991), S. 559 -562. 


C. Brecher, F. Klocke, M. Winterschladen und M. Heselhaus. 
“Ultraschallunterstiitztes Hartdrehen für die Fertigung von 
gehärteten Präzisionsstahlbauteilen”. wt Werkstatttechnik online 4 
(2006), S. 396-401. 


F. Schneider, J. Das, B. Kirsch, B. Linke und J. C. Aurich. 
“Sustainability in Ultra Precision and Micro Machining: A 
Review”. International Journal of Precision Engineering and 
Manufacturing-Green Technology 6.3 (2019), S. 601-610. 


E. Brinksmeier und R. Gläbe. “Elliptical Vibration Cutting Steel 
with Diamond Tools.” Proceedings of 14th Annual ASPE Meeting 
(1999). 

Innolite GmbH. IL 300. URL: https: / /innolite.de/wp- 
content/uploads/2018/pdf/Innolite_IL300_Brochure_en.pdf 
C. Wenzel, A. Stockfisch, T. Weske und D. Zontar. “Effiziente 


Produktion von komplexen Kunststoffoptiken”. wt 
Werkstatttechnik online 4 (2019), S. 388-393. 


M. Brozio, D. Sensen, C. Wenzel, S. Freutel und C. Brecher. 
“Surface adaptive fast axis ultra precision turning”. Proceedings of 
SPIE. 2018, 10544:1-6. 


The HDF Group. Hierarchical Data Format, version 5. 2019. URL: 
http:/ /www.hdfgroup.org/HDF5/ 


119 


[77] 


[78] 


[79] 


[80] 


[81] 


[82] 


[83] 


LITERATUR 


A. Beutler. “Flexible, non-contact and high-precision 
measurements of optical components”. Surface Topography: 
Metrology and Properties 4.2 (2016), 024011:1-12. 


A. Weckenmann und W. Hartmann. “Function-oriented method 
for the definition and verification of microstructured surfaces”. 
Precision Engineering 37.3 (2013), S. 684 -693. 


J. Stock, M. Beier, J. Hartung, S. Merx und H. Gross. “Simulation 
and analysis of optical imaging systems including real freeform 
components”. Advanced Optical Technologies 8.2 (2019), S. 111-117. 


F. Bornemann. Numerische lineare Algebra. Springer Spektrum, 
2018. 


M. Hemmerich. “Evaluation und Implementierung von 
Algorithmen zur Rückführug von Punktwolken in 
NURBS-Flächen zur Verwendung in lichttechnischen 
Simulationsprogrammen”. Masterarbeit. Hochschule Bochum, 
2018. 


G. Becker, M. Schäfer und A. Jameson. “An advanced NURBS 
fitting procedure for post-processing of grid-based shape 
optimizations”. 49th AIAA Aerospace Sciences Meeting Including the 
New Horizons Forum and Aerospace Exposition (2011). 


Y. Zhang, H.-N. Cheng, R. Wu und R. Liang. “Data processing for 
point-based in situ metrology of freeform optical surface”. Optics. 
Express 25.12 (2017), S. 13414-13424. 


N. Carlsson. “NURBS Surface Fitting with Gauss-Newton”. 
Masterarbeit. Lulea University of Technology, 2009. 


M. Randrianarivony und G. Brunnett. “Parallel Implementation 
of Surface Reconstruction From Noisy Samples”. TU Chemnitz - 
Sonderforschungsbereich 393 (2002). 


120 


[86] 


[87] 


[88] 


[89] 


[91] 


[92] 


[93] 


LITERATUR 


M. Randrianarivony und G. Brunnen. “Approximation by 
NURBS curves with free knots”. Proceedings of the Vision, Modeling, 
and Visualization Conference. 2002, S. 195-201. 


D. Brujic, I. Ainsworth und M. Ristic. “Fast and accurate NURBS 
fitting for reverse engineering”. The International Journal of 
Advanced Manufacturing Technology 54 (2011), S. 691-700. 


Y. Kineri, M. Wang, H. Lin und T. Maekawa. “B-spline surface 
fitting by iterative geometric interpolation /approximation 
algorithms”. Computer-Aided Design 44.7 (2012), S. 697 -708. 

S. Hu und J. Wallner. “A second order algorithm for orthogonal 
projection onto curves and surfaces”. Computer Aided Geometric 
Design 22 (2005), S. 251-260. 

X. Li u.a. “Convergence Analysis on a Second Order Algorithm 
for Orthogonal Projection onto Curves”. Symmetry 9 (2017), 
210:1-13. 

I. N. Bronstein, K. A. Semendjajew, G. Musiol und H. Mühlig. 
Taschenbuch der Mathematik. Verlag Harri Deutsch, 2001. 

C. Schwanengel. Comparison of techniques for measuring luminous 
intensity distribution overall and across segments. TechnoTeam 
Bildverarbeitung GmbH, 2010. 

D. Zimmermann und C. Neumann. “Simulation und Design 
mikrostrukturierter Linsen für automobile Scheinwerfer”. Lux 
Junior. 2019. 

K. Höschel und V. Lakshminarayanan. “Genetic algorithms for 
lens design: a review”. Journal of Optics 48.1 (2018), 134-144. 

R. L. Haupt und S. E. Haupt. Practical Genetic Algorithms. 
Wiley-Interscience, 2004. 

S. Luke. Essentials of Metaheuristics. 2. Aufl. Available for free at 
http://cs.gmu.edu/~sean/book/metaheuristics/. Lulu, 2013. 


121 


[97] 


[98] 


[99] 


LITERATUR 


H. Mühlenbein und D. Schlierkamp-Voosen. “Predictive Models 
for the Breeder Genetic Algorithm I. Continuous Parameter 
Optimization”. Evolutionary Computation 1.1 (1993), S. 25-49. 


V. Sanzharov, A. Gorbonosov, V. Frolov und A. Voloboy. 
“Examination of the Nvidia RTX”. Proceedings of the 29th 
International Conference on Computer Graphics and Vision. 2019, 
S. 7-12. 


NVIDIA. NVIDIA Turing GPU Architecture (WP-09183-001_v01). 
2018. 


122 


EIGENE VERÖFFENTLICHUNGEN 


[61] D. Zimmermann u.a. “A procedure for designing and 
manufacturing microstructured lenses used in automotive 
headlamps”. Advanced Optical Technologies 8.6 (2019), 483-489. 


[93] D. Zimmermann und C. Neumann. “Simulation und Design 
mikrostrukturierter Linsen für automobile Scheinwerfer”. Lux 
Junior. 2019. 


BETREUTE STUDENTISCHE ARBEITEN 


[851] M. Hemmerich. “Evaluation und Implementierung von 
Algorithmen zur Rückführug von Punktwolken in 
NURBS-Flachen zur Verwendung in lichttechnischen 
Simulationsprogrammen”. Masterarbeit. Hochschule Bochum, 
2018. 


123 


TABELLENVERZEICHNIS 


2.1 


Vergleiche der Leistung zwischen konventioneller und 


SIMD-basierter Auswertung von B-Spline-Flächen. Die 


Angaben zu den Implementierungen (Imp.) beziehen sich 


auf den Anhang ..... 2... 2.22... 


32 


5.1 


Vergleich von gesetzlichen Vorgaben fiir die LSV fiir ein 


Abblendlicht mit den Simulationsergebnissen. Die Angabe 


01_1.25,1] bezieht sich auf Abb.5.6|| ........ 2.2... 


125 


ABBILDUNGSVERZEICHNIS 


1.1 


BiLED-Scheinwerfermodul| . .... 2. 2.222 22222... 


1.2 Vergleich zwischen einer scharfen und einer weichen 


HDG!....... 


1.3 Bisher verwendete Prozesskette zur Herstellung 


mikrostrukturierte 


rLinsen.|........0..... 000084 


1.4 


Die zwei vertikalen und horizontalen Streifen stellen den 


lokalen Verlust der Mikrostruktur durch die Übertragung 


von Rhino zuCapadar,,,,,, a 


2.1 


Schematische Darstellung des Ray-Tracing.)......... 


10 


2.2 


Schematische Darstellung einer Hierarchie von 


Begrenzungsrahmen (Begrenzungsvolumen) einer Kurve 


(blau). Die vom Strahl getroffenen Begrenzungsrahmen sind 


grün markiert.|. . 


2.3 


Ausrichtung des Koordinatensystems.|............ 


13 
14 


2.4 


Vergleich von zwei B-Spline-Kurven mit den gleichen 


Kontrollpunkten u 


nd unterschiedlichen Knotenvektoren. 


Die äußeren Knote 


n sind vierfach in dem Knotenvektor 


vorhanden. ... 


16 


2.5 


Darstellung einer B-Spline-Fläche (rechts) und die 


dazugehörigen Kontrollpunkte und das Kontrollpolygon 


(links). Die Knotenvektoren in u- und v-Richtung sind [0, 0, 


0,0,%,%1,1,1,1]. 


17 


2.6 


Grafische Darstellung der Abhängigkeit gegeben durch Gl. 


(2.8) und Gl. (2.9) für kubische (p = 3) 


Basis-Funktionen. 


127 


18 


ABBILDUNGSVERZEICHNIS 


3.1 


Darstellung des Zusammenhangs zwischen der Lichtstärke 


I, und der Gradienten-Kurve G(ß) für eine scharfe (blau) 


und weiche (orange) HD", 


36 


3.2 


Schematische Darstellung des ersten Schrittes der 


Strukturierung: Generierung eines äquidistanten 


Punktgitters.).: »= +2: 6s eee eee Sea bea gabe 


38 


3.3 


Schematische Darstellung des zweiten Schrittes der 


Strukturierung: Projektion des Punktgitters.| agate eh ate. 


39 


3.4 


Schematische Darstellung des dritten Schrittes des ersten 


Verfahrens der Strukturierung: Verschiebung des 


Punktgitters entlang der Normalenvektoren.)........ 


39 


3.5 


Schematische Darstellung des vierten Schrittes des ersten 


Verfahrens der Strukturierung: Berechnung einer 


interpolierenden Fläche.| EE 


3.6 


Schematische Darstellung des fünften Schrittes der 


40 


Strukturierung: Zuschneiden der Fläche.| .......... 


41 


3.7 


Schematische Darstellung des zweiten Verfahrens. Links: 


Interpolierende Fläche des Punktgitters. Rechts: 


Strukturierung der Fläche durch Verschieben der 


Kontrollpunkte. | R E E E a A ee 


42 


3.8 


Einteilung der Gitterpunkte (vgl. Abb 3. 2) in Zellen der 


Größe 2x2 zur Zuordnung von Normalverteilungen. 


47 


3.9 


Streuverhalten strukturierter Platten für Zellen der Größe 


ER HEEN 


48 


3.10 _Strukturhöhen der Platten für Zellen der Größe 1x1. Die 


grauen Linien kennzeichnen die Zellengréfen.|....... 


49 


3.11 Streuverhalten strukturierter Platten für Zellen der Größe 


BEER 


128 


50 


ABBILDUNGSVERZEICHNIS 


3.12 Strukturhöhen der Platten für Zellen der Größe 2x2. Die 


grauen Linien kennzeichnen die Zellengrößen. Die 


Farbskala ist unterhalb von —1,5 um und oberhalb von 


1,5 pm logarithmisch.|.... 2.222222 


51 


3.13 Streuverhalten strukturierter Platten für Punktgitter mit 


unterschiedlichen Auflösungen in x- und y-Richtung.|. . . 


52 


3.14 Strukturhöhen der Platten für Punktgitter mit 


unterschiedlichen Auflösungen in x- und y-Richtung. Die 


grauen Linien kennzeichnen die Zellengrößen. Zur besseren 


Sichtbarkeit sind nur die Zellengrenzen in der Richtung der 


größeren Ausdehnung eingezeichnet.| ............ 


53 


3.15 Links: Werkzeughälfte mit Einsatz für die Lichtaustrittsseite 


der Linse sowie spritzgegossene Linsen. Rechts: 


Werkzeughalfte mit Lichteintrittseite der Linse.| ...... 


3.16 Werkzeugeinsatz mit Mikrostruktur.|. 2.2.2... 


54 
55 


3.17 Links: Im Projekt verwendete Drehmaschine „Innolite IL 


300” [73]. Rechts: Bauteil und Diamantwerkzeug auf der 


Drehmaschine Tal, 
3.18 Schematische Darstellung eines polaren Punktgitters.). . . 


56 
56 


3.19 Links: In dem Projekt verwendete Messmaschine Mahr 


MFU200 [77]. Rechts: Werkzeugeinsatz auf Messmaschine 


(Sr EEGEN 


57 


3.20 Bewegung des Messkopfes zur Justierung (A) und Messung 


= 


58 


3.21 Kontrollpunkte (rot) und Kontrollpolygon (blau) einer 


B-Spline Fläche. Die Lücke zwischen ug und un dient der 


Veranschaulichung.|. . 2... 0... a 


129 


60 


ABBILDUNGSVERZEICHNIS 


3.22 Ein bikubisches Kontrollpolygon [88]. Die schwarzen 


unterbrochenen Linien geben die Positionen der inneren 


Knoten an. Das schwarze Kreuz markiert den projizierten 


Punkt. Die blau markierten Kontrollpunkte beeinflussen die 


Fläche an der Stelle mit dem schwarzen Kreuz.|....... 


3.23 Schematische Darstellung der orthogonalen Projektion 


65 


(modifiziert von EN. . .. zz ea eed eee eee E A 


66 


3.24 Konvergenzverhalten des Algorithmus|3.1'für eine 


mikrostrukturierte Linse] ...............000. 


68 


3.25 Simulierte Gradienten-Kurven in Abhängigkeit eines 


künstlichen Messrauschens und der Auflösung der 


Ausgangsflache.|... 2... eoe n nennen 


4.1 Neue Prozesskette zur Herstellung mikrostrukturierter 


Linsen]. «Ania SB Rea we A ee a ben 


70 


72 


4.2 Schematische Darstellung einer 


Bildverarbeitungs-Anlage.|. .. 2.2... 2 2222.00. 


73 


4.3 Aufbau des BiLED-Moduls, das zum Testen der neuen 


Prozesskette verwendet wird.|...............00. 


74 


4.4 Simulation der LSV und der Gradienten-Kurve des 


BiLED-Moduls mit nicht-strukturierter Linse. Die 


gestrichelte Linie in der LSV gibt die Position der 


Gradienten-Kurvean.| .......0..0..0.0 00000004 


75 


4.5 Nominelle Mikrostruktur zur Aufweichung der HDG für 


den ECE-Bereich. Die Farbskala ist unter —0,5 um und über 


0,5umlogarithmisch. .......2 222220 


130 


76 


ABBILDUNGSVERZEICHNIS 


4.6 Simulation der LSV und der Gradienten-Kurve des 


BiLED-Moduls mit einer strukturierten Linse ausgelegt fiir 


den ECE-Bereich im Vergleich zu gemessenen 


Gradienten-Kurven von abgeformten Linsen. Die 


gestrichelte Linie in den LSVen gibt die Position der 


Gradienten-Kurvean.| ........0...0. 000000004 


77 


4.7 Abweichung zwischen Flächenrückführung und Messdaten 


für die zwei Werkzeugeinsätze mit Struktur für 


ECE-Gebiete. Links: Werkzeugeinsatz Nr. 1, Rechts: 


Werkzeugeinsatz Nr. 2... .. 2222222 


78 


4.8 Vergleich der Strukturhöhen zwischen den nominellen und 


den gemessenen Daten für Werkzeugeinsatz1.|....... 
4.9 Foto einer Linse mit der ECE-Struktur|............ 


79 
79 


4.10 Nominelle Mikrostruktur zur Aufweichung der HDG für 


den FMVSS-Bereich. Die Farbskala ist unter —0,5 um und 


über 0,5 um logarithmisch. . . 2.2.2.2 2222.20. 


80 


4.11 Simulation der LSV und der Gradienten-Kurve des 


BiLED-Moduls mit einer strukturierten Linse, ausgelegt für 


den FMVSS-Bereich im Vergleich zu gemessenen 


Gradienten-Kurven von abgeformten Linsen. Die 


gestrichelte Linie in den LSVen gibt die Position der 


Gradienten-Kurveanl,,,,, a a 


82 


4.12 Abweichung zwischen Flächenrückführung und Messdaten 


für die zwei Werkzeugeinsätze mit Struktur für 


FMVSS-Gebiete. Links: Werkzeugeinsatz Nr. 3, Rechts: 


Werkzeugeinsatz Nr. Ai. 


83 


4.13 Vergleich der Strukturhöhen zwischen den nominellen und 


den gemessenen Daten für Werkzeugeinsatz 3.|....... 


131 


83 


ABBILDUNGSVERZEICHNIS 


4.14 Foto einer Fresnel-Reflexion einer Deckenleuchte an der 


Lichteintrittsseite der Linse. Die Verzerrung der 


rechteckigen Leuchte entsteht durch den Einfall der 


nominell planaren Fläche.|..... 2... 222.220... 


84 


5.1 


Schematische Darstellung des Aufbaus zur Vorbereitung des 


Optimierungsprozesses aE & PS al EG ee EOE SEG 


86 


5.2 


Filterung der Strahlen: Es werden nur die Strahlen behalten, 


die in einen Bereich in der Nähe der HDG treffen (weißes 


Rechteck).| : + 222. : 2.0: ee ee ee ee m 


87 


5.3 


Schematische Darstellung des Ablaufes eines genetischen 


Algorithmus) m. vw werner. 


90 


5.4 


Vergleich zwischen der Soll- und den resultierenden 


Gradienten-Kurven für die drei Optimierungen mit 


unterschiedlich gefilterten Strahlen-Lichtquellen.| ..... 


93 


5.5 


Konvergenzverhalten des genetischen Algorithmus bei der 


Optimierung von stochastischen Strukturen für das 


BiLED-Modul. Links: Maximale Fitness pro Generation. 


Rechts: Durchschnittliche Fitness pro Generation.| .... . 


94 


5.6 


Simulierte LSVen für das BiLED-Modul mit jeweils der 


stochastischen Struktur mit der maximalen Fitness. Die 


gestrichelte Linie zeigt die Position, an der die Gradienten- 


Kurve bestimmt wurde le, 


95 


5.7 


Simulation des kompletten BiLED-Moduls mit einer 


Mikrostruktur erstellt mit den Parametern o: 1. Die 


gestrichelte Linie markiert die Position zur Auswertung der 


Gradienten-Kurve. Rechts: Vergleich zwischen Soll- und 


Ist-Gradienten-Kurve] . . 2.222 oo non 


5.8 


Aufbau des MLS-BiLED-Scheinwerfermoduls.|....... 


132 


96 
97 


ABBILDUNGSVERZEICHNIS 


5.9 Simulation des kompletten MLS-BiLED-Moduls mit einer 


Mikrostruktur ausgelegt durch einen Optikingenieur. 


Rechts: Vergleich zwischen der simulierten (griin) und der 


fiir die Optimierung verwendeten Soll-Gradienten-Kurve 


(OD) 2s giga oe eee aw eke han en ch ba ea 
5.10 Vergleich zwischen der Soll- und den resultierenden 


98 


Gradienten-Kurven fiir die zwei Optimierungen mit 


unterschiedlich gefilterten Strahlen-Lichtquellen.| ..... 


99 


5.11 Konvergenzverhalten des genetischen Algorithmus bei der 


Optimierung von stochastischen Strukturen für das 


MLS-BiLED-Modul. Links: Maximale Fitness pro Generation. 


Rechts: Durchschnittliche Fitness pro Generation.|.... . 


100 


5.12 Simulierte LSVen für das MLS-BiLED-Modul mit jeweils der 


stochastischen Struktur mit der maximalen Fitness. Die 


gestrichelte Linie zeigt die Position, an der die Gradienten- 


Kurve bestimmt wurde le, 


101 


5.13 Simulation des kompletten MLS-BiLED-Moduls mit einer 


Mikrostruktur erstellt mit den Parametern oz 2.59). 


Rechts: Vergleich zwischen Ziel- und Ist-Gradienten- 


Kurve) ereer ee ae wa 


101 


5.14 Vergleich zwischen der Soll- und den resultierenden 


Gradienten-Kurven fiir die zwei Optimierungen mit 


unterschiedlichen Freiheitsgraden der Mikrostrukturen.|. . 


5.15 Konvergenzverhalten des genetischen Algorithmus bei der 


103 


Optimierung von periodischen Freiform-Strukturen. Links: 


Maximale Fitness pro Generation. Rechts: Durchschnittliche 


Fitness pro Generation... 


104 


5.16 Simulierte LSVen mit jeweils der periodischen 


Freiform-Struktur mit der maximalen Fitness. Die 


gestrichelte Linie zeigt die Position, an der die Gradienten- 


Kurve bestimmt wurde le, 
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ABBILDUNGSVERZEICHNIS 


5.17 Oben & Mitte: Simulation des kompletten BiLED-Moduls 


mit den nicht-stochastischen Mikrostrukturen. Die 


gestrichelten Linien in den LSVen geben die Position des 


Auswertebereiches für die Gradienten-Kurve an. Unten: 


Darstellung der Strukturen... . 2... 22.22.20... 107 
6.1 Neue Prozesskette zur Herstellung mikrostrukturierter 
Linsen.|| si. 4 0004 hee ee Oe waren 109 
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ABKÜRZUNGEN UND SYMBOLE 


ABKÜRZUNGEN 


AK31 ArbeitsKreis 31 (alternative Bezeichnung für ,,Gradient”) 


AVX Advanced Vector Extensions 

ASCII American Standard Code for Information Interchange 
BMBF Bundesministerium fiir Bildung und Forschung 

CAD Computer-Aided Design 

ECE Economic Comission for Europe 


ePiTec effiziente Photonikproduktion durch intelligente 
Technologie 


FMVSS Federal Motor Vehicle Safety Standards 

HDF5 Hierarchical Data Format 5 

HDG Hell-Dunkel-Grenze 

IGES Initial Graphics Exchange Specification 

LED Light-Emitting Diode 

LSV LichtStarkeVerteilung 

NURBS Non-Uniform Rational Basis-Spline 

SIMD Single Instruction, Multiple Data 

STEP STandard for the Exchange of Product model data 
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ABKÜRZUNGEN UND SYMBOLE 


LATEINISCHE SYMBOLE UND VARIABLEN 


C(u) B-Spline-Kurve 

G Gradient bzw. AK31-Wert 

ly Lichtstärke 

J Jacobi-Matrix 

n Index des letzten Kontrollpunktes einer B-Spline-Fläche bzw. 
Kurve in u-Richtung 

Nip B-Spline-Basisfunktion vom Grad p 

p Grad der B-Spline-Kurve bzw. -Fläche in u-Richtung 

D Kontrollpunkt einer B-Spline-Fläche bzw. -Kurve 

q Grad der B-Spline-Flache in v-Richtung 

Qij Gitterpunkt 


S(u,v) B-Spline-Fläche 
R Richtungsvektor eines Lichtstrahls 


S-(u,v) Eine Koordinate der B-Spline-Flache mit c € x,y,z 


u; Knoten einer B-Spline-Kurve bzw. -Fläche 

U Erster Knotenvektor einer B-Spline-Kurve bzw. -Fläche 
Üi Knoten einer B-Spline-Fläche 

V Zweiter Knotenvektor einer B-Spline-Fläche 
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ABKÜRZUNGEN UND SYMBOLE 


GRIECHISCHE SYMBOLE UND VARIABLEN 


a Horizontaler Winkel 

b Vertikaler Winkel 

be Standardabweichung einer Normalverteilung 
e Qualitätsfunktion bzw. Fitness-Wert 

D Parameter einer Mikrostruktur 

a Parameter einer Mikrostruktur 


OPERATOREN UND MATHEMATISCHE SYMBOLE 


© Elementweise Multiplikation 


(A,B),  Frobenius-Skalarprodukt 
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Programmauflistung 6.1: C++ Implementierung des Algorithmus 2.1 


1 void NubsConventionalDiss::compUBasis(int span, double param, 


O oo NAURAN 


Pe 
mM © 


el o a 
no GF fF WN 


oo N ABAUN 


VectorN& basis) { 
ine at, Fe 
double temp, saved; 


basis[0] = 1.0; 


for ( i = 1; i < ndu; i++ ) { 


m_pLeft[i] = param - (*knot_u_ptr)[span+1-i ]; 


m_pRight[i] = (*knot_u_ptr)[span+i] - param; 


saved = 0.0; 
for ( j = 0; j 


<i; j++) { 


temp = basis[j] / ( m_pRight[j+1] + m_pLeft[i-j] ) 


7 


basis[j] = 


saved + m_pRight[j+1]*temp; 


saved = m_pLeft[i-j ]*temp; 


} 


basis[i] = saved; 


Programmauflistung 6.2: C++ Implementierung des Algorithmus/2.3 


__m256d CubicNubsSurface :: getBasisFunsU (int span, double param 


) const { 


const int spanMinusDegree = span-m_degree; 


const __m256d uVec = 


_mm256_set1_pd (param) ; 


// Contains [N_{i-3, 0}, N_{i-2, 0}, N_{i-1, 0}, N_{i, 
__m256d Nim3 = _mm256_setr_pd(0., 0., 0., 1.); 


// Calculate [N_{i-3, 


Uh, ROM =% 
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this ->getNextDegreeBasisFuns (spanMinusDegree , m_cllu, m_c2lu 
, m_c3lu, m_c4lu, uVec, Nim3); 


Hi Calllewillaiz® |[IN iS, Pal NEE A INL, AP INE, | 
this ->getNextDegreeBasisFuns (spanMinusDegree , mucl2u, m_c22u 
, m_c32u, m_c42u, uVec, Nim3); 


ji Cailletillenze [INL =e), Ste Nora, al, Nena, yl, N, 
this->getNextDegreeBasisFuns (spanMinusDegree, m_cl3u, m_c23u 


, m_c33u, m_c43u, uVec, Nim3); 


return Nim3; 


void CubicNubsSurface :: getNextDegreeBasisFuns ( 


int spanMinusDegree, const std :: vector<double>& cl, 

const std:: vector<double>& c2, const std:: vector<double>& c3 
const std:: vector<double>& c4, const __m256d& paramVec, 
__m256d &basisFuns) const { 


in Calculate [INE — 25 pi NEN ed NER NEE On p] 
__m256d Nim2 = _mm256_permute4x64_pd(basisFuns, 57); 

// Calculate [N_{i-2, pi, N_{i-1, p}, N_{i, p}, 0] 

Nim2 = _mm256_blend_pd(Nim2, _mm256_setzero_pd(), 8); 


// load constants 

const _ m256d clVec = _mm256_loadu_pd(&(c1[spanMinusDegree ]) 

De c3Vec = _mm256_loadu_pd(&(c3[spanMinusDegree ]) 

ee c2Vec = _mm256_loadu_pd(&(c2[spanMinusDegree ]) 

Er: A c4Vec = _mm256_loadu_pd(&(c4[spanMinusDegree ]) 
); 


// cl multiply u add c3 
const __m256d clmupc3 = _mm256_fmadd_pd(paramVec, c1Vec, 


c3Vec); 
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40 // c2 multiply u add c4 
41 const __m256d c2mupc4 = _mm256_fmadd_pd(paramVec, c2Vec, 
42 c4Vec); 


44 dE egen eet Wi, ING fia, Wi, INE, ir, INE, aly 


45 basisFuns = _mm256_add_pd(_mm256_mul_pd(clmupc3, basisFuns), 
46 _mm256_mul_pd(c2mupc4, Nim2) ) ; 


Programmauflistung 6.3: C++ Implementierung des Algorithmus]/2.2 


1 void NubsConventionalDiss::compUDers( int span, double param, 
int n, Matrix& derU ){ 


2 nenne: 

3 int si S2; 

4 int jk ‚uk; 

5 double temp, saved; 

6 double d; 

7 att ER 

8 

9 (*m_pBasisTmp) [0][0] = 1.0; 

10 for ( i = 1; i < ndu; i++ ) { 

11 m_pLeft[i] = param — (*knot_u_ptr)[span+1-i ]; 

12 m_pRight[i] = (*knot_u_ptr)[span+i] - param; 

13 saved = 0.0; 

14 or ( jj = j ei jaar i 

15 (*m_pBasisTmp)[i][j] = m_pRight[j+1] + m_pLeft[i-j ]; 
16 temp = (*m_pBasisTmp)[j][i-1] / (*m_pBasisTmp) [i ][j ]; 
17 (*m_pBasisTmp)[j ][i] = saved + m_pRight[j +1]*temp; 
18 saved = m_pLeft[i-j]*temp; 

19 } 

20 (*m_pBasisTmp)[i][i] = saved; 

2 | 

22 

23 for ( i = 0; i < ndu; i++ ) { 

24 derU[0][i] = (*m_pBasisTmp) [i ][ndu-1]; 

25 } 

26 


27 for ( j = 0; j < ndu; j++ ) { 
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= 0; 

= i; 

m_pColsTmp) [0][0] = 1.0; 
for ( k= 1; k <= n; k++ ) { 

el = 0; 

jk = j - k; 


uk = ndu -1 -k; 
jak) 
(*m_pColsTmp)[s2][0] = (*m_pColsTmp)[s1][0] / (* 
m_pBasisTmp) [uk+1][jk ]; 
d = (*m_pColsTmp)[s2][0] * (*m_pBasisTmp) [jk ][uk]; 
} 
i ( je Ss sil )) wil = ip 


else il = -jk; 
if ( j-1 <= uk ) i2 =k; 
else i2 = ndu - j; 


for ( al = lee) il 
(*m_pColsTmp)[s2][i] = ( (*m_pColsTmp)[s1][i] - (* 
m_pColsTmp)[s1][i-1] ) 
/ (*m_pBasisTmp)[uk+1][jk+i]; 
d += (*m_pColsTmp)[s2][i] * (*m_pBasisTmp)[jk+i][uk]; 
} 
ni j e a I 
(*m_pColsImp) [s2][k] = -(*m_pColsTmp)[s1][k-1] / (* 
m_pBasisImp) [uk+1][j ]; 
d += (*m_pColsTmp)[s2][k] * (*m_pBasisTmp) [j ][ uk]; 
} 
derU[k][j] = d; 


ı = oil? 
Sil = el: 
G2 = i? 
ndu - 1; 


for ( k = 1; k <= n; k++ ) { 
for ( i = 0; i < ndu; i++ ) { 


} 


derU[k][i] *= j; 
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Programmauflistung 6.4: C++ Implementierung des Algorithmus/2.4 


__m256d CubicNubsSurface :: getBasisDerivat(int span, 


double param m256d &firstDerivatives) const 


ee 


const int spanMinusDegree = span-m_degree; 
const _m256d uVec = _mm256_setl_pd(param) ; 


// Contains [N_{i-3, 0}, N_{i-2, 0}, N_{i-1, 0}, N_{i, Ou 
__m256d Nim3 = _mm256_setr_pd(0., 0., 0., 1.); 


fi Calculate Jett ët Wi, NE, I, INL, Wh, Ne, iyi 


this ->getNextDegreeBasisFuns (spanMinusDegree , m_cllu, m_c21lu 
, m_c3lu, m_c4lu, uVec, Nim3) ; 


ff <alenlase NETS 2), Nehl-2, 2, ING, 22 INL, 2 


this->getNextDegreeBasisFuns (spanMinusDegree, m_cl2u, m_c22u 
, m_c32u, m_c42u, uVec, Nim3); 


Vi Calculate BEE S NE EK a EE EARCH N a S 
ande (INA 3a ZU a ECO tL, 
a, NA, Bil 
this->getNextDegreeBasisFunsDerivative (spanMinusDegree, 
m_cl3u, m_c23u, m_c33u, m_c43u, uVec, Nim3, 


firstDerivatives); 


return Nim3; 


void CubicNubsSurface :: getNextDegreeBasisFunsDerivative ( 
int spanMinusDegree, const std:: vector<double>& cl, 
const std::vector<double>& c2, const std:: vector<double>& 
CS; 
const std :: vector<double>& c4, const __m256d& paramVec, 
__m256d& basisFuns, __m256d& basisDerivative) const 
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// Calculate [N_{i-2, p}, N_{i-1, p}, N{i, p}, N_{i-3, p}] 
__m256d Nim2 = _mm256_permute4x64_pd(basisFuns, 57); 

// Calculate [N_{i-2, pl, N_{i-1, p}, N{i, p}, 0] 

Nim2 = _mm256_blend_pd(Nim2, _mm256_setzero_pd(), 8); 


// load constants 


E 
const __m256d c3Vec = _mm256_loadu_pd(&(c3[spanMinusDegree ]) 
J 
const __m256d c2Vec = _mm256_loadu_pd(&(c2[spanMinusDegree ]) 
Ds 
const _ m256d c4Vec = _mm256_loadu_pd(&(c4[spanMinusDegree ] ) 
D: 


const _ m256d c1Vec = _mm256_loadu_pd(&(c1[spanMinusDegree ]) 


// el multiply u add c3 

const __m256d clmupc3 = _mm256_fmadd_pd(paramVec, c1Vec, 
c3Vec); 

// c2 multiply u add c4 

const __m256d c2mupc4 = _mm256_fmadd_pd(paramVec, c2Vec, 
c4Vec); 


// We calculate the derivatives before the next degree basis 
functions, because we need the previous basis functions 
whose registers are overwritten later 

basisDerivative = _mm256_add_pd (_mm256_mul_pd(c1Vec, 

basisFuns), _mm256_mul_pd(c2Vec, Nim2)); 

basisDerivative = _mm256_mul_pd(basisDerivative , 


mm256_setl_pd(m_degree)); 


// Calculate [IN {i-3, p+1}, N_{i—2, p+1}, N(i-1, p+1}, N {i 
o jor ily] 
basisFuns = _mm256_add_pd(_mm256_mul_pd(clmupc3, basisFuns) , 
_mm256_mul_pd(c2mupc4, Nim2) ) ; 
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Programmauflistung 6.5: Bisher verwendete C++ Implementierung der Gl. zur 
Auswertung einer B-Spline-Fläche. 


1 HPoint NubsConventionalDiss :: getPointFromUV ( 


2 const UVPoint& uv ) { 

3 int 1,k,indU,indV; 

4 int spanU,spanV; 

5 HPoint pt; 

6 

7 getSpanUV ( uv, spanU, spanV ); 

8 compUBasis ( spanU, uv.u(), *m_pBasisU ); 

9 compVBasis ( spanV, uv.v(), *m_pBasisV ); 

10 

11 for ( 1 = 0; 1 < ndv; l1++ ) { 

12 indV = spanV - ndv + 1 + l; 

13 

14 (*m_pMatrixTmp)[1][0] = 0.0; 

15 (*m_pMatrixTmp)[1][1] = 0.0; 

16 (*m_pMatrixTmp)[1][2] = 0.0; 

17 

18 for ( k= 0; k < ndu; k++ ) { 

19 indU = spanU - ndu + 1 + k; 

20 (*m_pMatrixTmp)[1][0] += (*m_pBasisU)[k] * (* 
c_points_x_ptr) [indU ][indV ]; 

21 (*m_pMatrixTmp)[1][1] += (*m_pBasisU)[k] * (* 
c_points_y_ptr) [indU][indV ]; 

22 (*m_pMatrixTmp)[1][2] += (*m_pBasisU)[k] * (* 
c_points_z_ptr) [indU ][indV ]; 

23 } 

24} 

25 for ( 1 = 0; 1 < ndv; I+) { 

26 pt[0] += (*m_pBasisV)[1] * (*m_pMatrixTmp) [1][0]; 

27 pt[1] += (*m_pBasisV)[1] * (*m_pMatrixTmp)[1][1]; 

28 pt[2] += (*m_pBasisV)[1] * (*m_pMatrixTmp) [1][2]; 

29} 

30 

3 return pt; 

32 } 
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Programmauflistung 6.6: C++ Implementierung des Algorithmus |2.5|zur Auswertung 


einer B-Spline-Fläche. 


1 void CubicNubsSurface :: getPointFromUV (const UVPoint& uv, 


Oo AN DHT FF Dä MM ra TU AN Oh OF QUN 


N N N 
mo Pr © 


23 


24 


25 


26 
27 
28 
29 
30 
31 
32 
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HPoint& retPoi 


nt) const{ 


int spanU = -1; 


int spanV = -1; 
__m256d m256_basisU, m256_basisV; 


double result! 


a] = 0 be 


this ->getSpanUV (uv, spanU, spanV); 


m256_basisU = 
m256_basisV = 


const Matrix& 
const Matrix& 
const Matrix& 


int indV = -1; 
int indU = 1F 


this->getBasisFunsU (spanU, uv.u()); 
this->getBasisFunsV (spanV, uv.v()); 


xMatrix = *c_points_x_ptr; 
yMatrix = *c_points_y_ptr; 
zMatrix = *c_points_z_ptr; 


indU = spanU - 3; 
indV = spanV - 3; 
double resultUnpacked[4] = { {0} }; 


__m256d N3j = 
_m256d N2j = 


__m256d N1j = 


m256d NOj = 


_mm256_permute4x64_pd (m256_basisU 


MM SHUFFLE(0, 0, 0, 0)); 


_mm256_permute4x64_pd (m256_basisU 


MM_SHUFFLE(1, 1, 1, 1)); 


_mm256_permute4x64_pd (m256_basisU 


MM_SHUFFLE(2, 2, 2, 2)); 


_mm256_permute4x64_pd (m256_basisU 


MM SHUFFLE(3, 3, 3, 3)); 


N3j = _mm256_mul_pd(N3j, m256_basisV) ; 


N2j 


_mm256_mul_pd(N2j, m256_basisV) ; 


NU = _mm256_mul_pd(N1j, m256_basisV) ; 
NOj = _mm256_mul_pd(N0j, m256_basisV) ; 


// Calculate X 
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__m256d N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd (&( 
xMatrix [spanU-3][indV]))); 
_m256d N2jP = _mm256_mul_pd(N2j, _mm256_loadu_pd(&( 
xMatrix [spanU-2][indV]))); 
__m256d N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd (&( 
xMatrix [spanU-1][indV]))); 
__m256d NOjP = _mm256_mul_pd(N0j, _mm256_loadu_pd (&( 


xMatrix [spanU-0][indV]))); 


__m256d resultX = _mm256_add_pd(_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 
resultX = _mm256_hadd_pd(resultX, resultX); 


_mm256_storeu_pd(resultUnpacked, resultX); 
result[0] = resultUnpacked[0] + resultUnpacked [2]; 


// Calculate Y 

N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd(&(y Matrix [spanU 
-3][indV]))); 
N2jP = _mm256_mul_pd(N2j, _mm256_loadu_pd(&(y Matrix [spanU 
-2][indV]))); 
N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd(&(yMatrix [spanU 
-1][indV]))); 
NOjP = _mm256_mul_pd(NO0j, _mm256_loadu_pd(&(y Matrix [spanU 

-0][indV]))); 


resultX = _mm256_add_pd(_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 
resultX = _mm256_hadd_pd(resultX , resultX); 


_mm256_storeu_pd(resultUnpacked, resultX) ; 
result[1] = resultUnpacked[0] + resultUnpacked [2]; 


Hi eegenen Z 

N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd (&(zMatrix [spanU 
-3]lindV ]))); 

N2jP = _mm256_mul_pd(N2j, _mm256_loadu_pd(&(zMatrix [spanU 
-2]lindV]))); 
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59 N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd (&(zMatrix [spanU 
-1][indV]))); 

60 NOjP = _mm256_mul_pd(NO0j, _mm256_loadu_pd (&(zMatrix [spanU 
-O][indV]))); 

61 

62 resultX = _mm256_add_pd(_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 

63 resultX = _mm256_hadd_pd(resultX, resultX) ; 

64 

65 _mm256_storeu_pd(resultUnpacked, resultX) ; 

66 result[2] = resultUnpacked[0] + resultUnpacked [2]; 

67 

68 retPoint.x() = result [0]; 

69 retPoint.y() = result[1]; 

70 retPoint.z() = result [2]; 

71 } 


Programmauflistung 6.7: Bisher verwendete C++ Implementierung der Gl. (2.22) zur 
Auswertung einer B-Spline-Fläche und Berechnung der ersten Ableitungen. 


1 void NubsConventionalDiss :: getPointFromUV ( double uParameter, 
2 double vParameter, double& pointX, double& pointY, 

3 double& pointZ, double& uTangentX, double& uTangentY, 

4 double& uTangentZ, double& vTangentX, double& vTangentY, 
5 double& vTangentZ ) 

6 { 

7 

8 int uSpan, vSpan; 

9 UVPoint uv (uParameter , vParameter) ; 

0 age 29 

1 int il Ze 

2 

3 getSpanUV (uv, uSpan, vSpan); 

4 compUDers(uSpan, uParameter, 1, *m_pDerivatU); 

5 compVDers(vSpan, vParameter, 1, *m_pDerivatV); 

6 

7 i2 = vSpan - (ndv-1); 

8 for (s=0; s<ndv; s++) { 

9 (*m_pMatrixTmp)[s][0] = 0.0; 
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(*m_pMatrixTmp)[s][1] = 0.0; 
(*m_pMatrixTmp)[s][2] = 0.0; 


il = uSpan - (ndu-1); 
for (r=0; r<ndu; r++) { 
(*m_pMatrixTmp)[s][0] += (*m_pDerivatU)[O][r] * 
(*c_points_x_ptr) [il ][i2]; 
(*m_pMatrixTmp)[s][1] += (*m_pDerivatU) [O][r] * 
(*c_points_y_ptr)[il][i2]; 
(*m_pMatrixTmp)[s][2] += (*m_pDerivatU) [O][r] * 
(*c_points_z_ptr)[il][i2]; 
il++; 
} 


i2++; 


pointX = 0.0; pointY = 0.0; pointZ = 0.0; 


for (s=0; s<ndv; s++) { 


pointX += (*m_pDerivatV)[0][s] * (*m_pMatrixTmp) [s ][0]; 
pointY += (*m_pDerivatV)[0][s] * (*m_pMatrixTmp) [s ][1]; 
pointZ += (*m_pDerivatV)[0][s] * (*m_pMatrixTmp) [s ][2]; 


} 


vTangentX = 0.0; vTangentY = 0.0; vTangentZ = 0.0; 


if (ndv >1 ) { 
for (s=0; s<ndv; s++) { 
vTangentX += (*m_pDerivatV)[1][s] * (*m_pMatrixTmp) [s 


[0]; 

vTangentY += (*m_pDerivatV)[1][s] * (*m_pMatrixTmp) [s 
111]; 

vTangentZ += (*m_pDerivatV)[1][s] * (*m_pMatrixTmp) [s 
DEI 


uTangentX = 0.0; uTangentY = 0.0; uTangentZ = 0.0; 
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56 227 (gh z 234 


57 i2 = vSpan - (ndv-1); 

58 for (s=0; s<ndv; s++) { 

59 (*m_pMatrixTmp)[s][0] = 0.0; 

60 (*m_pMatrixTmp)[s][1] = 0.0; 

61 (*m_pMatrixTmp)[s][2] = 0.0; 

62 

63 il = uSpan - (ndu-1); 

64 for (r=0; r<ndu; r++) { 

65 (*m_pMatrixTmp)[s][0] += (*m_pDerivatU)[1][r] * 

66 (*c_points_x_ptr)[il ][i2 ]; 

67 (*m_pMatrixTmp)[s][1] += (*m_pDerivatU)[1][r] * 

68 (*c_points_y_ptr)[il ][i2 ]; 

69 (*m_pMatrixTmp) [s][2] += (*m_pDerivatU)[1][r] * 

70 (*c_points_z_ptr)[il ][i2 ]; 

71 il++; 

72 } 

73 i2++; 

74 } 

75 for (s=0; s<ndv; s++) { 

76 uTangentX += (*m_pDerivatV)[O][s] * (*m_pMatrixTmp) [s 
110]; 

77 uTangentY += (*m_pDerivatV)[O][s] * (*m_pMatrixTmp) [s 
[1]; 

78 uTangentZ += (*m_pDerivatV)[O][s] * (*m_pMatrixTmp) [s 
112]; 

79 } 

80 } 

81 } 


Programmauflistung 6.8: C++ Implementierung des Algorithmus 2.5|zur Auswertung 
einer B-Spline-Fläche und Berechnung der ersten Ableitungen. 


1 void CubicNubsSurface :: getPointFromUV (double uParameter, 

2 double vParameter, double& pointX, double& pointY, 

3 double& pointZ, double& uTangentX, double& uTangentY, 

4 double& uTangentZ, double& vTangentX, double& vTangentY, 
5 double& vTangentZ) 

6 


{ 
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int spanU = -1; 


int span = -]; 

__m256d m256_basisU, m256_basisV, m256_basisDerivativU , 
m256_basisDerivativV ; 

double resultUnpacked[4] = { {0} }; 


this ->getSpanUV(UVPoint(uParameter, vParameter), spanU, 


spanV ) ; 


m256_basisU = this->getBasisDerivatU (spanU, uParameter, 
m256_basisDerivativU) ; 
m256_basisV = this->getBasisDerivatV (spanV, vParameter, 
m256_basisDerivativV) ; 


const Matrix& xMatrix = *c_points_x_ptr; 


const Matrix& yMatrix = *c_points_y_ptr; 


const Matrix& zMatrix = *c_points_z_ptr; 


int indV = -1; 
indV = spanV - 3; 


|| ========================================= 
// Point evaluation 
// ===== =====4==2 4422442444444 
m256d N3j = _mm256_permute4x64_pd (m256_basisU 
(0, 0, O, (UI: 
m256d N2j = _mm256_permute4x64_pd (m256_basisU 
(il, i, i, Wye 
m256d N1j = _mm256_permute4x64_pd (m256_basisU 
(2, Be Bo nu: 
_m256d NOj = _mm256_permute4x64_pd (m256_basisU 
(dr 3B, Bp Se 


N3j = _mm256_mul_pd(N3j, m256_basisV) ; 
N2j = _mm256_mul_pd(N2j, m256_basisV) ; 


Za Zi 
D 
ou 


_mm256_mul_pd(N1j, m256_basisV) ; 
_mm256_mul_pd(NO0j, m256_basisV) ; 


// Calculate X 
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__m256d N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd(&(xMatrix|[ 


spanU-3][indV]))); 


__m256d N2jP = _mm256_mul_pd(N2j, _mm256_loadu_pd(&(xMatrix [ 


spanU-2][indV]))); 


__m256d N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd(&(xMatrix [ 


spanU-1][indV]))); 


__m256d N0jP = _mm256_mul_pd(NOj, _mm256_loadu_pd(&(xMatrix [ 


spanU-O][indV]))); 


__m256d result_m256 = _mm256_add_pd(_mm256_add_pd(N3jP, N2jP 


), _mm256_add_pd(N1jP, NOjP)); 
result_m256 = _mm256_hadd_pd(result_m256, result_m256) ; 


_mm256_storeu_pd(resultUnpacked, result_m256) ; 
pointX = resultUnpacked[0] + resultUnpacked [2]; 


M @alleullatesx 
N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd (&(yMatrix [spanU 
-3][indV]))); 
N2jP = _mm256_mul_pd(N2j, _mm256_loadu_pd (&(yMatrix [spanU 
-2][indV]))); 
N1jP = _mm256_mul_pd (N1j, _mm256_loadu_pd (&(yMatrix [spanU 
-1][indV]))); 
NOjP = _mm256_mul_pd(N0j, _mm256_loadu_pd (&(yMatrix [spanU 

-0][indV]))); 


result_m256 = _mm256_add_pd (_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 
result_m256 = _mm256_hadd_pd(result_m256, result_m256); 


_mm256_storeu_pd(resultUnpacked, result_m256); 
pointY = resultUnpacked [0] + resultUnpacked [2]; 


// Calculate Z 

N3jP = _mm256_mul_pd (N3j, _mm256_loadu_pd (&(zMatrix [spanU 
-3]lindV]))); 

N2jP = _mm256_mul_pd(N2j, _mm256_loadu_pd (&(zMatrix [spanU 
-2]lindV]))); 
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63 


64 


65 
66 


67 
68 
69 
70 
71 
72 
73 
74 
75 


76 


77 


78 


79 


80 


81 


82 


83 


84 


85 


86 


87 


88 


89 
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N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd(&(zMatrix [spanU 
-1][indV]))); 

N0jP = _mm256_mul_pd(NOj, _mm256_loadu_pd (&(zMatrix [spanU 
-0]lindV]))); 


result_m256 = _mm256_add_pd (_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 
result_m256 = _mm256_hadd_pd(result_m256, result_m256); 


_mm256_storeu_pd (resultUnpacked, result_m256); 
pointZ = resultUnpacked [O0] + resultUnpacked [2]; 


NEE 


|| ========================================= 
Nä = _mm256_permute4x64_pd (m256_basisDerivativU 
MM SHUFFLE(0, 0, 0, 0)); 
N2j = _mm256_permute4x64_pd(m256_basisDerivativU , 
MM SHUFFLE(1, 1, 1, 1)); 
NI = _mm256_permute4x64_pd (m256_basisDerivativU 
MM SHUFFLE(2, 2, 2, 2)); 


NU = _mm256_permute4x64_pd(m256_basisDerivativU , 
MM_SHUFFLE(3, 3, 3, 3)); 


N3j = _mm256_mul_pd(N3j, m256_basisV) ; 
N2j = _mm256_mul_pd(N2j, m256_basisV) ; 
NI = _mm256_mul_pd(N1j, m256_basisV) ; 
NOj = _mm256_mul_pd(N0j, m256_basisV) ; 


// Calculate X 
N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd(&(xMatrix [spanU 
-3]lindV]))); 
N2jP = _mm256_mul_pd (N2j, _mm256_loadu_pd (&(xMatrix [spanU 
-2][indV]))); 
N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd (&(xMatrix [spanU 
-1][indV]))); 
N0jP = _mm256_mul_pd(N0j, _mm256_loadu_pd(&(xMatrix [spanU 
-0]lindV]))); 
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result_m256 = _mm256_add_pd(_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 
result_m256 = _mm256_hadd_pd(result_m256, result_m256); 


_mm256_storeu_pd(resultUnpacked, result_m256); 
uTangentX = resultUnpacked[0] + resultUnpacked [2]; 


A4 Calerie N 

N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd(&(yMatrix [spanU 
-3][indV]))); 
N2jP = _mm256_mul_pd(N2j, _mm256_loadu_pd(&(yMatrix [spanU 
-2][indV]))); 
N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd (&(yMatrix [spanU 
-1][indV]))); 
NOjP = _mm256_mul_pd(N0j, _mm256_loadu_pd (&(yMatrix [spanU 

-0][indV]))); 


result_m256 = _mm256_add_pd (_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 
result_m256 = _mm256_hadd_pd(result_m256, result_m256); 


_mm256_storeu_pd(resultUnpacked, result_m256) ; 
uTangentY = resultUnpacked[0] + resultUnpacked [2]; 


A4 (Caillevilenre 74 

N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd (&(zMatrix [spanU 
-3]lindV]))); 
N2jP = _mm256_mul_pd(N2j, _mm256_loadu_pd (&(zMatrix [spanU 
-2]lindV]))); 
N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd (&(zMatrix [spanU 
-1]lindV]))); 
N0jP = _mm256_mul_pd(NOj, _mm256_loadu_pd (&(zMatrix [spanU 

-0]lindV]))); 


result_m256 = _mm256_add_pd(_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 
result_m256 = _mm256_hadd_pd(result_m256, result_m256); 


_mm256_storeu_pd(resultUnpacked, result_m256); 
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uTangentZ = resultUnpacked[0] + resultUnpacked [2]; 


Aessen 
// vTangent Evaluation 

|| ======================================== 

N3j _mm256_permute4x64_pd (m256_basisU, MM SHUFFLE(0, 0, 


0, 0)); 
= _mm256_permute4x64_pd (m256_basisU , 
il, Ws 

NI = _mm256_permute4x64_pd(m256_basisU , 


MM_SHUFFLE(1, 1, 


MM_SHUFFLE(2, 2, 


Dy Ann: 

NU = _mm256_permute4x64_pd (m256_basisU, _MM SHUFFLE(3, 3, 
By 8) 

N3j = _mm256_mul_pd(N3j, m256_basisDerivativV) ; 

N2j = _mm256_mul_pd(N2j, m256_basisDerivativV) ; 

NI = _mm256_mul_pd(N1j, m256_basisDerivativV) ; 

NOj = _mm256_mul_pd(N0j, m256_basisDerivativV) ; 


// Calculate X 
N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd (&(xMatrix [spanU 
-3]lindV]))); 
N2jP = _mm256_mul_pd (N2j, _mm256_loadu_pd (&(xMatrix [spanU 
-2][indV]))); 
N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd(&(xMatrix [spanU 
-1][indV]))); 
N0jP = _mm256_mul_pd(NOj, _mm256_loadu_pd (&(xMatrix [spanU 
-0]lindV]))); 


result_m256 = _mm256_add_pd (_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 
result_m256 = _mm256_hadd_pd(result_m256, result_m256) ; 


_mm256_storeu_pd(resultUnpacked, result_m256) ; 
vTangentX = resultUnpacked[0] + resultUnpacked [2]; 


// Calculate Y 


N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd(&(y Matrix [spanU 
-3][indV]))); 
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153 
154 
155 
156 
157 
158 
159 


160 


161 


162 


163 
164 


165 
166 
167 
168 
169 } 
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N2jP = _mm256_mul_pd(N2j, _mm256_loadu_pd (&(yMatrix [spanU 
-2][indV]))); 

N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd (&(yMatrix [spanU 
-1][indV]))); 

NOjP = _mm256_mul_pd(N0j, _mm256_loadu_pd (&(yMatrix [spanU 
-0][indV]))); 


result_m256 = _mm256_add_pd (_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 
result_m256 = _mm256_hadd_pd(result_m256, result_m256); 


_mm256_storeu_pd(resultUnpacked, result_m256); 
vTangentY = resultUnpacked[0] + resultUnpacked [2]; 


// Calculate Z 

N3jP = _mm256_mul_pd(N3j, _mm256_loadu_pd (&(zMatrix [spanU 
-3]lindV]))); 
N2jP = _mm256_mul_pd(N2j, _mm256_loadu_pd (&(zMatrix [spanU 
-2][lindV]))); 
N1jP = _mm256_mul_pd(N1j, _mm256_loadu_pd (&(zMatrix [spanU 
-1][indv]))); 
N0jP = _mm256_mul_pd(NOj, _mm256_loadu_pd (&(zMatrix [spanU 

-0]lindV]))); 


result_m256 = _mm256_add_pd(_mm256_add_pd(N3jP, N2jP), 
_mm256_add_pd(N1jP, NOjP)); 
result_m256 = _mm256_hadd_pd(result_m256, result_m256); 


_mm256_storeu_pd(resultUnpacked, result_m256) ; 
vTangentZ = resultUnpacked[0] + resultUnpacked [2]; 
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